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CHAPTER I 
INTRODUCTION 

It is remarkable that interplanetary space provides simul- 
taneously one of the best vacuums available to man, a marvelous labora- 
tory for the study of collisionless plasma, and phenomena that can be 
represented by a continuum model. The reason for the latter is that 
the solar wind, a tenuous high-speed outflow of fully ionized and 
magnetized gas from the Sun, pervades interplanetary space and behaves 
in many ways in conformity with the equations of ordinary .continuum 
niagnetohydrodynamics and gasdynamics. This is true with respect to 
both large-scale irregularities in the solar wind itself and to the 
interaction of the solar wind with planetary or lunar-sized objects. 

This thesis will be concerned primarily with the development of an 
appropriate theory based on fluid concepts that accurately represents 
the large-scale features of the interaction between the solar wind and 
the planets Earth, Mars, and Venus and that advances the present 
understanding of such phenomena. The basic approach of the presenta- 
tion followed throughout will be first to review and discuss observed 
phenomena that are pertinent to this topic , and then to formulate and 
develop a mathematical model that will adequately represent the ob- 
servations . 

Extensive theoretical studies, specific calculations, and 
detailed comparisons with observations have shown that the interac- 
tion of the solar wind and the Earth may be represented satisfactorily, 
insofar as gross features of the flow such as density, velocity, tem- 
perature, and magnetic field are concerned, by the standard continuum 
equations of magnetohydrodynamics (see Spreiter £t al. [1968] or 
Spreiter and Alksne [1969, 1970] for recent summaries). This thesis 
extends the same line of analysis to objects in the solar system 
which have no significant intrinsic magnetic field, but a sufficiently 
dense and electrically conducting ionosphere to stop and deflect the 
solar wind before it is absorbed by the planetary surface or atmosphere. 
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Mars and Venus are known examples of’ such objects. It is the plan of 

the present investigation to base the mathematical analysis of these 
< 

interaction phenomena with nonmagnetic planets on the equations of 
magnetohydrodynamics representing the steady flow of a perfect gas 
having infinite electrical conductivity as well as both zero viscos- 
ity and thermal conductivity and then to proceed to an approximate 
solution through simplifications that are well known and tested in 
previous applications to the interaction of the solar wind with the 
Earth. Modifications are introduced only as required because it is 
the ionosphere rather than the geomagnetic field through which the 
pressure of the planet is communicated to the solar wind. 

In broad outline, the elements of the theory are as follows. 
The incident solar wind is considered to be steady and supersonic, and 
to flow in accordance with the equations of magnetohydrodynamics for 
a perfect gas in which all dissipative processes, such as those asso- 
ciated with shock waves and boundary layers, take place in the inter- 
ior of thin layers idealized as magnetohydrodynamic discontinuity 
surfaces of zero thickness. The ionosphere, or at least the outer 
part of it that participates in the interaction with the solar wind, 
is idealized as spherically symmetric and hydrostatically supported 
plasma having infinite electrical conductivity. Since the two bodies 
of plasma are of different origin, and have different properties, 
they must be considered to be mutually impenetrable in the idealized 
hydromagnetic representation, and to be separated by a tangential 
discontinuity surface. This surface is called the ionopause, since 
it marks the outer boundary of the ionosphere. Once the problem is 
formulated, the solution proceeds in two steps, just as in the analo- 
gous problem for the Earth. First the shape and location of the iono- 
pause is calculated for selected values of the parameters characteriz- 
ing the solar wind and the ionosphere. Follov7ing that, the location 
of the how wave and the properties of the flow field are determined 
from the equations of gasdynamics to which those of magnetohydrody- 
namics have been simplified. It is found. that the results for a wide 
range of ionospheric pai'ameters may be brought into close correspond- 
ence with those for solar -wind flow past the Earth's magnetosphere by 
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application of a simple geometric transformation of the coordinates. 

Finally, these theoretical results based on gasdynamic theory are 

« 

compared with the observations made by Mariner 5 as it flew by Venus 
and Mariners 4, 6, and 7 as they flew past Mars, and implications of 
the points of agreement and disagreement are discussed, 

A more refined analysis of this interaction of the solar wind 
with the Earth, Mars, and Venus in v?hich the complete equations of 
magnetohydrodynamics are solved without any simplification is then 
carried out for the case of the magnetic field aligned with the flow 
direction. Furthermore, -for Mars and Venus the exact- shape and loca- 
tion of the ionopause is calculated using the hydromagnetic theory 
without introducing any approximations. Solutions of this improved 
theory are computed, and the location of the bow wave and the proper- 
ties of the flow field are determined. For high Alfven Mach number 
this refined hydromagnetic solution, shows that the bow-wave location 
and the flow-field properties are very similar to those of the 
simpler gasdynamic theory, but for low Alfven Mach num.bers they are 
strikingly different. This unusual bow-wave position for flow past 
the Earth has recently been confirmed by spacecraft measurements. A 
detailed comparison of these two models is carried out to. determine 
the accuracy and usefulness of the simpler gasdynamic theory. Lastly, 
the results of the hydromagnetic solution are sho’S'jn to be in very good 
accordance with the observations made by Mariner 5 as it flew past 
Venus although certain differences near the theoretical location of 
the ionopause suggest the presence of a thick boundary layer . 

The present version of the gasdynamic model for the interac- 
tion with Mars and Venus is based on a preliminary analysis originally 
reported by Spreiter, Summers, and Rizzi (1970b). In addition, be- 
cause of the similarity of the flow past Mars and Venus with that 
past the Earth as indicated. by the correspondence rule and also be- 
cause our more accurate hydromagnetic solution has direct application 
to the interaction with the Earth, we choose, for the sake of unity, 
to parallel our gasdynamic analysis for the interaction v?ith Mars and 
Venus with the gasdynamic investigations already carried out in the 
past years for the case of the Earth. For this purpose we borrow 
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heavily from the work of Spreiter ^ (1966a, 1968) and Spreiter 

and Alksne (1969, 1970). 

However, before embarking on the major topic of this thesis, 
let us reviev7 the basic properties and idealizations of the solar 
wind itself in order to provide a background knowledge of the inter- 
planetary medium and also to illustrate the basic coherence of the 
entire theoretical description of the solar wind and its interaction 
with the planets. Much of this is repeated from Spreiter et al . 
(1968) . 
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CHAPTER II 
THE SOLAR WIND 

Although geomagnetic evidence (Chapman and Ferraro, 1931) 
had long before led to the belief that ionized gas was projected all 
the way from the Sun to the Earth and beyond, both occasionally from 
flares and more steadily in localized 'M-region' beams, the modern 
discussion of the solar plasma begins with Biermann's (1951, 1957) 
analysis of the behavior of ionized comet tails. An outstanding 
photograph of a more recent comet, Ikeya-Seki, taken at the Mojave 



Figure 2-1. Comet Ikeya-Seki, October 29, 1965 
(Spreiter ^ , 1968) . 


head of the comet near the horizon at the bottom of the photograph 
and the extremely elongated tail. Although it had become customary 
to ascribe the antisolar direction of cometary tails to radiation 
pressure of sunlight, Biermann shov?ed this explanation to be inade- 
quate in many respects, and proposed in its place the concept of ar 
approximately radial high-speed flow of ionized gas from all parts 
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the Sun at all times. It is not difficult to accept this explanation 
of a photograph like that of figure 2-1, since the general effect is 
very remindful of smoke from a smudge pot or of a wind sock in a 
strong and steady wind. 

Although the solar corona is not normally visible because of 

the overwhelming brilliance of the light from the photosphere, it is 

clearly visible, even to the unaided eye, at the time of a total 

eclipse of the Sun by the Moon. As a result, much of what is kncxm 

about the solar corona has been learned from observations made during 

the few-minutes duration of such events. In more recent years, such 

observations, particularly of the parts closer to the Sun, have been 

augmented by those taken at other times with the aid of a coronagraph, 

* 

a telescope system fitted with special obscuring equipment invented 
by Lyot in 1937. Figure 2-2 shows a photograph of the corona taken 
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Figure 2-2. Solar corona during eclipse of May 30, 1965, 
from NASA GV-990 aircraft at 38,000 feet 
(Spreiter ^ , 1968) . 


during the eclipse of May 30, 1965, from a NASA Convair 990 jet air- 
liner at 38,000 feet altitude (Smith ^ 1967). Clearly visible 
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beyond the disk of the Moon are the polar rays, similar in shape to 
the magnetic lines of force at the poles of a magnetized sphere. 

More interesting from the present point of view is the large outward 
extension of the corona from the equatorial and middle latitudes. 
Although these plumes appear to extend only a few solar radii in this 
photograph, other photographs taken during other eclipses have sho’t-jn 
rather similar features visible to distances as much as five or ten 
times farther from the Sun. Observations indicate that the coronal 
gas in these plumes is both somewhat hotter and many times denser 
than over the poles , 

Beginning effectively with the flight of Mariner 2 to Venus 
in 1962, direct measurements betoeen 0.7 and 1.5 astronomical units 
(a.u.) have demonstrated beyond all doubt the existence of a continuous 
flux of magnetized plasma from the Sun, essentially as predicted theo- 
retically by Parker (1958) and termed by him the solar wind. However 
even this detailed information from the spacecrafts characteristically 
possesses substantially greater uncertainties than would usually be 
experienced in laboratory measurements of the same properties of a 
plasma. In particular, the rather low resolution of directional, 
temporal, and energy aspects of the plasma measurements resulting from 
weight, space, power, and telemetry limitations imposed by the space- 
craft tends to limit the quality of the determination of the velocity 
distribution of the particles and ultimately, therefore, the know- 
ledge of the bulk properties, such as density, velocity, and tempera- 
ture, of the plasma. Further uncertainties arise from the non- 
Maxwellian and anisotropic character of the velocity distribution of 
the random motion of the particles. Nevertheless, a good understand- 
ing of the broad features of the properties of the solar wind is now 
in hand, and further refinements can be anticipated in the near future. 

As a review, some 'of the basic properties and idealizations 
of the solar wind pertinent to our study are summarized in figure 2-3, 
which was originally presented by Spreiter et a^. (1968). First of all 
it should be recognized that it is the velocity of the solar wind and 
not the planet's motion in its orbit around the Sun, that provides the 
dominant contribution to the relative motion of the planet and the 
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interplanetary plasma. In many ways we may think of the solar wind 
as a hypersonic gale that blows through interplanetary space. From 
an overall point of view, the solar wind may be considered to result 
from forces in the interior of the Sun which are too large to be con- 
strained by the feeble pressure that can be imposed by the tenuous 
cool gas of the interstellar medium, ks the immediate source of the 
solar wind, we may look to the unsteady conditions in the relatively 
cool visible surface, or photosphere, of the Sun. Here are to be 
found a wide spectrum of acoustical, hydromagnetic , and internal 
•gravity waves which grow -tremendously in amplitude as they propogate- 
upward into an atmosphere of diminishing density. The dissipation of 
these waves, either directly or indirectly after their culmination in 
shockwaves, results in a substantial deposition of energy generally 
believed to be sufficient to account for the high temperature of the 
solar corona relative to the photosphere. Since the energy is being 
deposited -in the form of heat at a greater rate than can be transported, 
away by radiation or thermal conduction, convection in the form of the 
solar wind ensues. This process is believed to account for both the 
2 X 10^ °K temperature of the corona, compared with the 5800 °K 
temperature of the photosphere, and the outward convection that leads 
ultimately to the supersonic solar wind at greater distances from the 
Sun. The average flow starts gradually with insignificant velocities 
of the order of 1 km/sec or less in the lower corona, and increases 
steadily with increasing distance from the Sun until its ultimate 
velocity is attained beyond a few tens of solar radii. The expanding 
coronal gas carries with it the solar magnetic field, stretching the 
lines of force outward through the solar system and enhancing the 
strength of the interplanetary magnetic field there v/hile generally 
.maintaining direct connection to the Sun. Because of the rotation of 
the Sun, the magnetic lines of force wind in spirals in the plane of 
the ecliptic making, on the average, an angle of about 45° with the 
direction from the Sun at the orbit of the Earth. 

Representative values for the observed properties of the solar 
wind' in the vicinity of the Earth's orbit are outlined in the box in 
the center of figure 2-3. The bulk velocity V of the solar wind has 
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Figure 2-3. The solar -wind, a hypersonic gale in 

interplanetary space (Spreiter ^ al , 
1968) . 


been observed to vary between about 275 and 850 cm/sec , zeith the loxeer 
values usually associated with the quiet solar wind and the higher 
values xeith disturbed conditions in space. The number density fl 
of the protons in the solar wind, which is substantially more diffi- 
cult to determine from the plasma-probe data than V j is usually 
found to be between 1 and 10 cm“3 in the quiet solar wind and as much 
as an order of magnitude greater during disturbed times. The intensity 
of the magnetic field B is usually betxjeen about 3 and 7 ^ { I'S ~ 
i0~5 gauss), but may also be substantially greater during disturbed 
times. All of these quantities have been monitored with varying de- 
grees of exactness and completeness for several years, and are knox-m 
to display substantial variations with time and location. The chemical 
composition is principally ionized hydrogen with a small mixture, 
usually about 4 percent, but sometimes more than 20 percent, of helium 
nuclei (llirshberg et al., 1970). 

At great distances from the Sun, the solar x^ind must merge 
somehox'7 into the interstellar medium. Although the spatial distribu- 
tion of interstellar gas, and its degree of ionization, is both very 
non-uniform and uncertainly determined, temperatures and densities of 
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the order of those indicated in figure 2-3' are frequently quoted (see, 
e.g. Kaplan, [1966]) . It is generally presumed that the region of super- 
sonic solar-x^ind flow terminates with a shock wave and that the flow 
beyond is subsonic. The location of this terminal shock wave is usually 
estimated to be between 50 and 100 a.u. , although Brandt (1962) has 
-suggested, on the basis of an analysis of the directions of comet tails, 
that it may be as near as 2 a.u. (see Betton ^ [1963] and Brandt 

[1967] for later comments). One aspect of this terminal shockwave 
that does not appear to have been adequately considered in these dis- 
cussions is that the temperature attained by the gas behind the terminal 
shock wave is far in excess of the temperatures quoted above for the 
interstellar medium. Whether the gas cools gradually as it continues 
to flovj subsonically away from the Sun by means of some process such 
as radiation not normally included in the analysis of the solar wind 
or remains at high temperature throughout a spherical shell terminated 
at still greater heliocentric distances by a second shock advancing 
through the interstellar medium or behaves in some other manner alto- 
gether remains essentially unkno™. It may be, on the other hand, that 
the supersonic solar wind terminates without the intervention of a 
shock wave in the manner suggested by Faus (1966) through the action • 
of non-adiabatic transfer of kinetic energy to the interstellar medium 
in a roughly constant -pres sure process that takes place between about 
1 and 3.5 a.u. Corresponding uncertainties plague any discussion of 
the interaction of the solar wind with’ objects beyond the orbit of Mars. 

Ultimately, at greater distances from the Sun, the subsonic 
solar wind must interact in some way with the interstellar gas . Al- 
though non -uniform! tie s’ in the spatial distribution and uncertainties 
in the measurements are great, fluid considerations suggest the develop- 
ment of an interface, the heliopause, separating the solar and inter- 
stellar plasma. However, so little is knoxm about the nature of the 
interaction and the properties of the gases involved that even the 
existence of the heliopause must be regarded as hypothetical. 

In spite of the obvious deviations from spherical symmetry in 
the solar corona, and presumably in the solar wind as well, the theory 
of the solar wind is generally based for reasons of mathematical 
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tr.actability on the assumption of radially symmetric flow. This 
theory was initially presented by Parker ( 1958 ) and summarized in 
more developed form by Parker ( 1963 , 1965 ), Liist ( 1963 , 1965 ), 

Scarf ( 1964 , 1970 ), Dessler ( 1967 ), Hundhausen ( 1968 , 1970 ), and others 
Perhaps the most fundamental concept is that of treating the solar 
wind, even as it becomes extremely rarefied far beyond the orbit of 
the Earth, as an ordinary continuum gas. Thus to a theoretician, a 
very significant finding of plasma and magnetic-field measurements in 
space is the degree to which not only the broad features of the solar- 
wind flow, but also much of the small-scale features, can be repre- 
sented by the continuum equations of fluid dynamics . Although it is 
apparent that not all matters of significance can be treated in this 
way (for example the anisotropy of the velocity distribution), the 
amount of understanding that is possible is truly astonishing, par- 
ticularly in view of the fact that the solar-x-7ind plasma is one of the 
most rarefied gases nov? accessible to direct measurement. Hox^ever, 
it X'jas alx-jays knoxra that the fluid description could not be justified 
in terms of ordinary Coulomb-collision effects alone. In his earliest 
papers Parker pointed out that as the gas streams axjay from the Sun 
into the supersonic region, the collisional mean free path must rapidly 
increase so that near the Earth it is on the order of 1 a.u. Neverthe- 
less, in these initial discussions Parker also provided a general jus- 
tification for the use "of fluid equations beyond the dense collisional 
region. He pointed out that some microscopic plasma instabilities 
must naturally develop and that associated xjave-particle scattering 
then produces a thermalization mechanism leading to a finite mean 
free path and effective fluid transport characteristics. It is nox\r 
quite clear that plasma turbulence generates partial thermalization at 
1 a.u. Apparently the microscopic scattering mechanism produces ran- 
domizing and isotropizing effects so that the gasdynamic equations for 
conservation of mass and momentum remain appropriate for descriptions 
of the main flox-/ properties, at least out to 1 or 2 a.u. Hoxjever, 
microscopic interactions may play a significant role in energy trans- 
port and energy balance (see Scarf [ 1970 ] for a reviex^ of microscopic 
phenomena in the solar wind) . Clearly these phenomena are important 



in the study of the microscopic structure, but the equations of 
fluid dynamics appear to be appropriate for the description of the 
macroscopic features of the solar wind and its interaction with the 
planets Earth, Mars, and Venus. 



CHAPTER III 


BROAD FEATURES OF THE OBSERVATION! 

Data acquired in recent years in the vicinity of the Earth, 

Moon, Mars, and Venus have ’disclosed three essentially different types 
of interaction with the solar wind as illustrated in figure 3-1, For 
the Earth (Ness ^ 1966; Spreiter et ^l. , 1968; Binsack, 1968), the 

geomagnetic field prohibits the solar wind from approaching nearer than 
about 10 Earth radii under ordinary circumstances. The shielded region, 
the magnetosphere, acts as an obstacle in the supersonic flow of the 
solar wind, and a bow wave forms a few Earth radii upstream of the mag- 
netosphere boundary, or magnetopause. Neither Venus nor Mars has a 
significant magnetic field, but the electrical conductivity and density 
of the upper ionosphere of both Mars and Venus are sufficiently great 
that the solar wind is prevented from flox^ing directly into either the 
planetary surface or the lower absorbing levels of the atmosphere 
(Fjeldbo and Eshleman, 1968, 1969). The solar wind is thus deflected 
around the ionosphere, and a boxj wave is formed upstream of the planet, 
similar in many ways to that associated x-?ith the Earth. Aside from 
evident differences in the underlying physical processes at the boundary 
surface betx^een the ionosphere and the solar x^ind, the ionopause, the 
principal difference betxJeen the flow fields around Mars or Venus and 
Earth is the size of the cavity. As illustrated in figure 3-1, the 
ionopause is wrapped much closer around Mars and Venus than the mag- 
netopause is around the Earth, the nose being at an altitude of about 
500 km for Venus and probably about 155 ’to 175 km for Mars, compared 
x-7ith about 60,000 km for the Earth. 

Interaction of the solar wind and the Moon is very different 
because of the lack of either a sensible magnetic field or ionosphere 
tq deflect the incident flox^ (Lyon et al. , 1967; Ness et , 1967; 
Spreiter et , 1970a). The particles of the solar wind thus proceed 
unchecked until they contact the lunar surface, x^here they are effectively 
removed from the floxxr. As a consequence, no box7 xxrave forms, and the 
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principal features of the flox^ field are associated with the closure 
of the lunar -wake, or cavity in the solar wind (Colburn et al. , 1967), ' 
that extends domstream from the Moon as illustrated. 
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CHAPTER IV 

FORMULATION OF FLUID REPRESENTATION 

Just as in the theory of the solar wind itself, tne funda- 
mental assumption underlying nearly all theoretical analyses of large- 
scale -features of the interaction of the solar V7ind v?ith the planets 
is that the average bulk properties of the flow can be described ade- 
quately by 'the continuum -equations of fluid dynamics. That thrs should 
be so is not at all obvious since the mean free path of the particles, 
based on Coulomb interactions, is of the order of half the distance 
to the Sun at the orbit of the Earth, Indeed, theoretical justifica- 
tion is essentially qualitative at the present time, and usually con- 
sists of referring to randomizing and isotropizing effects of small 
irregularities, much like a -form of plasma turbulence, which are always 
present in the solar wind. The details of these processes are not well 
understood at this time, and it is fair to state that the real support 
for the use of the continuum- fluid model is provided by the outstanding 
agreement between results calculated in this way and those actually 
measured in space. Moreover, the difficulties in solving the nonlinear 
continuum equations and the richness of the phenomena contained within 
this theory are sufficiently great so that it is worthwhile to explore 
the consequences of this model before introducing additional complica- 
tions . 

If an ionized gas like the solar wind moves in a magnetic 
field, electric fields are induced in it and electric' currents flow. 

The magnetic field' exerts forces on these currents which may consider- 
ably modify the floxj. In turn the currents themselves modify the mag- 
netic field. Thus the interaction between the magnetic and fluid 
dynamic phenomena is complex, and the flov? must be examined by combin- 
ing the magnetic-field equations with those of fluid dynamics. Hence 
the appropriate mathematical model capable of description of this phe- 
nomena is provided by the standard equations of magnetohydrodynamics 
(see, e.g. Landau and Lifshitz [i960]). 
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Eor explicit application to the study of the solar wind's 
interaction with planets, these equations can be simplified by treat- 
ing the flox-jing plasma as a steady, calorically perfect, and dissi- 
pationless conducting gas. The wind is knoxm to fluctuate in time, 
but very fex-7 theoretical analyses of these effects have been carried 
out as yet. There is little doubt that such studies must be made if 
many important dynamical properties of the interaction are to be under- 
stood. However, in at least some instances, observations (Spreiter 
and Alksne, 1969; Binsack and Vasyliunas, 1968) and theory (Spreiter 
and Summers, 1965) suggest that many features of the flow are able to 
respond rapidly enough to achieve an essentially steady-state equilib- 
rium V7ithin a fex^ minutes after a sudden change in the solar wind. 
Consequently, the steady-state solution for conditions existing at a 
given time should adequately represent the flox^ properties of the wind 
at that time. Similarly, the assumption that the solar wind is per- 
fectly conducting, inviscid, and nonheat conducting is also not strictly 
valid, but the full dissipative equations are much more complicated, 
and as yet a solution cannot be achieved. Hence, in order to fully 
develop and examine the fluid theory, our model is restricted to the 
ideal equations in xxrhich all dissipative processes, such as those as- 
sociated x^ith shock x-yaves , boundary layers, and current sheets, take 
place in the interior of thin layers idealized, as magnetohydrodynamic 
discontinuity surfaces' of zero thickness. Finally,' even though some 
of its properties are knox-m to be anisotropic, the solar x-7ind is taken 
to be a calorically perfect isotropic gas primarily because its high- 
temperature and lox7 pressure satisfy the conditions for a perfect gas, 
but also because a more detailed theory of the thermodynamic equation 
of state for a txro-species ionized' g'as embedded xjith a magnetic field 
has not been developed yet. The theoretical justification for these 
three assumptions is essentially heuristic in nature, and, like the 
•original pos.tulate of representing the solar X’j'ihd by a fluid model, 
the real support derives -from the fact that extensive theoretical 
studies, specific calculations, and detailed comparisons xji'tK observe- . 
tions have shox7n that the interaction of the solar xxrind and the Earth 
has been represented satisfactorily, insofar as gross features of the 
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flow are concerned, by this model (see Spreiter ^ al. [1968] or 
Spreiter and Alksne [1969] for recent summaries). Moreover, recent 
investigations (see Spreiter et al . [l970a]for an extensive account) 
have shoTim that the same is true for the Moon, provided that boundary 
conditions are applied at the lunar surface that correspond to the 
Moon absorbing or neutralizing all particles of the solar wind on 
impact . 
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CHAPTER V 


FUEDAMEKTAL EQUATIONS 

The following fundamental differential equations of magneto - 
hydrodynamics for the steady, flow of an isotropic and calorically 
perfect gas having infinite electrical conductivity as well as both 
zero viscosity and thermal conductivity (Spreiter ^ al. , 1966a) 


V'PV = 0 

P(V*V)V 4 - V p “ ~ ^ xcuriB 

ro pi/ ‘ ! i ^ rv; /nj 

Curl (B xV )=0 divB =0 (^-d 

rNJ rvJ 

(V-V)S = 0 S - So = Cv I'/) -T/,%, 

~ ~ (P/Po> 


provide the basis for all theoretical matters to be discussed in thi 

thesis. In these and the equations to follow, the S 3 mibols yO j P 

S 5 and y refer to the density, pressure, entropy, and' velocity 

of the gas; and B refers to the magnetic field expressed in terms 
1 ^ ' ■ 

gaussian units. Cp and are constants representing the specif 

heats at constant pressure and constant volume. The quantity = 
Cp/Cv is the ratio of specific heats and is equal toN/(N+2) 


The use of this standard notation contrasts with the practice 
in the much older science of geomagnetism of using the symbol H to 
represent the Earth’s magnetic intensity. .Discussion of this property 
does appear in some space-research papers but causes no confusion 
because the magnetic permeability of the solar wind is very 'nearly ■ 
unity and hence in gaussian units B = H . These two quantities will 
be used interchangeably in this thesis in order to avoid making any 
alteration of original sources. . . 
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where N is the number of degrees of freedom of the gas particles. 

On first thought the choice of the value 2 for the ratio of 
specific heats may seem justified in a variety of ways, generally 
centering around the presumed two-degrees-of-freedom nature of the 
interactions of charged particles in a magnetic field. This argument 
weakens, however, when consideration is given to the extremely irreg- 
ular character of the magnetic fields observed in space, particularly 
do\mstream of the bow shock wave. In fact, the supposition of apply- 
ing hydromagnetic theory to the flow of solar plasma around the planets 
involves the assumption of an isotropic pressure, and the whole con- 
cept appears more internally consistent if the particles are considered 
to behave as if they have three rather than two degrees of freedom. 

Thus for a monatomic gas, such as either of the two leading constitu- 
ents of the solar wind, atomic hydrogen and helium, N = 3 and = 
5/3 . The subscript o refers to conditions at an arbitrary refer- 
ence location in the flow field. 

Important auxiliary relations for the temperature T , speed 
of sound d , internal energy Q , and enthalpy h are as follows: 


p = pRT/m = nkT 

d =V9P/0P = '\ITPIp — V/RT/n^ (5-2) 

e = CvT h = CpT = e + p/p 

where R = { Cp,— CY)m = 8.314 x 10^ erg/®K is the universal gas 
constant, and iTi ,is the mean molecular weight nondimen sionali zed so 
that m = 16 for atomic oxygen. For fully ionized hydrogen plasma, 
rn is thus 1/2 , and H = number of particies/cm^ = 2Hp where ri[> is 
the number of protons /cm^ , and k = Boltzmann’s constant = 1.38 x 10“^^ 
erg / °K . The relationship between 0 and H is easily seen to be 
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p — rnj> Hji where trip = 1.67 x 10“^^ gm is the mass of a proton. 

It is also fully equivalent and convenient for some purposes to re- 
place the entropy equation (V •■V)-S 0 in equations (5-1) by the 

following energy equation: 



(5-3) 


Implicit in the relations of equations (5-1) are two impor- 
tant nondimensional parameters that characterize the flow at any point, 
the Mach number M and the Alfven Mach number defined by 

where c3 represents the speed of sound defined in equations' (5-2), 

A the speed of an Alfven wave which will be discussed later, and 
V — |V| represents the speed of the gas. 

Although only first derivatives appear in the system of 
differential equations (5-1) , the neglected dissipative terms are 
described by second derivatives. Their neglect requires that the 
gradients be small. In-magnetohydrodynamics , as in gasdynamics, how- 
ever, compressions tend to coalesce and steepen into shock waves of 
such small thickness that they can be considered virtual discontinui- 
ties for many purposes. In this vray gradients tend to become very 
large , and continuous motions tend to break dovm some place if they 
Involve compressions. In addition, attraction between like electric 
currents tends to cause distributed currents to collapse into thin 
sheaths, across which the magnetic field changes nearly discontinuously 
Mathematically, continuous solutions of the dissipationless equations 
cease to exist beyond the point of breakdown, and the flow is no longer 
governed solely by the differential equations (5-1). Mass, momentum, 
magnetic flux, and energy must still be conserved, however, and con- 
servation relations provided by integrals of equations (5-1) and (5-3) 


Ma 


V 


V 




(5-4) 
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must hold between quantities on the two sides of any such discon- 
tinuity. The conservation equations for magnetohydrodynamic discon- 
tinuity surfaces are (see, e.g. Landau and Lifshitz [i960]) 


PVn 


0 



(5-3) 


PVn(h+vV2)‘+ VnB]/47T -Br,y-'B/47T 



Unit vectors normal and tangential to the discontinuity surface are 
designated by n and t , and components of y and B in these 
directions are indicated by subscripts H and t . The square 
brackets are used to indicate the difference between the enclosed 
quantities on the two sides of the discontinuity, as in [q] = Q,- Qq 
where subscripts O and 1 refer to conditions on the upstream and 
dowstream sides of the discontinuity. These relations are frequently 
supplemented by the statement that there is a current sheet flowing 
along the discontinuity surface and that the value O of this cur- 
rent per unit width is given by 






curl B 


(5-6) 


It is now well knovm (Landau and Lifshitz, 1960; Jeffrey and Taniuti, 
1964; Jeffrey, 1966) that five classes of discontinuities are described 
by the conservation equations (5-5) , but for the sake of completeness 
and unity we repeat some of the mathematical development establishing 
the basis for the existence of these different classes. The particular 
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version presentea cioseiy toiiows tnat ot ypreiter et (.iybba) . 

Friedrichs and Kranzer (1958) originally developed an alterna- 
tive form for the conservation equations that affords a quick survey 
of the possible types of discontinuities and some of their properties, 
and is particularly useful in the determination of asymptotic proper- 
ties of weak discontinuities. I-t is based on extensive use of mean 
values = (Qp+ Q, )/Z 5 and the new variables W = \ jp and 

hn = representing the specific volume and the mass flux across the 

discontinuity. In terms of these variables, the relations given in 
equations (5-5) become 


m[w]-[v] = 0 


m 


m 


m 




4ff 

<b)M - bJvI = 0 


(5-7) 


[b„] = 0 


The last of these can be replaced by 


- 0 


m 


e + (p) W 






NN= 


0 


(5-8) 


If the mean quantities are considered known, there are thus eight 
equations in eight scalar variables [y/]‘ / [p] , , [b] , 

and l^c] since the condition Bj-| = 0 always holds. 

Since 'e appears in only the last relation of equations (5-7), 
v?e can consider alternatively the seven scalar equations in seven 
variables defined by the first three relations of .equations ( 5 - 7 ). 

They are all linear homogeneous equations, and solutions exist only if 
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their determinant vanishes. Equating the determinant to zero yields 
the following equation for the mass flux m 


<W)m(w)iTf- 



(5-9) 


"With the density p^' upstream of the discontinuity given, this equa- 
tion can also be considered as an equation for the normal component 
of the incident velocity • 

Immediately we see that tn = O is a solution to equation 
(5-9) and that three additional plus-minus pairs of real roots exist 
for which tn ^ O . To follov? Landau and Lifshitz (1960) and Spreiter 
et al . (1966a) , discontinuities that lie along streamlines ( m = 

— O) are called tangential discontinuities or contact discontinuities 
according to whether or not the normal component of the magnetic field 
Bp) vanishes. Discontinuities across which there is flow ( tfi ^ O ) 
are divided into categories. Those associated with the roots On — 

± <CV\C>]'^” are called rotational discontinuities, although they 

are frequently termed transverse or intermediate shock waves by many 
authors. The term shock wave is reserved here, however, for the dis- 
continuities associated with the four remaining roots. 


A. Tangential Discontinuities 

Tangential discontinuities are defined as those in which both 
tn and Bf) vanish. For these conditions, equations (5-5) and (5-7) 
yield the following relations between the quantities on the two sides 
of the discontinuity surface 


Vn=Bn=0 [Mt]=0 [gt]=.0 


[P] ^ 0 


p + b/stT 


= 0 


(5-10) 
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Although the velocity and magnetic field are required to be parallel 
to a tangential discontinuity, these relations shw^ that the density 
and the tangential components of the velocity and magnetic field .may 
have jumps of any magnitude. The other thermodynamic quantities, such 
as temperature, entropy etc., are also discontinuous in accordance 
■with their definitions provided in' equations (5-2). The sum of fluid 
pressure p and magnetic pressure ©Tt must , however , be con- 

tinuous across a tangential discontinm'tv. 

B. Contact Discontinuities 

The properties of contact discontinuities are given by solv- 
ing equations (5-5) and (5-7) with hi = O and 0 . They are 

as follows:' 


Vn^O Bn^O [v] =[B]=[p] = 0 ■ 0 (5-n) 


These relations show that the fluid not only flo-^^s parallel 
to the contact discontinuity in a steady flow, but that the velocity, 
pressure, and magnetic field must be continuous. The density, and 
therefore the temperature, entropy, and other thermodynamic variables, 
may have any discontinuity however. 

C. Rotational Discontinuities 

The properties of rotational discontinuities are determined by 
solving equations (5-5) or (5-7) with m = • 

They are 




1/2 


[yj - [Bj/(4jt^) 


1/2 


(5-12) 


[P] = [P] 


V, 


[B„] = [B^] = 0 


These relations show that only the tangential component of the velocity 
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and the direction of the magnetic -field vector change across a rota- 
tional discontinuity. The flow-velocity component Vr, normal to a 
rotational discontinuity is equal to the phase velocity of an Alf\ren , 
wave in the direction of the wave normal, and is moreover independent 
of the strength of the discontinuity. There is thus no tendency for a • 
series of weak rotational discontinuities in a homogeneous medium 
either to overtake one another and coalesce into a shock wave or to 
spread out as in an expansion fan. Although there has been considerable 
discussion of possible effects of traveling rotational or Alfven x^aves 
in the magnetosphere and in space (see, e.g. MacDonald [1964] for a 
review), and it has been suggested (Levy ^ 1964) that the mag- ■ 

netosphere boundary may be resolved into a rotational discontinuity 
and a slox 7 expansion fan, the role and importance of stationary rota- 
tional discontinuities in the steady-state flox^ of solar plasma past 
the planets has not been established. 


D. Shock Waves 

The mass flux tTl through magnetohydrodynamic shock waves 
satisfies the equation 


(W)m^ 




[p]Byif7f[w] = 0 


(5-13) 


obtained by setting to zero the last factor of equation (5-9). 

Spreiter £t ad. (1966a) shov?ed that equations (5-7), together with the 
further requirement stemming from the entropy considerations that all 
shock X'javes are compression x^aves, yields the folloxvfing relations be- 
txireen the quantities on the txTO sides of such a shock wave 




2 (B>-[§; 


[B„] = 0 



m[w]Bh<gt)/4'?r 

(VV)m^ ~ BV4T( 
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M=-N 


m 


[w]<0 


[p] > 0 


(5-14) 


where In represents any of the four roots of equation (5-13). The 

first of these, they pointed out, indicates that the sum 

and difference [B-j.] of the tangential components of the magnetic 

fields on the two sides of the shock wave lie in either the same or 

/% 

opposite directions. This requires the shock normal H and the 
magnetic -field vectors on the two sides of the discontinuity to he 
coplanar. The velocity relations also show that the tangential com- 
ponent of the velocity discontinuity Vf — Vf “ Vf^ is in the 
same or opposite directions as and therefore in the same 

plane as Bj and Bt> 

Spreiter £t (1966a) discussed the properties of shock 
waves characterized by the various roots for Tfi by rewriting equation 
(5-13) to read 


m 4 - [pj/[vv||(w)m^-Bp/4"^^ ~ 


(5-15) 


Since the roots iTl ~ fn^' ~ were considered before , 

it follows that , and conversely 

that tris < -[p] / [w] . if . shock xvaves character- 

ized by the larger value flip for the mass flux are called fast, and 
those by the smaller value slow. The mass flux across either 

class of shock wave must therefore satisfy the inequalities 


> mt 


:> ms :> -[p]/[W] > m^ (5-16) 


It follows immediately from equations (5-14) and (5-16) that 
and increase through a fast shock wave and decrease through a 

slow shock wave. 

Solutions of equations (5-5) or (5-7) are determined by 
straightforward but lengthy algebraic manipulations, details of which 
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are available in several books (see, e.g. Anderson [1963]' or Jeffrey 
and Taniuti [1964] for extensive accounts). It is important to recog- 
nize, however, that these equations possess extraneous solutions that 
cannot occur in nature. These solutions were originally referred to 
as unstable, but this term has now been generally abandoned since it 
has become known that this is not an ordinary instability growing 
exponentially with time, but rather a sudden disintegration of the 
shock wave. Alternatively, as noted by Spreiter (1966a), 

a solution containing an extraneous shock wave does not have neigh- 
boring solutions corresponding to a small change in the boundary con- 
dition which requires an arbitrarily small angle or rotation of the 
plane containing the magnetic field. Identification of the physically 
relevant solutions cannot be made on the basis of entropy considera- 
tions alone, as in ordinary gasdynamics. It is necessary in' magneto - 
hydrodynamics to consider how a shock wave could evolve through waves 
of small amplitude overtaking one another and coalescing. The appro- 
priate requirements for a physically relevant solution can be stated 
in general mathematical terms by a pair of evolutionary conditions (see 
e.g. Jeffrey and Taniuti [1964]). A more physical description that 
leads to the same conclusions has been given by Kantrowitz and Petschek 
(1964), From either point of view, it is found that no additional con- 
ditions need be imposed for fast shock waves and that the conserva- 
tion equations (5-5) possess extraneous solutions only- for slow shock 
waves. Fortunately, the extraneous solutions can be easily recognized 
by the fact that they indicate the tangential component of the mag- 
netic field to be directed oppositely on the two sides of the. discon- 
tinuity, a physical impossibility in all but certain degenerate cases. 
Since the first of equations (5-14) shows that the tangential component 
of the magnetic field cannot reverse direction through a fast shock 
wave, one has the general' result that the tangential components of the 
magnetic field on the two sides of any physically relevant shock wave 
must always lie in the same direction. Spreiter (1966a) pointed 

out the further consequence that the flow velocity must be greater than 
the rotational x^ave speed on both sides of a fast shock xxrave and less 
than the rotational v?ave speed on both sides of a sloxv shock wave. 
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Thus a fast shock wave overtakes a slow shock wave, and conversely a 
slow shock wave cannot overtake a fast shock wave . In application 
to the supersonic flow of the solar wind past either the Earth or Mars 
or Venus, we therefore anticipate that the bow wave would be a fast 
magnetohydrod 3 naaniic shock wave. However, slor? shock waves, as well as 
additional fast shock saves, can be expected in the flow downstream 
of the bow V7ave. 


E. Weak Shock Waves 

The relations given in equations (5-13) and (5-14) simplify 
for weak shock waves to the point where they can be expressed explicitly 
in terms of the sound speed c3 , the Alfven speed A , and its com- 
ponents Af\ and ‘Af normal and tangential to the plane of the dis- 
continuity surface. The latter are defined as 


A = -\/iyi^ 


An = A 


COS 9 


Ah- 




A 


s I n 9 


(5-17) 


where 0 refers to the angle between the magnetic field B and the 
shock normal D . For such shock waves, the jump [q] in every 
quantity Q is much smaller than the quantity itself, and the desired 
simplification can be accomplished by simply removing the mean value 
brackets and replacing, for clarity of expression, the difference 
brackets with tf , i.e. ^9^ — Q and [oj — cJ Q • Equation (5-13) 
may then be solved to yield the following expression for the normal 
velocity component Vn of the incident stream; 

( \ 


Vn = 


3 ^ + 


1/2 


(5-18) 
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With Vn taken positive, the larger of the two values obtained by 
using the plus sign before the inner radical pertains to the fast 
wave and the smaller to the slow wave. 

Equation (5-18) coincides with the corresponding expression 
for the velocity component normal to the characteristic surfaces of 
equations (5-1) and can be recognized, upon changing to a reference 
frame fixed in the fluid upstream of the .shock wave, as the appro- 
priate expression for the phase velocities of fast and slow traveling 
hydromagnetic plane waves. The propagation of such waves is ’conveni- 
ently described graphically by the familiar normal-speed or Friedrichs 
I diagrams 5 examples of which are shown in figure 5-1 originally pre- 
sented by Spreiter £t al. (1966a) for ci = A , c3 = A » 
cj hj^ ♦ In these diagrams is plotted as a function of 

angle 6 between the shock normal and the magnetic field H (here 
taken parallel to the horizontal axis) . The outer curve in each dia- 
gram represents the results for the fast wave, and the inner curves 
the results for the slow wave. 



Figure 5-1. Normal speed or Friedrichs I 

diagrams for propagation velocity of 
weak plane waves (Spreiter _et _al. , 1966a) 

Also included are dotted curves representing the normal velocity of a 
rotational discontinuity defined by equation (5-12) . For ease of 
illustration, the diagrams v?ere drawn for the plane containing the 
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•wave no3nnal and the magnetic- field vectors. The corresponding results 
for more general orientations are represented by the corresponding 
three-dimensional surfaces formed by rotating the plane curves of 
figure 5-1 about the 6=0 axis. These plots illustrate that, for 
any direction 0 , the speed of a rotational wave is intei'mediate 

between that of the fast and slow waves. Further, Vn for the fast 
wave is equal to the larger of dl and A when the wave normal is 
parallel to H (0 = 0), and to (c3^+ when . the wave normal 

is perpendicular to H (0 = -?-) • For the slow wave, Vn vanishes 
when 8= 7f/2 5 and is equal to the smaller of c3 and A when 
0=0. If A approaches zero vyhile c3 remains finite, the sur- 
faces representing the slow and rotational waves disappear and the 
surface representing the fast wave approaches a sphere of radius c3 . 
Under these circumstances , the equations of magnetohydrodynamics ap- 
proach those of gasdynamics, and calculations based on aerodynamic 
methods should provide a good approximation for most applications. If, 
on the other hand , c3 approaches zero while A remains finite , . the 
surfaces representing the slow wave disappear while the one represent- 
ing the fast wave approaches a sphere of radius A , and those repre- 
senting the rotational discontinuity remain spheres of radius A/ 2 
as always. Although equation (5-18) clearly shows that the normal 
velocity of fast and slow waves are invariant with interchange of 
values for c3 and A and that the former become independent of di- 
rection 0 as c3 vanishes, the rotational discontinuities retain a 
.finite normal velocity, and neither the differential equations nor the 
shock relations of magnetohydrodynamics approach those of gasdynamics 
in this limit (except in a special and restricted sense for flows in 
which the velocity and magnetic -fie Id vectors are parallel at all 
points [Kanti*owitz and Petschek, 1964]). 

Spreiter ^ (1966a) indicated that the jump relations 

given in .equations (5-14) similarly reduce in the limit of weak shock 


waves to 
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where the quantity in the last bracket of the expression for 
gives the direction and sign of that quantity. The corresponding ex- 
pressions for the change in pressure p and entropy S across a 
weak shock wave are (see, e.g. Jeffrey and Taniuti [1964]) 
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The latter relation shows that the change in entropy through a weak 
shock wave is proportional to only the third power of 6p j p , 
and hence it is vanishingly small for weak waves. Physically, this 
indicates that infinitesimal expansions, as well as compressions, can 
be considered without vio.lation of the entropy requirement. Thus 
can be either positive or negative. Moreover, equation (5-20) shows 
that although the change in pressure through a weak discontinuity is 
proportional to (S'Pf P the first term in which -the magnetic field 
appears is proportional to • {cfP/pY • The changes in entropy and 
pressure through magnetohydrodynamic expansions or compressions are , 
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thus the same as in. gasdynamics until the strength of the discontin- 
uity is sufficiently great that third -order terms must be retained. 

As indicated by equation (5-19), however, this statement does not ex- 
tend to other quantities such as the velocity or the magnetic field. 

These results for weak compression and expansion waves are 
useful for describing conditions at great distances from the Earth or 
other disturbing obstacles where it- can be safely assumed that 
<< 1 . They are not sufficient for the discussion of the entire 

bow -wave problem, however, because typical conditions in the solar 
wind close by the Earth are such that ^p/p easily exceed unity 

near the nose. The maximum value for the density ratio across a 


hydromagnetic shock is finite, however, and given simply by 
(5 + I just as in gasdynamics. 
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CHAPTER VI 

APPLICATION TO PLANETARY INTERACTION 

The concepts and equations of magnetohydrodynamics which 
were reviewed in the preceding chapter can now be applied to the 
specific study of the interaction of the steady solar wind and the 
planets Earth, Mars, and Venus, 

Extensive theoretical investigations, specific calculations, 
and detailed comparisons with observations have been carried out in 
the past by many workers (see Spreiter ^ ail. [1968] and Spreiter and 
Alksne [1969, 1970] for recent summaries) in -the study of the gross 
features of the interaction of the solar wind with the Earth. We 
have recently extended and modified the hydromagnetic theory of the 
solar-wind flow past the Earth so that it is applicable to nonmag- 
netic planets, such as Mars and Venus, that have a sufficient iono- 
sphere to deflect the solar plasma around the planet and its atmos- 
phere. (A preliminary version of this work was reported by Spreiter 
^ aJ.. [1970b]). Here, because of the similarity of these two theo- 
retical analyses, we choose next to review those elements of the 
theory of the interaction with the Earth which are necessary for com- 
parison with and completeness of this more comprehensive account of 
solar -wind interaction with nonmagnetic planets. 

Of the five classes of discontinuities described in the last 
chapter, only the tangential discontinuity has properties compatible 
with those found in space (Spreiter and Alksne, 1969, 1970) describ- 
ing the boundary surface that separates the planet from the flowing 
solar plasma. Specifically, for the Earth the condition Bh = 0 
complies with the finding that there is no connectivity between the 
geomagnetic and interplanetary fields (Behannon and Fairfield, 1969). 
For planets such as Mars and Venus that have no significant intrinsic 
-magnetic field, but a sufficiently dense and electrically conducting 
ionosphere to stop and deflect the solar wind before it is absorbed 
by the planetary surface or atmosphere, B-, •= 0 states that at the 
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boundary the interplanetary field is parallel to it. In both cases, 
the condition Vh = 0 prohibits flow across the surface, which thus 
can be identified as a streamline and therefore a boundary of the flow 
Arbitrary differences in the density and the velocity and magnetic - 
field components parallel to the boundary are allowed, but the sum of 
the gas pressure p and magnetic pressure must be the same 

on the two sides of this surface. 

The properties of contact discontinuities appear to be appro- 
priate for the boundary of the distant wake far downstream from the 
planets. There, any current system or velocity discontinuity, such 
as exists on the for\^ard part of the Earth’s magnetopause, has ample 
time to diffuse or decay no matter hoxj slight the departures from 
perfect electrical conductivity or inviscid flow. Under these cir- 
cumstances the 'external magnetic field penetrates the discontinuity 
and ultimately at great distances doxmstream may be expected to tra- 
verse the x\?ake completely v7ithout distortion (Spreiter et , 1966a) . 
Thermodynamic properties of the wake remain different from those of 
the surrounding flox^, hox^ever , because of the different previous his- 
tories of the gases. 

As pointed out before, too important parameters that charac- 
terize the flow at any point are the Mach number M and the Alfven 
Mach number defined by equations (5-4) . The former is the ratio 

of the flox7 velocity to the speed of sound c3 , and. the latter is 
the corresponding ratio of velocity to the speed A of a rotational 
or Alfven x^ave propagating along the direction of the magnetic field. 
Values for c3 and A for conditions typical of those encountered in 
the solar X’yind as it flox^s past Venus, Earth, and Mars are illustrated 
in figure 6-1. The solar wind is knoxm to vary substantially v?ith 
time., but number densities of the order of 2,5 to 25 protons /cm^, 
magnetic fields of 3 to 10 if , and temperatures of 50,000 to 100,000 
‘'K may be considered representative. Since the velocity of the in- 
cident solar x<7ind ranges from about 300 to 800 km/sec, it is evident 
that the free-stream Mach number Mco and Alfven Mach number 
are generally much greater than unity. And for normal conditions, the 
interaction of the solar x^ind and the planets is related to highly 
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Figure 6-1. Speed of aound c3 and of Alfven waves A 
in fully ionized hydrogen for conditions 
representative of those in the solar wind 
(Spreiter and Alksne, 1970). 

supersonic and; superAlfvenic maghetohydrodynamic flows in which M 
and are of the order of 10. 

Even though it is possible for weak disturbances to travel as 
fast as { c3^-h in a direction normal to the magnetic field, 

the. maximum propagation speed is much less than the flow speed. Con- 
sequently, it can be inferred that a standing shock vjave must exist 
in the flow around these planets. Of the five classes of discon- 
tinuities outlined above, Spreiter ^ (1968) pointed out that 

only the fast shock x^ave can be used to represent this standing box? 
wave because the solar x\’ind approaches these planets x^ith' a mass -flu? 
Pco^co that greatly exceeds that of a rotational discontinuity, anc 
therefore also of a slox^ shock xxrave. 


A. Asymptotic Directions of Shock Waves and Wakes 

The Friedrichs X diagrams of figure 5-1 describe the propa- 
gation speeds of x^eak plane x^aves as vlex^ed in a coordinate system ir 
xvhich the undisturbed plasma is at rest. In applications to steady 
floxjs, the corresponding diagram for the infinitesimal standing wave 
pattern of a point disturbance is often of greater interest. The 
directions of such x-zaves are of significance for finite disturbances 
as x7ell because they coincide x-jith the as^miptotic directions of the 
standing wave pattern at- great distances from the obstacle. • So when 
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}V5j^ and given together with the directions of the velocity 

and magnetic-field vectors Vqo and Bro in iite undisturbed incident 

r~* 

solar wind, the asymptotic directions of the planetary tail (or wake) 
and -shock waves at great distances from an obstacle can be determined 
immediately by a simple geometrical construction based on these 
norraal-wave-speed diagrams and originally devised by Spreiter e_t al . 
(1966a) . The results apply equally whether the obstacle is in effect 
the confined geomagnetic field as in the case of the Earth's magneto- 
sphere, or the highly conducting ionosphere of either Mars or Venus. 

This construction is illustrated in figure 6-2(a) for a case 




Figure 6-2. Sketches showing how asymptotic directions of 
shock waves, rotational discontinuities, and 
wakes can be determined by use of Friedrichs I 
or II diagrams. (a) Using Friedrichs I diagram 
for plane waves, (b) Using Friedrichs II diagram 
for waves from a point disturbance (Spreiter 
et al . , 1966a) . 

in which the magnetic-field vector ^cd is inclined 45“ from the 
free-stream-veloclty vector , and Meo . Mco i® 

chosen to be 2 , abnormally small values for the solar wind, in 
order to better show the construction details. For and 

of the order of 10 , the Vcb vector would be much longer, and all 
waves would make much smaller angles with the direction of flow. For . 
ease of illustration, only the traces of the various three-dimensional 
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surtaces in the plane containing the velocity and lield vectors are 
‘shown. In this plot the flow is to the right and parallel to the 
horizontal axis. The obstacle is supposed stationary at the center 
of the Friedrichs diagram which has been rotated by the 45~degree 
angle between ,§oo • With the center at the tip of the 

vector V CO and 9 measured from the direction of , equation 

(5-18) is used to construct curves representing the normal velocities 
of weak shock waves and rotational discontinuities . Except for a 
translation of the origin and rotation of the axis, the latter are of 
course the same as shown in figure 5-1 for cS / A = 'JZ . The asymp- 
totic directions of the various discontinuities at great distances 
from the obstacle can now be determined by application of the theorem 
that states that the two lines drawn from any point on a circle to the 
extremities of a diameter intersect at a right angle. The required 
construction is thus performed simply by drawing a circle of radius 
V 00/2 with center at and connecting the origin of Vco to 

the various points of intersection with the normal speed curves by 
straight lines as indicated in figure 6-2(a). These lines are not 
tangent to the Friedrichs I diagram, as in the simpler isotropic case 
of gasdynamics, because this diagram represents the speed of plane 
waves and is very different from the Friedrichs II diagram for the 
speed of waves from a point source. The asymptotic directions can 
also be constructed by drawing straight lines from the origin of the 
vector Vco tangent to the Friedrichs II diagram as illustrated in 
figure 6-2(b). This procedure often requires more effort because the 
expression describing the Friedrichs II diagram is much more complex 
than equation (5-18), The construction employed in figure 6-2(a) is 
more direct and fully equivalent for the present purposes. 

Also indicated in figure 6-2 is the aS37mptotic direction of 
the wake or tail, which as ‘noted previously would be represented by 
either a tangential discontinuity or a contact discontinuity depend- 
ing upon whether or not the normal component of the magnetic field 
vanishes. In either case Spreiter ^ al. (1966a) noted that equations 
(5-10) and (5-11) state that Vd ~ 0 indicating that the vzake.must 
be aligned parallel to Vco independently of the direction of the 

f>J 
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magnetic field. The orientations of the other discontinuities illus- 
trated in figure 6-2 change, however, as the direction of the magnetic- 
field vector is altered. The extent of these changes can be readily 
visualized, moxover, since a change in the direction of the magnetic 
field leads to a rotation of the Friedrichs speed diagrams about the 
end of the fixed vector , and this in turn leads to associated 

changes in the intersection points of figure 6 -2(a) or the tangency 
points of figure 6-2(b). 

B. Relation Between Hydromagnetic and Gasdynamic Flows 

Although the diagrams of figure 6-2 are drawn for the case 
in vjhich the free-stream Mach number is 2 and the speed of sound is 
times the Alfven speed, that is for (c^/a)^ re 
= '/2 , the qualitative character of these diagrams remains the same 

for all Mach numbers and ratios greater than unity. If 

Aco should become small relative to , however, as might readily 

occur if the magnetic field should diminish in intensity, the inner 
loops representing the propagation speeds of rotational and slow waves . 
become small relative to the outer ovals representing the fast waves. 
Finally, if the magnetic field approaches zero, the inner loops shrink 
toward a point at the end of the vector while the outer oval 

approaches a circle (or sphere in the corresponding three-dimensional 
representation) of radius . In this way Spreiter ^ (1966a) 

indicated that the fast hydromagnetic wave degenerates to the Mach 
wave of ordinary gasdynamics, and the rotational and slow waves lose 
their physical significance. The fluid and electromagnetic aspects 
of the flow thus decouple, and the fluid motion is described entirely 
by the equations of gasdynamics. The approach to the limiting case 
appears, moreover, to be free of singular behavior, and gasdynamic 
theory should thus provide a useful approximation to hydromagnetic 
flows when (cS/a)^^ is substantially greater than unity. In such 
approximations the gasdynamic Mach number should be associated with ■ 
the Mach number M of the corresponding hydromagnetic flow, rather 
than the Alfven Mach number as was done at times in the past. 

.If on the other hand the Alfven speed is substantially greater 
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than .the sound " speed , as would occur if the magnetic field is suffi- 
ciently strong, the inner loops representing the slow wave speed again 
shrink toward a point at the end of the vector, while the outer 

oval representing the fast wave speed approaches a circle of radius 
Aq 3 . The resulting asymptotic directions for fast waves in a hydro- 
magnetic flow of Alfven Mach number thus approach those of 

gasdynamics for Mach number Mco ~ ^"^Aco * described by Sears 

(1960) and others, however, the flows about a given obstacle may be 

/ 

quite different. For certain combinations of , and 

angles between Vco and the physically significant waves are 

those that extend upstream rather than doi;TOStream from the disturbance. 
When such conditions prevail, hydromagnetic flow about a given obstacle 
may tend to resemble gasdynamic flow about the same obstacle, but with 
the flow direction of the related incident stream reversed. For other 
combinations of these quantities, hydromagnetic flow about a given ob- 
stacle bears little, if any, resemblance to its gasdynamic counterpart. 

C . Shape of the Boundary 

In order to calculate the shape of the boundary for a given 
planet, the physical criterion that determines which of the three 
classes of interactions sketched in figure 3-1 occurs must be suitably 
expressed so that the jump conditions for a tangential discontinuity 
can be applied. To stop the solar wind above the subsolar point, as 
in the flows past the Earth, Mars, and Venus, the sum of the planetary 
gas pressure p and the magnetic pressure B’/S'Tf' must be sufficient 
to balance the stagnation pressure of the solar xvind at a height 

great enough so that the ionosphere gas can be considered to be ef- 
fectively collisionless. For Earth, all estimates and observations 
of the magnitudes of the gas pressure p and the magnetic pressure 
BVeiT in the outer magnetosphere lead to the conclusion that B^/ 8"?T 
p , x^hereas for nonmagnetic planets like Mars and Venus, 
spacecraft measurements have shown that 8^ « p . Hox<rever, for 
the Moon the solar wind vias found to flox-/ directly into the surface 
which indicates -that 'P + 'BVS'^ • Intermediate cases are easily 
conceived, but observations have not yet indicated a need for their 
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consideration. 

]. , Location of the magnetopause 

For the Earth, then, exceeds p throughout all 

the magnetosphere above a few hundred kilometers. So far as the flow 
outside the magnetosphere is concerned, the magnetosphere boundary 
may thus be represented by the limiting case of a tangential discon- 
tinuity in which there is a vacuum (p= 0^ P=0) oiie side, and 
across which the sum of the gas and magnetic pressure of the exterior 
•flow must be balanced entirely by the magnetic pressure of the magnet- 
osphere (Spreitex ^ al. , 1966a) . With the neglect of small effects 
of electrical currents in the ionosphere, the dominant features of the 
magnetic phenomena that the impinging solar wind produces in the outer 
magnetosphere are represented by the terms that remain in equations 
(5-1) when p and p are equated to zero, namely 

ciiv B = 0 curl B=0 (6-i) 

The form and size of the magnetosphere boundary are therefore effec- 
tively independent of the details of the flow within the magnetosphere 
and- determined primarily by the interaction between the magnetic 
field within the outer magnetosphere and the hydromagnetic flow around 
the exterior. There is still flow to be expected within the magneto- 
sphere, but it would be required to move about as if it were confined 
in a container described by the magnetopause . 

To complete the analysis, a mathematical representation for 
the intrinsic magnetic field of the Earth is required. This is provided 
with sufficient accuracy for the present purposes by a magnetic dipole 
at the center of the Earth having a strength such that |^©j ~ Bec^ 

= 0.312 gauss at the geomagnetic equator, and oriented so that the 
north geomagnetic pole is at 78.6^ North latitude and 70.1° West longi- 
tude (Chapman and Bartels, 1940), Its extension into space according 
to equations (6-1) is described by 
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B©= -Beq(ae/r)^(9sm 0 + r2cos 0) ( 6 - 2 ) 

8 

in which c3e = 6.37 x 10 cm is the radius of the Earth, P is the 

geocentric distance, and © is the polar angle measured with respect 

A A 

to the north geomagnetic pole, and P and 0 are unit vectors in the 
P and 0 directions. Thus, the mathematical problem of determining 
the shape of the magnetopause' and the distortion of the geomagnetic 
field confined therein is now completely reduced to finding the solu- 
tion to equations (6-1) that has a dipole singularity given by 
equation (6-2) at the origin and that simultaneously matches with 
the solution for the exterior flow through the tangential discon- 
tinuity conditions Bn = 0 and “ (Pd-B^87() gxt 

the unknOTm location of the magnetopause. 

The resulting free-boundary problem is, however, extremely 
complex because the location and shape of the magnetosphere boundary 
is not known a priori , but must be determined as part of the solution 
much as in the classical theory of free- streamline flows in hydro- 
dynamics. The problem as such is too difficult to solve, and reason- 
able approximations must be sought , At the magnetosphere boundary 
one simplification of great practical utility can be made by introduc- 
ing an approximation to the tangential discontinuity conditions . The 
basis for this is that space observations and theoretical estimates 
indicate that the pressure p of the solar wind tends to be much 
larger than the magnetic pressure bVSIT in the exterior flow around 
the forward part of the magnetopause (Spreiter ^ j 1966a). Con- 
sequently, as far as the flow exterior to the magnetosphere is con- 
cerned, the equations for the discontinuity at the magnetopause can be 
satisfactorily approximated by those of the limiting case of a tan- 
gential discontinuity in x\?hich there is a vacuum on one side and no 
magnetic field on the other side, i.e., 

Br, = 0 ’ (BV8Tf),n-|- = (P)ext 

However, since both and remain unknown and must satisfy 
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equations (6-3) at an unknown location, the problem as such continues 
to be intractable. 

Two additional approximations finally reduce the problem to 
amenable form. The first, used in all existing analyses, is the in- 
troduction of the Nex-?tonian approximation for the pressure on the sur- 
face of an obstacle in a high-Mach-number flow. Thus Pg^-j- 
sumed given by 


Pext = ^ = K^„vJcos2 Y C6-4) 

in which ^ is the angle betx^een the outx^ard normal to the magneto- 
pause and the incident free-stream-velocity vector , subscript 

GO refers to conditions in the solar wind upstream of the bow wave, 
and K is a constant equal to 0.88 for high-Mach-number flow of a 
monatomic gas although usually taken as unity in most applications 
(Spreiter ^ al., 1966a). The Newtonian approximation greatly sim- 
plifies the problem because it decouples the calculation of the shape 
of the magnetosphere boundary from the detailed analysis of the sur- 
rounding floxj. The final approximation assumes that the unknovm mag- 
netic field at the magnetopause, x^hich results from the dis- 

torting of the geomagnetic field by the solar x^ind, is simply equal 
to E f times the tangential component of the Earth's dipole field, 
xvThere f is a constant usually taken to be unity (Spreiter and Briggs, 
1962). Its use, together x^ith the Nex^tonian pressure approximation, 
in the discontinuity conditions (6-3) leads to a partial differential 
equation for the geocentric distance P to the magnetopause. Figure 
6-3 illustrates the solution of that equation for the equatorial plane 
(Beard, 1960), the noon meridian plane (Spreiter and Briggs, 1962), and 
the remainder of the magnetosphere boundary (Briggs and Spreiter, 1963) 
for the case in xdiich the dipole axis is perpendicular to the solar- ‘ 
wind floxj. Higher-order solutions which do not make use of the last 
approximation for B|p|f have been given more recently (Mead and Beard, 
1964; Olson, 1969), but the results differ from each other and from 
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those of the first-order solution by only about the width of the 
gray lines of figure 6-3. As may be seen by comparison with figure 



Figure 6-3. The- boundary of the geomagnetic field 
(Spreiter £t ad.. , 1968). 

3-1, the neutral points indicated on figure 6-3 have physical signifi' 
cance as the points on the magnetopause from which extend the only 
magnetic field lines to connect directly with the Earth. 

2 . Location of the ionopause 

For Mars or Venus the magnetic field is weak or nonexistent 
(B^8'iT«p) , but the gas pressure p of the ionosphere is suffi- 
cient to stop the solar wind well above the planetary surface. So 
in this case, the ionosphere rather than the planetary magnetic field 
communicates the presence of the planet to the solar xvind. The iono- 
sphere, or at least the outer part of it that participates in the in- 
teraction with the flowing solar .plasma, is idealized as spherically 
symmetric and hydrostatically supported plasma having infinite elec- 
trical conductivity. Since the two bodies of plasma are of different 
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origin and have different properties, they must he mutually impene-. 
trable in the idealized hydromagnetic representation and be separated 
by a tangential discontinuity surface. This surface is called the 
ionopause because it marks the outer boundary of the ionosphere 
(Spreiter eh , 1970b). 

More explicitly, the assumption of hydrostatic support is 
equivalent to assuming that all motions of the gas within the iono- 
sphere are sufficiently small with respect to the planetary body 
that equilibrium exists between the pressure gradient and the force 
of gravity, thus 


dp 

dr 



(6-5) 


where p and p are the gas pressure and density, r is the radial 

distance from the center of the planet, and is the acceleration of 

gravity. Values for the latter are inversely proportional to P ^ 

Thus g rr 93 ( Ps/P 3 where subscript S refers to values at 

the surface of the planet. Values for g^ for Venus, Mars, and 

2 

Earth are about 870, 375, and 982 cm/sec , and those for Pg are 

8 

about 6.1, 3.4, and 6.4 x 10 cm, respectively. The pressure is as- 
assumed to be related to the density by the perfect gas lax<7 

p = nkT = ^RT/m . 

which xms originally stated in equations (5r-2). The mean molecular 
mass Ffl is equal to 1/2 for fully ionized hydrogen plasma, and 1 and 2 
for singly ionized molecular hydrogen and helium. The density may be 
eliminated from equation (6-5) by introduction of equation (6-6), and 
the result integrated to yield 



(6-7) 
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which shows that H represents the height interval in which the pres- 
sure decreases by a factor 6 . If the variation of g with P is 

disregarded over the range of application of equation (6-9), T is 
also constant, and the density varies with height in the same way as 
the pressure, thus 


^ ~ P 

'^R Pr. 



( 6 - 10 ) 


In the upper atmosphere, where little mixing would be presumed 
to occur, diffusive equilibrium may be’ considered to "prevail among the 
various ionic constituents, at least in an idealized sense. As em- 
phasized recently by Bauer (1969), however, this equilibrium is more 
complicated in an ionosphere than in a neutral atmosphere for which 
Dalton's law of partial pressures is customarily invoked. This is be- 
cause an electric polarization field, dependent on the mean ionic mass 
and charged -particle temperatures, acts on all the ions and couples 
the diffusive equilibrium distribution of all of them. In most cir- 
cumstances, however, the lightest ionic constituent, ionized atomic 
hydrogen, would emerge as dominant a't great altitudes. At lower and 
intermediate attitudes, more specific knowledge, such as is now, re- 
sulting from space experiments conducted by the USSR and USA, is 
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needed to specify the variation of pressure with radius with any 
degree of precision. 

Figure 6-4 shows the variation with distance from the center 
of Venus of the electron number density, as deduced from the 



igure 6-4. Profiles of electron concentrations in the night 
and day ionospheres of Venus, as deduced from the 
dual-frequency-radio -occultation measurements by 
Fjeldbo and Eshleman (1969) . The dashed lines 
represent the average profile independent of 
strong fluctuations in the data. 
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dual-frequency-radio-occultation measuraments by Fjeldbo and Eshleman 
(1969). For a plasma with singly ionized ions, charge neutrality re- 
quires that these curves also describe the ion number density, summed 
over all constituents. These measurements indicate that the nose of 
the ionopause is at a distance of about 6500 km from the center of 
Venus, or about 450 Ion above the surface. At high altitudes, there 
is some uncertainty in the data associated with the reference level 
of interplanetary ionization so that the nightside profile may be 
somewhere between the two indicated curves. Between 6300 and 6600 km, 
however, the experimenters have reported that the nightside ionosphere 
probably consists of ionized helium at a temperature between 330 and 
710 ‘‘K. And although there is some ambiguity in the data above 6600 km, 
they suggest that at high altitudes the dominant constituent may be 
either ionized atomic or molecular hydrogen x^ith a plasma temperature 
of either 675 or 1350 °K. 

The proper interpretation of these data, and also those from 
the Mariner-5 -Lyman-alpha measurements of Barth ^ (1967), has, 

hov?ever , been the subject of considerable discussion. Barth (1968) 
and Barth ^ al. (1968) conclude that the entire upper atmosphere of 
Venus above about 6500 km can be represented best by an atmosphere con- 
sisting primarily of molecular hydrogen at 650 i 50 °K. On the other 
hand, Wallace (1969) and Donahue (1969) have concluded from their 
analysis of the Mariner -5 -Lyman-alpha data that the dominant constitu- 
ent above 6500 km is deuterium at about 650 McElroy and Strobel 

(1969) have compared results from a number of models of the topside 
nighttime ionosphere of Venus 'xjith observations of Mariner 5 and con- 
cluded that the primary ionized constituent is either helium or molecu- 
lar hydrogen. More recently, Whitten (1970) has performed a similar 
study for a number of models of the upper daytime ionosphere and con- 
cluded that a dominant amount of ionized helium also provides a suit- 
able explanation of the spacecraft's findings. 

In viexj. of these uncertainties, and even greater ones x-jith 
respect to Mars at the time of this x^riting, x^e put aside any further 
discussion of refinements of the variation of H x^ith altitude even 
though there are obvious violations of the conditions for spherical 
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S 3 aninetry 5 constant scale height, and hydrostatic pressure support, 
particularly near the dawn meridian of the planet and at great alti- 
tudes where H cannot possibly be constant because of diminishing 
g . However, one inadequate result of this is that equation (6-9) 
then predicts the implausible occurrence of the planetary pressure be- 
coming vanishingly small at only a distance of several planetary radii 
on the nightside of the planet. This prediction directly conflicts 
with the analysis of the occultation data from Mariner 5 suggesting 
that there may be a wake of planetary ionization extending a great 
distance in the antisolar direction (Fjeldbo and Eshleman, 1969). 

To reconcile this difficulty equation (6-9) is slightly modified so 
that instead of approaching zero, the planetary pressure now varies 
exponentially to a small, finite value p, , that corresponds to the 
very low pressure in the tail of the ionosphere at great distances 
in the antisolar direction. Thus we alter equation (6-9) to read 


P 








( 6 - 11 ) 


- Because nothing about the structure within the tail is knoxjn and 
also, because at very great distances this tail must somehow merge 

* V 

with the solar wind, is taken to be 'a constant equal to 

^ / 8 'TT the sum of the gas pressure and magnetic pressure in 

the undisturbed solar wind. The value' of is about two orders of 

magnitude smaller than the maximum pressure on the ionopause. If 
the constant is taken to be the pressure at the subsolar point 

of the ionopause, where the maximum pr'essure occur's, then the coeffi- 
cient of the exponential term differs- from p by only about 1 per- 
cent. Therefore, the exponential variation of p with T is main- 
tained on the daytime side of the ionosphere, and the small constant 
py^ does not bear effect in equation (6-11) until r reaches a value 
of a few planetary radii and the exponential becomes much smaller 
than p^ . Because the ionopause is located so near the dayside of 
the planet that r can only become this large in the antisolar direc- 
tion, equation (6-11) satisfactorily predicts the approximate constant 
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pressure for the ionized planetary wake, which ultimately 

balances the sum of the magnetic and gas pressure of the solar wind. 
Since equation (6-11) corresponds at least generally with the little 
amount of information so far obtained from spacecraft measurements, 
we adopt equation (6-11) as adequate to represent the variation of 
p with r in the upper atmosphere of either Venus or Mars. 

Results will be presented for several values for H , how- 
ever, in order to provide an indication of the variations to be ex- 
pected betx^een ionospheres of different chemical composition. Ror 
Venus, with Tp equated to 6500 Icm and T to 700 "K, H is approx- 
imately 1500, 760, 380, and 35 km for ionized atomic hydrogen, singly 
ionized molecular hydrogen, helium, and carbon dioxide. The corre- 
sponding values for the nondimen sional ratio 5 in which Pq , 

the distance from the center of the planet to the nose of the iono- 
pause, is equated to 6500 km are 0.23, 0.115, 0.058, and 0.0052. 

For Mars, the smaller radius and acceleration of gravity and lower 
ionospheric temperature of about 200 °K (Kliore 1965; Fjeldbo 

and Eshleman, 1968) lead to values that are about 1.2 times larger 
• than for Venus, namely about 0.28, 0.138, 0.069, and 0.0062. The 
calculations of these values makes use of the result, indicated by 
the data from Mariners 4, 6, and 7 together with the present theory, 
that the nose of the ionopause is at an altitude between 155 and 175 
km, and that P^ therefore scales nearly in proportion to the plane- 
tary radius. 

In a manner exactly analogous to the mathematical formulation 
of the Earth's magnetopause, th'e boundary of the ionosphere is repre- 
sented, under the assumption that the upper ionosphere can be treated 
as a perfectly conducting fluid effectively bound to the planet and 
incapable of mixing with the solar -wind plasma, by a tangential dis- 
continuity. For Mars and Venus the planetary magnetic field is weak 
or nonexistent, and as far as the flow around the ionosphere is con- 
cerned, the ionopause can be modeled by the limiting case of a tangen- 
tial discontinuity in which there is a perfectly conducting iono- 
.sphere ' on the interior side, and across which the sum of the gas and 
magnetic pressure of the exterior flow must be balanced entirely by 
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the pressure of the upper ionosphere. For this situation the tangen- 
tial discontinuity equations become 


Vn= Bn=0 


(P)mt=(P + 


and relate conditions on the interior and exterior sides of the iono- 
pause . Therefore the mathematical problem of calculating the shape 
of the ionopause is now completely reduced to finding the solution 
for the exterior flow that simultaneously matches with equation (6-11) 
through the tangential discontinuity conditions of equations (6-12) at 
the unknown location of the ionopause, v?hich has to be determined as 
part of the solution. 

This free -boundary problem, like that for the determination 
of the magnetopause, is too difficult to solve as stated. In order to 
proceed, we introduce, as in the corresponding application to the 
Earth's magnetosphere, too important simplifications in these rela- 
tions. In the present application, p tends to be much larger than 
B/s^ on the exterior side of the ionopause. Therefore the pressure- 
balance relation of equations (6-12) may be reduced to 


Pint Pext ( 6 - 13 ) 


a simple balance between the ionosphere pressure, which is given tJith 
sufficient accuracy for the present purposes by equation (6-11) , and 
the exterior pressure of the flox^ing solar plasma adjacent to the 
ionopause. However, p , still remains unknown until the exterior- 
flow solution is determined, but the problem can be reduced to tract- 
able form, as was done for the corresponding determination of the 
Earth's magnetopause, by 'the introduction of the Newtonian approxima- 
tion for the exterior pressure on the ionopause (see equation [6-4]). 
The Newtonian approximation p^xt , however, can be 

improved because it underestimates the pressure ori the flanks of the 
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ionopause ('Spreiter et , 1966a) . Behind the planet where the 
boundary surface becomes parallel to the flow, reaches 90 degrees, 
and the Newtonian expression indicates that the plasma pressure ap- 
proaches zero vjhich is physically implausible. A more reasonable 
limit for the exterior pressure on the flanks of the ionopause is the 
sum of the static magnetic and gas pressure of the solar plasma 
p^ + Bra/SIT because very far behind the planet the pressure must 
return to free-stream conditions. Therefore, a better approximation 
for the exterior pressure on the ionopause is 

Pg^t= {g^-p^-B^8T)C0S^)^ -f Pco~!-^y^^ (6-14) 


where p ^ /^co Vco is the stagnation pressure of the solar wind 

S 1 

exerted on the nose of the ionopause (Spreiter ^ j 1966a). This 

expression together with equation (6-11) and the relation p = 

t \v 

p^+ Baj/S'Tt” from above reduces equation (6-13) for the pressure 
balance to 

= (K^»V|-1^)cOs2Y 


p - p lexp 


r-r 


R 


H 


(6-15) 


So in exactly the same way as for flow past the Earth, use of the 
foregoing approximations greatly simplifies the problem because they 
decouple .the calculation of the shape of the ionopause from the de- 
tailed solution of the surrounding flow and permit the coordinates of 
the ionopause to be computed without determining any further properties 
of the’ exterior flow field; 

It is convenient in the calculation of the shape of the iono- 
pause to let > the distance from the center of the planet to 

the nose of the ionopause. Since cos^"^ = 1 at this point, equation 
(6-15) simplifies .to 
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Pr = Pq = (6-16) 

at the ionopause nose, and to 

COS^'i^ — e ^o)/h (6-17) 


elsewhere along the ionosphere boundary. If, for example, the solar 
wind is considered to approach Venus with a number density of 3 pro- 
tons /cm and a bulk velocity of 5 x 10 cm/sec, as the data reported 

by Bridge ^ (1967) show occurred at the time of the Mariner ~5 

"8 

measurements, equation (6-16) indicates that = 1.1 x 10 
dynes/cm^. Equation (6-6) shows that such a value corresponds, as- 
suming T = 700 ®K, to an ion or electron number density Hj =• D= = 
n/2 at r = Pq of about 5.7 x 10 particles/cm . This value is 
somewhat greater than that indicated by the data of figure 6-4 for 
r = 6500 km, but is not unreasonable if the averaged nature of the 
Mariner -5 data involved in the estimate and the uncertainties attend- 
ant with all plasma measurements in space are considered. In any 
case, it may be seen from figure 6-4 that the essential condition that 
the ion-electron number density required to stop the solar wind be 

exceeded at some altitude is satisfied by a considerable margin, since 

5 3 

the peak daytime value of 6 x 10 electrons/cm is larger than the 
required value by a factor of about 10.’ Similarly for Mars, the values 
Vco = 3.3' X 10 cm/sec, Hco =f 0.8 protons/cm , and T = 200'^K 
inferred from the measurements of Mariner 4 (Lazarus £t , 1967 ; 
Kliore et , 1965) indicate that an electron number density of 
2.6 X 10 electrons/cm is required at the subsolax point for the 
ionosphere to stop the solar wind. This value is exceeded by a fac- 
tor of about 6 by the peak electron density of 1.5 x 10 electrons/ 

3 

cm deduced by Kliore £t aT. (1965) from Mariner 4, and subsequently, 
confirmed by similar experiments with Mariners 6 and 7 (Fjeldbo ^ aJ., , 
1970) . On the basis of the further observations that the altitude of 
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the peak electron density is between 120 and 130 km above the Martian 
surface and that the electron scale height there is between 20 and 
25 Ion, we may estimate from the above relations that the subsolar 
point of the ionopause is at an altitude between 155 and 175 km. for 
these reasons, we feel that the applicability of the theory to both 
Venus and Mars is supported by existing data, even though those for 
Mars are restricted in quantity and inferential in nature at the 
present time. 

To proceed, we must express COS^'f' in terms of r and 
0 , where r(0) represents the, coordinates of the boundary and 

© is the angle measured at the center of the planet with respect 
to a line that extends directly upstream. Upon carrying out this step 
in the way made clear by the illustration in figure 6-5, and substi- 
tuting the result into equation (6-17), we find 

(rd0cosQ + dr Sin Q)^ 

Ids) 'dr^ -i- (r d0)“ 



E 


(6-18) 


The numerical solution of this differential equation is facilitated 
by solving for dp/pd© to obtain 


dr _ sm26± E ^ 

rd0 2(E - sin^ 0) 


(6-19) 


At the ionopause nose 0 = -0, P=Po , and dp/pdO = 0 with either 
choice of sign. The proper choice of sign is dictated by the following 
considerations based on the assumption that P increases monotonically 
from pQ to oo , and hence that E diminishes monotonically from 
1 to 0, as 0 increases from 0 to 180°. Since SI 11^0 = 0 at ©■ = 0, 
increases to unity at 0 '= 90° , and then returns to 0 at 180°, the 



Figure 6-5. View of element of ionopause and 

coordinates used in equation (6-18) 
(Spreiter, Summers, and Rizzi, 1970b), 


denominator of equation (.6-19) must vanish at one or more values for ' 
0 , and dr/rd 0 would be infinite unless the numerator vanishes 

simultaneously. If the critical values for 6 and E are designated 


by the subscript CP , we have 
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and also that 


S! n. 20cr - 2 


2 

E^cr^^crj -0 


(6-21) 


Substitution, of the former into the latter to obtain 


s! n 2 ©cp — 2 


SI n^Ocp - SI n"^0cp| 


i/2 


( 6 - 22 ) 


= 2 sin ©cpcos 0CP •— 2 sinGcpCOsOcp 


0 


shows that the minus sign must be used in equation (6-19). That the 
resulting indeterminate form actually leads to a finite value for 
dr/d0 at © ©CP may be confirmed by application of L.'Hospital’ s 
rule. The shape of the ionopause may therefore be determined directly 
by integration starting from the boundary condition that r Qj 
at 0 = 0 . 

Results obtained by numerical integration are presented in 

terms of cylindrical coordinates X — P COS^O and P =P ©IH© for 

several values for H/r^ from 0.01 to 1 in figure '6-6. As noted 

previously, this range is ample to include all likely possibilities 

1 

for both Venus and Mars. Also included on figure 6-6 is a dotted 


Although the displayed results are independent of the alti- 
tude of the ionopause nose, a planet silhouette has been added to the 
plot to indicate the radius f'v of Venus for the ratio’ 'V/'i “ 

6050/6500 = 0.93 . The corresponding ratio Pm/H> for Mars is not 
knoTO so definitely, but the estimate given above indicates that it 
is about 0,95. 



0 


- 1.0 

k/x^o 


- 2.0 


Figure 6-6. Calculated location of ionopause for various 
H/tQ . The corresponding results for the 
shape of the Earth’s magnetopause nondimen- 
sionalizcd so that the magnetosphere nose is 
at x/ro = 1 , is included for purposes of 
comparison. 
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line indicating the coordinates of the magnetic equatorial trace of 
the magnetosphere boundary, as determined by Beard (1960) and Spreiter 
and Briggs (1961, 1962), that has been used extensively in the calcu- 
lation of solar -wind' flow past the Earth. These coordinates have been 
nondimensionalized by dividing by the distance from the center of the 
Earth to the magnetosphere nose. We shall continue to call this dis- 
tance pQ , even though it is much larger than for Venus or Mars, 
and evaluated differently using the expression P^ — 

^(BsqjSirK in which = 6.37 x 10^ cm is the radius of 
the Earth, and = 0.312 gauss is the average intensity of the 

geomagnetic field at the geomagnetic equator . It may be seen that 
this curve is very similar to that for the ionosphere boundary- for 
H/Qj = 0.2 . The size, with respect to the planet, of the cavity 
carved in the solar v?ind is very different in the two applications, 
however, since P^ is only a few percent greater than the planetary 
radius for Venus or Mars, whereas it is usually of the order of 10 
Earth radii for the magnetosphere (see figure 6-3). 
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CHAPTER VII 


NUMERICAL SOLUTION OF THE APPROXIMATE EXTERIOR FLOW 


A . Simplification, to Gasdynamles 

We have just shown in the preceding chapter that the exterior 
and interior regions of the complex free-boundary problem represent- 
ing the solar -x^ind interaction xjith planets can be disjoined and solved 
separately and that the interior boundary shape can be calculated inde- 
pendently of the exterior-flow solution. We have determined the shape 
and location of both the magnetopause and ionopause, and xirhat remains 
is to solve the exterior -flox>7 problem represented by equations (5-1) 
satisfying the remaining tangential discontinuity relations Vn = Bh'= 

0 of equations (6-12) at the location of the previously determined 
boundary . These reduced tangential discontinuity conditions are com- 
pletely equivalent to the boundary conditions that must be satisfied 
on a solid body, and the problem may be thought of as a hypersonic, 
but x^eakly magnetized plasma streaming past a solid body having the 
shape of either the magnetopause or ionopause. In many x<rays it is 
very similar to the aerodynamic problem of a body in a supersonic 
stream. 

Hox^ever, this magnetohydrodynamic-f lox>7 problem is far too com- 
plicated to be solved exactly, and resort must be had to approximations. 
For Earth, Mars, and Venus, an important simplification of the magneto- 
hydrodynamic equations (5-1) can be made on the basis of the combined 
effects of the large magnitude of the Alfven Mach number and the 
strong interaction nature of the flo\^ (Spreiter ^ , 1966a). This 

is so because the order of magnitude of the inertia term ^{V • V) V 
in the differential equation for the momentum is related to the mag- 
netic term (l/4ii')B XCUpI B by the square of .the Alfven Mach number. 
When the latter is large, therefore, the magnetic terra can usually be 
safely dropped from the equation xxfith little loss of accuracy, just as 
discussed in section B of chapter VI. In order for this approximation 
to be valid, hoxvever , the square of the Alfven Mach number must remain 
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large at all points in the exterior flow. For normal conditions of 
the solar v?ind upstream of the how vrave, this term is about 100, but 
observations made in space indicate that it decreases somevjhat in the 
disturbed region of the flow downstream of the bow wave. Lees (1964) 
has suggested that a small region may exist near the stagnation point 
in which the magnetic field may become sufficiently enhanced under 
certain circumstances to produce a substantial reduction in the gas 
density. The region affected is very small for weak magnetic fields, 
however, and is anticipated to have little effect on the gross fea- 
tures of the flox^ so that the simplification is justifiable throughout 
the exterior -floxj region. Although factors of 2 appear v?hen similar 
comparisons are made of the magnetic and inertia terms in the conser- 
vation equations for the momentum and energy, the general conclusion 
that the magnetic terms may be dropped from these equations when the 
Alfven Mach number is large still holds. An immediate consequence is 
that both the differential equations and the conservation equations 
defining the floxv are decoupled from those involving the magnetic 
field. The differential equations (5-1) for the fluid motion thereby 
reduces to 


\7 . pV - 0 

r~^ I rJ 

/?(v ♦ V) V Vp = 0 


(7-1) 


(V*V)S = 0 


S 3 q 


P/Pq 

\PlPof 


and the conservation equations (5-5) to 


[pVn]=0 
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+ 


pn 


= 0 


/flVn(h-4v2) 


0 


(7-2) 


Thus, all the properties of the flow except the magnetic field B 
can be determined by solving the equations of gasdynamics (7-1, 7-2). 
The properties of the magnetic field can then be determined subse- 
quently by solving the remaining equations 


curl (V X B) r: 0 


dlV B rr 0 



n 



Bt V 


n 





(7-3) 


using the values for V already computed from the preceding step. 
These relations, or equivalently, the following equations derived 
from them 



Drl/^ I 



(7-4) 


in which is an arbitrary surface moving with the fluid and 
D/or= |r+(y0y =|y-v) V is the substantial derivative for 
steady flovj, indicate (L'audau and Lifshitz, 1960) that the magnetic 
flux passing through 51 is conserved, and moves with the fluid. In 
more picturesque language, the magnetic field is frozen in the fluid, 
the .general effect being essentially similar to that of a line of 
smoke or dye released at some instant into a medium flowing at speeds 
high relative to the diffusion rate of the contaminant. 
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B . Description of Numerical Method 

The preceding discussion has shovjn how the nonmagnetic proper- 
ties of the exterior flow described by equations (5-1) can be repre- 
sented without undue loss of accuracy by purely gasdynamic equations. 

In this manner, the flow of the solar wind past the Earth, Mars, and 
Venus is exactly equivalent to the familiar, but complex, problem of 
the external aerodynamics of round-nosed bodies in a supersonic stream. 
Although this represents an enormous simplification, the mathematical 
complexity associated with the non-linear and mixed elliptic -hyperbolic 
character of the governing partial differential equations requires that 
the solution be sought by numerical methods. A variety of numerical 
procedures can be used to solve the equations of gasdynamics for 
flow about a given body in a supersonic stream, but the difficulti.es 
are sufficiently great that at the present time solutions can only 
be computed for round-nosed, axis 3 nranetric bodies. This restriction 
poses no difficulty for the calculation of flow past the ionopause 
which is a surface of revolution, but it remains necessary to approx- 
imate the magnetopause with an axisymmetric shape. In nearly all 
cases for which calculations have been carried out, the shape selected 
is that obtained by rotating the magnetic equatorial trace of the 
boundary illustrated in figure 6-3 about its center line. Inspection 
of the magnetopause coordinates shows that this approximation is 
reasonable, except possibly where the boundary is dented inward near 
the neutral points. Even this exception may be of little importance, 
however, because these dents result from the use of the Nex^tonian 
pressure approximation and XTOuld disappear in a more accurate analysis 
based on fluid concepts (Spreiter and Summers, 1967). In our calcula- 
tions x7e use the shape obtained by rotating the magnetic equatorial 
trace of the boundary given by Spreiter and Briggs (1961, 1962) . 

Several methods are available for the numerical solution of 
the equations of gasdynamics for flox7 about a given round -nosed body 
of revolution in a supersonic stream. The method employed in the 
calculations for our x^ork here is that developed by Inouye and Lomax 
(1962) and Lomax and Inouye (1964) x^hich is founded on the basic 
approach of Van Dyke (1958) and Van Dyke and Gordon (1959) as modified 
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by Fuller (1961). As a representative example, figure 7-1 illus- 
trates the gasdynamic equations (7-1) and (7-2) used to calculate 
the flow around the magnetosphere and the cylindrical coordinate 
system originally presented in figure 6-6, As indicated, the method 


DIFFERENTIAL EQUATIONS 
V*pV=0 • 
p(V*V}V+Vp = 0 


(V-V)S = 0 
S“ So =Cy %n 


P/Pq 



CONSERVATION EQUATIONS 

[pVr,V + pn]=o 
[pVn(h^+v^/2)]=0 



WHERE [q] = Q,-Qq 


Figure 7-1. Gasdynamic equations for steady flow of 
a dissipationless perfect gas (Spreiter 
et , 1968) . 


used for the subsonic region. near the magnetosphere nose and the 
immediately adjoining portion of the supersonic region is an indirect 
one in which the location of the bow-shock wave and the free -streamline 
conditions are assumed knoxm and the associated flow field and body 
shape are found as part of the solution. The desired solution for 
the specified body, either the magnetopause or the ionopause, is then 
found by iteration following judicious selection of the initial 
trial shape for the bow wave based on experience with a vast number 
of cases of aerodynamic .interest , The solution for the remainder of 
the supersonic region is accomplished directly by using as indicated 
the method of characteristics in the manner described explicitly for ■ 
this application by Inouye et al. (1965). As a background to our 
later modification of these schemes and discussion of the computed 
results, it is appropriate to present at this time a detailed review . 
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of the mathematical and numerical aspects, such as convergence, 
stability, and accuracy for both of these methods. The following 
discourse is essentially repeated from Lomax and Inouye (1964) and 
Inouye ^ (1965). 


r. Method for subsonic-transonic region 

In forming difference equations from equations (7-1), a more 
useful form of the energy equation can be developed from the fact 
that entropy is constant along streamlines . For a gas in equilibrium, 
the pressure is determined from the equation of state using any two 
thermodynamic variables, say S) Then using the chain 

rule of differential calculus we find 


but 


Dp _ 

or 


'^1 ^ 

, a Pis or 


4- 


ps 

las)/?0T 


since ps =0 
or 


we finally obtain 


DP _ ^2 D_^ _ 
Dr DT 


(7-5) 


where c3 = is the local speed of sound. In terms of the 

cylindrical coordinate system {XiP) described above, the components 
of the first two vector equations of (7-1) and scalar equation (7-5) 
are 


ar 




q. 


a X 


a X 


+ 42=0 

r 


■ 


a X 


0 


Although we have already specified the equation of state to 
be that of a perfect gas, the explicit form of that equation is inde- 
pendent of both the form of the difference schemes and the numerical 
procedures. We therefore use the more general expression for the 
thermodynamic equation of state p = P(^^S) 
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(7-6) 
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where U and "tr are respectively the velocity components in the 
X and r directions. These four partial differential equations 
must be solved simultaneously for -the four dependent variables p 
p j U 5 and U . The thermodynamic relationship 


a .= d{|0,p) 


(7-7) 


which is derived from the definition of a by using the equation of 
state provides -the fifth equation connecting the five unknowns p 
p , U j If , and a . 

In the subsolar region of the planetary boundary, equations 
(7-6) exhibit different character; namely, the equations are elliptic 
in the subsonic region, parabolic on the sonic line, and hyperbolic in 
the supersonic region. Despite these complications, Lomax and Inouye ■ 
(1964) have perfected an inverse method in which a shock shape is as- 
sumed and the equations are integrated numerically by a finite-dif fer- . 
ence method to determine the corresponding body shape . Their particular 
scheme has proven to be accurate and efficient in solving such flow 
fields, and we have used the computer program which is comprised of a 
main program and 26 subroutines written by them at Ames Research Center. 

The complexity in the details of programming difference equa- 
tions vJith more than one independent variable depends critically both 
on the nature as well as the manner of treating the boundary conditions. 
Tor this reason the governing equations (7-6) are transformed to a co- 
ordinate system that provides the simplest application of the boundary 
values. Because the initial-boundary conditions are specified along 
the shock, the axis of one set of coordinates coincides as shoiiTn 
in figure 7-2 with the shock, and the coordinates are those lines 
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placement of the shock in 
the free- stream direction. 
The other set of coor- 
dinates is parallel to the 
free-stream direction. 

This coordinate system is 
obviously convenient for 
starting the inverse ,problem 
although not optimum for" 
studying body shapes . The 
coordinates are equally 
spaced in both directions , 
and their intersections 
define the points for the 
difference mesh. If the 
analytic expression for the 
shock shape is X =X(r) , 

then these new coordinates 
S and + are defined by 

s = x-X(f) 
t - f 



Figure 7-2. Sketch of flow field 
around blunt-nosed 
body and coordinate 
system' (Lomax and 
Inouye, 1964). 


(7-8) 


and by the usual application of calculus, the X and p partial 
differential operators x^ith respect to X and P become 

S ^ ^ — xr*'' ^ 

~ as a r as at' 

Equation (7-6) can then be transformed and expressed in matrix form 


[a] • [B] = [-C] 


(7-9) 
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where 



0 u -z V ' p pi! 

1 0 /O(u-Zv) 0 
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Equatioa (7-9) can be inverted, and the partial derivatives with 
respect to S of the four dependent variables can be expressed as 


€ (f P P U IT ^ ^ ^ 

t, r ? ^ st ^ ^ at 


3y 
5 3t 


^ - 
ss 


f 2 [t, P,/°, 


3 /O 

’ at 


'bv \ 

at •) at j 


s 


f3(t,P,/’, u, V, 
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at 
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, at J 


(7-11) 


a V 
as 


= f4(t P> 
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5U 
, at 



where > fa , f 3 , and represent complicated expressions 

involving the indicated variables and are the four elements of the 
matrix derived by performing the operation 

For a given shock shape .and free-stream conditions, the values- 
of P , j U , and V just behind the shock wave are calculated 

from the discontinuity relations (7-2). Lomax and Inouye (1964) then 
divide the numerical calculations into two parts. First, the deriva- 
tives with respect to t of P , P , U , and Lr are calculated 
numerically from their known values along a t coordinate. They are 
determined by a standard five-point central difference method except 
at the upper end of a t coordinate where a skewed five-point differ- 
ence scheme is used. This information together with equation (7-11) is 
then used to advance the solution in the 5 direction by a predictor— 
corrector procedure. A second— order method is used and is illustrated 
for a typical flow variable P as follows: 

(1) Numerically differentiate the flow properties 
for the i th step to obtain 

W+ ii 

(2) Calculate from equations (7-11) 



(3) Predict new value f^ 4 .|~ P]i v-j 2 A S ^ 




5t 


(5) Calculate 


f i + i 
5P ' 


ss 


il-H 


(6) Correct new value Pl^v-l “ pj^ 4- 


■AS 


) obtain 

( 

(7-11) 
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Usli ^ 

Us 


(7-12) 


In this manner the kno^m values of the flow properties are advanced 

2 

by marching step-by-step towards the body. The stream function 

is calculated for each point, and a body is determined as the 
locus of points where the stream function vanishes. This body shape 
is then compared with the desired shape, and if the fit is not satis- 
factory, then the shock shape X(T) is varied slightly, and the 
whole procedure is repeated. In practice several such iterations 
yield the desired body shape to a very accurate degree. 

a. Convergence and stability . The question of conver- 
gence is considered by Lomax and Inouye (1964) in two parts. First, 
as the mesh size is reduced, do the difference equations converge to 
the differential equations, and second, if they do, does the calcula- 
tion procedure itself converge throughout the region of application? 
The complexity of the governing equations makes the answer to these 
questions impossible from a purely mathematical approach. From a 
physical viewpoint, however, some insight to them can be achieved. 

With regard to the first question, if the difference equa- 
tions do not reduce to the specified differential equations , then they 


2 

The solution is started by a first-order predictor followed 
by two second-order correctors. 
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must reduce to some set of differential equations which violate one 
or more of the conservation lav?s. To check whether or not these 
lax\TS have been violated, however, can be done independently of a 
mathematical analysis of the difference equations. For example, 
the body surface is located by the condition of mass conservation and 
if either momentum or energy is not conserved, the entropy and total 
enthalpy calculated at the body location X'yould not be constant. 

These conditions x^ere checked, and for the solutions presented here, 
the entropy along streamlines xjas constant to an accuracy consistent 
V7ith the calculations, and the total enthalpy v?as constant throughout 
the floxv7. These are independent checks because entropy and enthalpy 
are not used in calculating the flov? field, and they insure that the 
results of the finite difference methods employed do correctly repre- 
sent solutions to equations (7-6) and (7-7). The second question of 
whether or not the calculation procedure itself converges is treated 
by the follox\?ing argument. If the calculations give ansx^ers that are 
physically consistent i.e,, satisfy the independent physical checks 
and are not significantly affected by changes in either the mesh size 
or mesh ratio, then we assume they have converged and represent the 
correct ansxjer . 

Hoxsrever, considerable difficulty is encountered in the study 
of the flox7 properties behind a prescribed shock with the flox7 proper- 
ties ahead of it knoxra. The difficulty arises because analytically 
our approach is only well-posed for an- initial -value or Cauchy-type 
problem, and x-je are attempting, in mathematical terms, to solve a 
boundary-value problem by an approach suitable for an initial-value 
problem. The difficulties arising from such an attempt lead to some 
of the most fundamental problems in the numerical treatment of both 
partial and total differential equations. 

In the solution of equations (7-6) and (7-7) by means of 
equations (7-12), numerical instabilities are easy to detect once 

i 

they have started. A convenient demarcation of their onset is the 
locus that separates regions x^here the variables axe at least physi-' 
cally possible from regions xvhere they fluctuate beyond all reason 
predicting negative pressures, densities, etc. Unfortunately, the 
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complexity of the problem makes complete mathematical rigor in these 
studies practically impossible, and one is forced to rely on experience 
with linearized equations and familiarity with the physical problem 
for help in making the arguments plausible. Nevertheless, Lomax and 
Inouye (1964) have been able to classify all the difficulties encoun- 
tered in this study with one of the following categories (see Fox 
[1962] 'for terminology and further background): 
i. Inherent instabilities 

(a) Due to ill conditioning (nonessential) 

(b) Due to singularities (essential) 
ii. Induced instabilities 

Inherent instabilities are brought about because (a) the differential 
equations themselves contain an unstable (exponentially growing) 
solution in the direction of one independent variable while remaining 
bounded in the other, or (b) because the differential equations con- 
tain a solution that is singular. Such instabilities are not caused 
by the finite differencing. Induced instabilities are brought about 
by the particular numerical techniques employed, mesh size, degree of 
truncation, implicit or explicit methods, etc., and can lead to im- 
plications quite spurious with regard to true solutions to the partial 
differential equations. ^ 

Following Lomax and Inouye (1964) we will refer to the two 
different kinds of inherent instabilities mentioned above as essen- 
tial and nonessential. The terminology is, perhaps, not apt because 
what is referred to as an essential instability is actually caused by 
the appearance of a singularity or group of singularities in the flow 
region between the shock and body. The singularities are Invariant 
to coordinate transformations and represent the locations of sources 
or sinks that would appear in an exact analytic solution behind an 
analytic shock. Nevertheless, they refer to the numerical behavior 
caused by them as an instability because of its similarity in ap- 
pearance to other numerical phenomena which are identified by that 
terra. The important distinction between an essential and a nonessen- 
tial instability is that the latter would not occur in an exact analytic 
solution. Nonessential instabilities are started by round-off, 
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truncation, or end -of -array inaccuracies due entirely to the fact 
'that numerical methods are employed. In most of the literature per- 
taining to the stability of partial differential equations, the state- 
ment is made or implied that initial-value data are always unsuitable 
for elliptic equations (see Hadamard [192^). In the present termin- 
ology it is the existence of nonessential instabilities that leads 
to this conclusion, a conclusion that is valid if such instabilities 
cannot be controlled. 

A great deal has been x-Tritten on the effect and control of 
what we call nonessential instabilities. So far as xre know, most of 
the reported discussion attempting to identify’ these instabilities 
with the local form of the governing equations (elliptic, parabolic, 
or hyperbolic) have been limited to the elliptic or subsonic region. 
Although nonessential instabilities could be detected in the subsonic 
portion of the flow field if the numerical differentiation in the 
"i' direction was of very low order and no smoothing was used, Lomax 

and Inouye (1964) encountered 
no difficulties in this area. 
According to them the really 
critical area for stability 
was located in the supersonic 
portion of the flow, roughly 
in the shaded region of fig- 
ure 7-3. In fact, in this 
region they were able to 
isolate all three of the 
major types of instabilities 
listed above . And although 
the difficulties inherent in 
the asymmetric differentia- 
tion formulas used at the 
upper end of the arrays often 
accentuate' the breakdown, 
they are not the cause of 
these instabilities. 



Figure 7-3. Sketch of flow field 

showing critical region 
of instability problems 
(Lomax and Inouye, 

1964) . 
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Consider the flow behind the shock in -the area where it is 
supersonic. As we inarch inward from one t coordinate to the next, 
xo'e are perpetrating an initial-value problem. We use the term 
"properly set" in the usual way; that is, Cauchy data (function and 
derivative) are properly set for the wave equation along the + = 0 

axis, and Neumann or Dirichlet data (function or derivative) are 

% 

properly set for Laplace's equation. Usually, stability proofs for 
hyperbolic equations stem from the assumption that initial-value or 
Cauchy data are properly set. Hadamard (1923) pointed out, and it is 
well known in studies of supersonic wing theory, that according to 
hoxv the line carrying the initial data cro . >'.s each characteristic, 
Cauchy data may or may not be properly set for a hyperbolic equation. 

By Hadamard' s terminology, a data line is "duly inclined" if Cauchy 
data are properly set, and "nonduly inclined" if they are not. One 
can easily show that a shock is a nonduly inclined surface in the 
region inhere the flow is supersonic behind it (see point A, figure 
7-3), and Cauchy data are, therefore, nox^here properly set along it. 

In fact, the t coordinate continues to be nonduly inclined until 
the upgoing characteristics lie to its right as is shoxvn at point B 
in figure 7-3. 

The principal point made by Lomax and Inouye (1964) in this 
regard is that the nonlinear equations governing the flox-j field may 
have inherent instabilities in the supersonic as x^ell as in the sub- 
sonic regions. In fact much of the numerical discussion, as x^ell as 
most of the criticism, of the inverse method concerns the existence 
and manner of treatment of the nonessential instabilities. The defense 
usually invoked by those who use the inverse method is, simply, that 
the gas layer between the shock and the body is thin enough that the 
errors caused by numerical calculations cannot grox^ sufficiently 
large to invalidate the first few significant digits in the results. 

We adopt the hypothesis of Lomax and Inouye (1964), that this argument 
is sufficient to validate those cases for x-ihich the solution passes 
the consistency checks on total enthalpy and entropy as discussed 
above . 

Absolute reliance upon this hypothesis is not alxiiays 
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(a) Pressure distribution. 



(b) Fourier coef fieients. 


Figure 7-4. 


Harmonic analysis of 
pressure distribution 
just behind bow shock 
wave (Lomax and 
Inouye , 1964) . 


satisfactory in application, 
and is certainly not satis- 
factory from a theoretical 
point of view. Basically, 
the validity of such a - • 

hypothesis is coupled with 
the word size (number of 
significant digits carried 
in each arithmetic operation 
and stored in memory) avail- 
able in the computing machine . 
It is not always wise to 
demand larger computing capa- 
city merely to push ahead a 
few more steps before ex- 
ponentially growing instabili- 
ties started by numerical trun 
cations swamp the first few 
significant digits in the cal- 
culations. A more sensible 
approach is to face the 
problem with analysis and 
attempt to suppress non- 


essential instabilities by appropriate numerical methods. One such 
method that has been used by several authors (Fuller, 1961; Lomax and 
Inouye, 1964) is to "smooth" or filter the data along each + coor- 
dinate as the computations proceed. Because of the important role 
this method plays in controlling instabilities, we shall present some 
of the discussion of Lomax and Inouye (1964). 

A representative variation of pressure along a 'b coordinate 
• and behind a curved chock is sho\m in figure 7 -4 (a). If we expand the 
pressure distribution in a Fourier series betxreen the maximum absolute 
values of i" used in the calculations, we obtain 
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Eei = f hcosMl 
p <i_i n f . 

max n=o 


-■^max< + <-Ux . 


The magnitudes of the coefficients are shown in figure 7"4(b). We 
see that the first few terms dominate the expansion, and higher-order' 
terms, corresponding to higher frequencies, are negligible. Next we 
assume that a Fourier expansion for any of the dependent variables 
over the same interval along any ' 1* coordinate between the shock and 
the body would exhibit the same general behavior, that is, could be 
expressed with acceptable error by the first few- terms in a Fourier 
expansion. 

That assumption certainly warrants some discussion. However, 
to attack it simply on the basis that it arbitrarily prohibits high- 
frequency terms in the true solution is not justifiable, since any 
method using finite-difference techniques is subject to such a 
criticism per se . In fact, as pointed out by Lomax and Inouye (1964) 
it is a fundamental theorem in communication theory (Shannon and Weaver, 
1949) that the highest resolvable frequency in any finite trigonometric 
series to be represented by a discrete number of points is related to 
the number of points themselves. Consequently, they reasoned that if 
more points are used than are required to resolve the highest fre- 
quency, the remaining points carry redundant information. The real 
assumption, then, is not that high-frequency terms are arbitrarily ex- 
cluded, which is the case for numerical calculations in any event, 
but rather that enough mesh points are taken to make data' contained 
in them largely redundant in the sense just mentioned. Under these 
circumstances high-frequency terms, appearing from more or less ran- 
dom errors brought about by numerical truncation in calculations made 
at discrete points, fall above the frequency range required to ex- 
press the true solution and can be excluded if the proper numerical 
filter is applied. The following discussion briefly describes such a 
filter used by Lomax and Inouye (1964) . 

Consider an even function F(x) that can be represented by 
a Fourier cosine series in the interval O^X < L as 
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F(x) = ^ 




nmrx 

[T 


(7-14) 


and define the operator fj by the equation 


fj(F) = F(x+Jdo) + F(x-Jdo) 


(7-15) 


where is the spacing of the points in the X direction. With 

W I defined as arbitrary weighting factors, the sequence of opera- 
tions ^ applied to equation (7-14) results in the expression 
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where 
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If we define a polynomial by 
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(7-17) 


then equation (7-16) becomes 
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And 


n repetitions of the operation yield the expres 


Sion 


F (x) 


CO 


m=o 
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L 


(7-18) 


Thus, Lomax and Inouye (1964) concluded that if the function F(x) 
is replaced by certain weights of its average at equally spaced inter-' 
vals to the right and left of X (i.e., applying the operator 

"fj ) , the coefficients in its Fourier expansion are multi- 
plied by which is independent of X and depends only on the ar- 
bitrarily assigned weights Wj and the factor . Furthermore, 

if the operator is applied H times, each coefficient is multiplied 
by the H th power of P.<s . • ‘ 


If the terms 


can be chosen such that 


J , Wj , and 
is near unity for values of m less than some arbitrary fre- 
quency j and close to zero for higher values of fTl greater than 
some Mg and falls in the range 0 < Rp < I for intermediate m 
where < tn < , then will be suppressing the higher-fre- 

quency terms and be acting like a filter. We can therefore easily 
see that the corresponding weighted-average operation will leave the 
low harmonics for which 0 ^ hi ^ unchanged and x<?ill destroy all 
the higher harmonics for which in>M 2 . Harmonics for which Mi< 
iTt <C Mg \rLll be distorted according to the nature of the poly- 
nomial Pp> . To illustrate that P^ actually does possess such 
filtering properties, consider equation (7-17). An examination of 
this equation reveals it to be a truncated Fourier cosine series 

with summation index J . If we choose Wo = 1/2 and Wj = 

a 5m Jit /2, ^ , 

tor J “ 1, 2, J... , then rp is simply a finite-term 

Fourier expansion in variable III do approximating a positive square 

wave which clearly has the desired characteristics of a filter. To 

improve the approximation and make the filter sharper, either more 

terms of the series can be taken, or, as in our case, higher powers ’ 

of P^c corresponding to repeated averaging can be used. 

In our, application , do is unkno^m. We simply assume that 
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an L exists for which our numerical results can he represented with 
acceptable error by harmonics below M| and for which our numerical 
procedures will introduce errors consisting of harmonics above 
(but below those very high frequencies for V7hich fT) > )• By 

numerical experimentation with the spacing do and the weights Wj , 
Lomax and Inouye (1964) have verified that this is so. The actual 
polynomial which they found appropriate and which we have used in our 
calculations is 
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(7-19) 


Each application of Rp" amounts to moving the central point in a 
group of five to the curve obtained by a fitting a least-squares 
quadratic to them. 

The use of numerical filters of this type can be contrasted 
to the use of higher-order numerical differentiating processes. Con- 
sider, as sho™ in figure 7 -5 (a), a set of data all zero except for 
one "bad” point. These data were numerically differentiated by Lomax 
and Inouye (1964) using central-difference schemes and the results 
harmonically analyzed. Their results given in figure 7-5 (b) show that 
with higher -order difference methods the amplitudes of higher- frequency 
terms are increased and the over-all maximum amplitude itself is in- 
creased. On the other hand, a five-point difference scheme followed 
by repeated filtering, of the type just described, has the opposite 
effect. This does not mean, of course, that relative to low-order 

schemes, high-order differentiating processes are inferior. Their 

* / 

value comes in reducing the initial error. But once this error has 
been committed, the higher-order methods are of no value in control- 
ling its effect on the stability. 

To demonstrate the growth of an inherent instability and the 
effect of numerical filtering, we present some results from the study 
by Lomax and Inouye (1964). Four curves are sho'sm in figure 7-6, one 
corresponding to each of the four possible combinations of tv?o spacings 
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(a) Set of data with one "bad" point. 
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(b) Harmonic analysis of data after 
differentiation and filtering. 


Figure 7-5. Effect of differ- 
entiation and filter- 
ing on data (Lomax and 
Inouye , 1964) . 
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Figure 7 -6 . Effect of interval 
size on determinant 
[a] in equation (7-10) 
(Lomax and Inouyg, 
1964) . 


for AS and At . The quan- 
tity plotted is the magnitude 
of the determinant of [A] from 
equation (7-10), the value of 
which is, of course, critical 
in evaluating the S deriva- 
tives of the dependent vari- 
ables. (Actually, any one of 
the dependent variables could 
be used to display the same 
result.) Plotted along the 
abscissa are the points located 
along a 1" coordinate in the 
supersonic region above the 
body and doT-Tn stream from the 
shock. Notice that the in- 
stabilities have no correla- 
tion with the ratio of At/AS 
(all other conditions were 
held constant) , .but rather 
with the absolute value of At . 
For the larger At no fluctua- 
tions at all are observed, 
whereas, for the smaller At , 
unstable oscillations begin at 
a t of around 0.6 , and the 
distance between successive 
maxima and minima is about five 
times the At interval size. 
Lomax and Inouye (1964) ex- 
plained this behavior by es- 
timating the magnitude of the 
polynomial p.p' . With the 
smaller "t spacing correspond- 
ing to a smaller value for 
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do > equation (7-19) indicates that is about 0.99 for one 

application of the filter and about 0,90 for 10 successive applica- 
tions. Hence, the amplitude of a term vjith the frequency correspond- 
ing to the oscillating curves presented in figure 7-6 is only slightly 
damped' by the filter represented by . A similar study for the 

larger t spacing in figure 7-6 shows that the same frequency t-jould 
be heavily damped, “especially for multiple applications of the filter. 

The extent of the curves shorn in figure 7-6 is well within the region 
x<Fhere the h coordinates in the supersonic floxj are nonduly inclined, 
a region defined by the inequality det^Aj!>0 . This example illus- • 
trates a typical case of a nonessential instability and the manner' by 
which it can be controlled . 

In their report Lomax and Inouye (1964) also examine and dis- 
cuss an example -of an essential instability which is caused by a line 
singularity in the flow field. They found that this singular behavior 
corresponding to mass absorption xjas directly caused by the use of an 
elliptic shock -xcfave shape S(r) that exactly elliptic shocks 

in the nose region simply do not occur ahead of blunt bodies in 
source-free flov?. The difficulty is easily overcome, however, simply 
by starting with an analytic shock shape that is not elliptic and 
never becomes tangent to a free -stream Mach line. Actually, the shock 
shape used in their work asymptotically approaches tx-rice the slope 
of a free-stream Mach line, and x^ith its use, essential instabilities 
such as they discuss disappear. 

The question of just hox7 close the shock can come to the free- 
stream Mach line and just hov3 it should behave at infinity 'does not 
have to be settled for most bodies, because xve need only to enter the 
supersonic region far enough to provide information x^hich can be suc- 
cessfully continued by the method of characteristics. Usually this means 
that the body shape need conform to its desired value only up to a 
point x-7here the local Mach number is around 1.05 . Past this point 
the given body shape can then be used explicitly because the method of 
characteristics is a direct method. The exact shape of the body past 
•this point, however, does not affect the solution in the nose region 
because no disturbance in the supersonic region can propagate upstream 
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into the subsonic part. Thus, the solution in the nose region is 
quite independent of the details of the flow further do^mstream in the 
supersonic part. The real importance of this whole concept lies only 
in the fact that analytic shock equations for such bodies with source - 
free flow ahead of the -limiting characteristic are not difficult to 
construct. 

The instabilities discussed so far existed because of the 
nature of the differential equations. They would actually have existed 
in a frictionless flow if the boundary conditions were exactly repro- 
duced. Another kind of instability arises which is caused by the par- 
ticular choice of numerical procedure. This kind of instability is 
referred to as an induced instability, and it can presumably be elimi- 
nated by the proper choice of coordinates, mesh size, and differencing 
technique. Lomax and Inouye (1964) examined this aspect of the problem 
as follows. 

The study of the stability of difference methods, as they 
apply to hyperbolic equations, usually pertains to the continuation of 
Cauchy data given along a duly inclined line. The criterion given for 
the stability of the simplest' explicit difference scheme in application 
to the linear two-dimensional wave equation (Fox, 1962) is illustrated 
in figure 7 -7 . The scheme is 
based on a three-point central- 
difference formula written about 
the point m,n for the first 
derivatives in both directions. 

Given equally spaced data along 
column n , the method is stable 
if, as the calculation proceeds 
from point a to point b in 
figure 7-7(a), the columns are 
spaced by A X j , and unstable if 
spaced by .AX 2 . The condition 
for stability is related, simply, 
to whether the point at the m th 
row and (n+l) th column does or does not lie in the shaded area bounded 



Figure 7-7. Sketch of conditions 
for stability (Lomax 
and Inouye , 1964) . 
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by the characteristics shown’. If this criterion were to apply to our 
case (when the n th column became duly inclined) , point b in figure 
7 -7(b) would have to fall in the shaded region shoxjn. Actually our 
case is much more complicated. A five-point difference scheme is used 
along the columns and, after smoothing, a second-order corrector equa- 
tion is applied. 

Coupled with the ever present nonlinear form of our basic 
equations, these complications make a stability analysis of the situa- 
tion quite involved. However, Lomax and Inouye (1964) have observed 
that whenever the flow crosses into a region where a t coordinate is 
duly inclined, the calculation becomes unstable. When the i" coor- 
dinate changes from a nonduly to a duly inclined line, the value of 
det [a] in equation (7-10) passes through zero. Since det [A] is the 
denominator in the equation for the fon-?ard -marching derivatives, this 
can lead to difficulties in numerical methods even though, of course, 
the numerators also must vanish. As was the case for the end-of-array , 
asymmetric, numerical differentiation, this occurrence may help start 
the growth of instability, but it is not the cause of it. In their 
solutions , this crossing only occurs in the high supersonic region so 
no attempt was made by Lomax and Inouye (1964) to control it. The 
line along which the crossover occurs simply becomes a boundary past 
which calculations were not continued. This way of treating such an 
instability is only an expedient. Blunter bodies and solutions at 
lower Mach numbers demand that these instabilities be controlled. Con- 
sequently, at the present time, this particular computer program is 
incapable of solving such cases, fortunately, however, the bluntness 
of the magnetopause and ionopause shapes and the high free-stream Mach 
number that characterize our applications fall within the range for 
which accurate solutions can be computed by this program so that fur- 
ther development is not required. 

b. Relating the shock and body shapes . Even with all 
the stability and convergence problems controlled where necessary, the 
real success of the inverse method still depends upon the ability to 
find some relationship between families of shock and body equations. 

In other words, it depends upon the ability to formulate a shock 
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equation that can be systematically modified to a point where it yields 
a given body to a prescribed accuracy. 

The fruition of the shock-shape iteration requires an effi- 
cient way of comparing the calculated body shape with the desired shape 
of the magnetosphere or ionosphere boundary and a means of minimizing 
the difference. However, the fact that the magnetopause and ionopause 
shapes were numerically calculated and are not expressed by an analytic 
expression makes a systematic comparison very complicated. To reduce 
these complexities, it was found that the nose of both the magneto- 
pause and the ionopause can be represented extremely accurately by a 
fitted ellipsoid. This simplification made the comparison straight 
forward. The following method was then used to relate the associated 
bow shock wave with this boundary. 

As Lomax and Inouye (1964) pointed out, elliptical or spheri- 
cal shock shapes introduce singularities into the flow field which can 
cause essential instabilities to occur in the low supersonic region near 
the body. In order to avoid these difficulties, they attempted to rep- 
resent the general shock shape X{r) by a simple ratio of polynomials 
in r . Since the body shape and size corresponding to the given shock 
are not known at the start of the computation, they found it convenient 
to reference the X and P coordinates to the radius of curvature 
of the shock wave at p = 0 . 

The general equation for an ellipsoid, which is to be fitted 
to the subsonic region of the given planetary boundary, can be written 
in terms of the present coordinate system as 
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(7-20) 


The sjraibol A is the shock- standoff distance, R^, is the radius of 
curvature of the body at P = 0 , and is the body bluntness para- 

meter obtained ’from the fit to the specific planetary boundary under 
consideration and defined by B^^={b/R)^ where b and R 


are 
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ellipsoidal dimensions illustrated 
'in figure 7-8. 

In relating the body shape to 
the shock shape, two factors must be 
considered. First, the ratio Rb/*^s 
which relates the size of the calculated 
body to the shock shape must be deter- 
mined, and, second, the closeness of 
fit of the calculated body points to 
the desired body shape characterized 
by must be determined. The pro- 

cedure devised by Lomax and Inouye (1964) 
and used in the present application is 
to first assume values for the shock 
shape parameters and then calculate 
the flow field including the body shape. 

The body coordinates for surface Mach numbers between 0.5 and about 
1.05 are then substituted into the left side of equation (7-20) and 
plotted against (x~AyR _5 These points are fitted with a horizontal 



Figure 7-8. 


Sketch of el- 
lipsoid (Lomax 
and Inouye , 
1964) . 


straight line that minimizes the sum 


E, 


of the absolute differences 


between the points and the line. The ordinate of the line is the 
ratio Rb/Rj > ^4 the sum of the absolute differences is a measure 
of the closeness of fit to the desired planetary boundary. If neces- 
sary, any given parameter in the shock equation can be changed by a 
specified increment, and this entire procedure is then repeated until 
is minimized. Each iteration requires approximately 15 seconds 
on an IBM 7094 computer. 

Lomax and Inouye (1964) have experimented using the above 
procedure with a polynomial ratio for X(r) • After extensive 


The body coordinates for Mach numbers less than 0.5 are 
ignored because they generally lie on a circular arc and because of 
the scatter introduced by small values of the denominator x — A 
in the first term of equation (7-20). 
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investigation they have discovered the rather remarkable fact that 
ellipsoidal bodies are produced to a high degree of accuracy for values 
of the body bluntness parameter up to 2.25 by the simple, one-para- 
meter equation 
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(7-21) 


where As^ is the shock-wave parameter (illustrated in figure 7-9) . 
This equation is valid for free -stream Mach number greater than 5. 

The denominator provides that for large distances from the nose, the 
shock shape is twice the Mach line slope. This condition was imposed 
to avoid essential instabilities although it does not necessarily pre- 


vent the shock angle from becoming 
intermediary region. The use of 
the one-parameter equation '(7-21) 
greatly simplifies the body- 
optimization procedure. In 
our planetary applications, the 
body parameter is always 

less than one and the corres- 
ponding shock-wave parameter 
A5 ■ is between 0.2 and 0.7 
as illustrated in figure 7-9. 

For these values equation (7-21) 
yields body shapes that match the 


smaller than the Mach angle in some 



magnetopause and the ionopause to 
a very high degree of accuracy. 
The typical value for the total 
absolute error in the body 

coordinates is less than 0.1% . 


Figure 7-9. Shock-wave parameter 
A5 for ellipsoidal 
bodies (Lomax and 
Inouye, 1964). 


Values of A_j which Lomax and Inouye (1964) found for ellipsoids of 
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varying bluntness and for Mach numbers = 6 and 10 are illustrated 

in figure 7—9. 

For our purposes the inverse method just outlined satisfac- 
torily computes the solution for the flow field in the subsonic -transonic 
region. The computed flow properties are interpolated along a line 
joining the shock and body in the supersonic region. These data are 
then used as input for the method of characteristics program to con- 
tinue the calculation do’imstream. 


2 . Method for supersonic region 

A computer program developed by Inouye ^ (1965) at 

Ames Research Center and based on the method of characteristics is 
used to determine the flow field in the supersonic region. This program 
is comprised of a main program and 33 subroutines, for which Inouye 
ejt (1965) have described and also displayed flow charts. Since 

this method is both analytically and numerically stable for supersonic 
flov7 and since specific accounts of the method of characteristics are 
given in numerous textbooks (see,e.g., Courant and Hilbert [1962], 
Garabedian [1964], Hayes and Probstein [1966]), we do not present a de- 
tailed discussion of this method, but rather only outline the general 
procedure for the calculations . 

Along the characteristic or Mach lines defined by = 

tc3n(0— /-^) 3 the first three partial differential equations of (7-6) 

reduce to the following ordinary differential equations (Hayes and 
Probstein, 1966) 


cofju 


dp 


-!- de 


sin 0 dr 


(7 -22a) 


on the left -running characteristic line 


dr 

dx 
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on the right -running characteristic line n ( 6 ~/-i) where 

0 denotes the angle between the streamline and the free-stream 

direction X , and M “ SIH*— r is the local Mach angle. Equations 

' M 

(7-22) are solved in conjunction with the energy equation in integrated 
form (see chapter IX for derivation) 


h + - constant (7-23) 

the conservation of entropy along streamlines 

S = SW (7-24) 


and the equation of state in the form 


h = h{p,S) 


(7-25) 


Briefly, the method con- 
sists of starting with flov7 
properties along a non-char- 
acteristic line between the 
body and the shock wave, as 
determined from the solution 
for the subsonic-transonic 
region, and then integrating 
the equations downstream 
along the Mach lines. The 
stepwise procedure is illus- 
trated 'by the typical char- 
acteristic mesh sho’tm in 
figure 7-10. Beginning with 
Icnown data on the starting 


Shock 


Sterling 
line 



Current dale is stored 
in orroy, P(J,K), 
v/here 

J - identities flow 
variables (x,y,V,-- etc) 

K - identifies the iocotion 
of the field points 
(numbered os shown) 

Direction of colculotion 


Typicol choracteristic mesh 


Figure 7-10. The characteristic mesh 
used in the supersonic 
region (Inouye et al. , 
1965). 
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line the calculation proceeds to the body along a right -running 
characteristic 5 and then back to the next starting point or shock 
point as the calculati_ons advance doxjnstream (see Inouye et al . 

[1965] for a flow chart which illustrates this part of the program 
logic.) In figure 7-10 the previously calculated (or input) data 
points are identified by small circles , and the point currently being 
calculated is identified by the shaded symbol. Only the numbered 
points are available in computer memory at this time since the re- 
maining circled points have been written out previously. The stored 
data points are contained in a two-dimensional array, P ( J,K ), 
in which the index J identifies the various flow variables, and the 
index K identifies the location of the point. 

for the calculation of an interior mesh point C in figure 
7-10, the flow properties at three adjacent points are required. 

These points are labeled A, D, and B in. figure 7-10, and cor- 
respond, in the example sho^m, to the points K = 2, 3, and 4 in the 
P array. The calculation of data at the new point is .carried out 
with the use of equations (7-22) through (7-25) and a standard pre- 
dictor-corrector procedure which averages the coefficients of the 
differentials. The procedure is started vjith a crude predictor (i.e., 
that conditions at C equal those at B ) and is followed, therefore, 
by at. least two correctors. This is in contrast to the method used 
for the sub sonic -transonic solution which makes use of only one cor- 
rector, but which uses a second-order predictor. 

For- the calculation of a mesh point located on the shock wave, 
the flow properties at. only two mesh points, the adjacent upstream 
shock point and the nearest interior point, are required because the 
upstream segment of the right -running characteristic line does not 
exist at this point due to the presence of the shock wave. The cal- 
culation of data at the new shock point proceeds with the use of 
equations (7 -22a) through (7-25) along with the shock discontinuity 
equations (7-2) which ha^e replaced equation (7-22b). The procedure is 
then started and refined in the same manner as for an interior point. 
In this way the shock shape and location in the supersonic region is 
calculated as part of the solution ‘and is no longer determined by 
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equation (7~21) as was the case in the nose region. 

For the calculation of a mesh point located on the body, 
again the flow properties at only two mesh points, the adjacent up- 
stream body point and the nearest interior point, are required be- 
cause the upstream segment of the left -running characteristic line 
does not exist at this point due to the presence of the body. For 
this case, the calculation of data at the new body point is carried 
out' with the use of equations (7-22b) through (7-25) along with the 
specification of the flow angle 0 which is determined for that poin_ 
from the known shape of the planetary boundary around which the flow 
is being computed. This procedure is also started and refined in the 
same manner as was done for an interior point. This method thus allows 
the direct use of the previously calculated magnetopause and ionopause 
shapes and contrasts with the inverse method which requires that these 
boundary shapes be fitted with an analytic equation for an ellipsoid. 
Consequently, this method computes the shock shape and flow properties 
for boundary shapes exactly as they dre displayed in figure 6-6. 

In calculating rotational supersonic flow by the method of 
characteristics, Inouye ^ (1965) found it convenient to introduce 

entropy as a flow property since it remains constant on streamlines. 
Their particular version of this method which we used assumes that the 
entropy varies quadratically between streamlines. To illustrate this 
procedure, consider four points in the flow field (see figure 7-10), 

A, D, and B where the flow properties are known and C where they 
are to be determined. The entropy at C can be calculated using the 
flow propertie's at A, D, and B along with the assumption that the 
entropy varies quadratically along the normal to the streamlines be- 
tween A and B . This assumption is valid provided that the higher- 
order derivatives of entropy betoeen A and B remain small. Al- 
though serious errors may occur in the flow-field calculations if this 
condition is not realized, Inouye ^ (1965) have found that these 

errors are the order of the cube of the mesh size and, for not-too- 
coarse meshes, that this scheme computes accurate solutions. 

As in the case of the inverse method, a check on the overall 
accuracy of the method of characteristics is based on physical rather 
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than mathematical arguments. In the subsonic -transonic region the 
conservation of total enthalpy ever 3 rt^here and entropy along stream- 
lines were adopted as independent checks on the physical consistency 
of the solution. However, for the' present method these two require- 
ments are used in the numerical calculations and therefore are not 
independent and cannot be adopted as checks on the overall accuracy. 
For the method of characteristics, instead, Inouye ^ (1965) 

calculate the mass -flow balance between the shock and the body at 
several positions along the body. This calculation is independent of 
the numerical computation of the flow field and can be adopted as an 
independent check on the overall accuracy. For the results presented 
in this thesis, the mass -flow balance xi'as constant to an accuracy! 
consistent with the calculations throughout the entire supersonic 
flow region. This check helps to insure us that our computed results 
do correctly represent solutions to equations (7-6) and (7-7). 
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CHAPTER, VIII 

FEATURES OF THE APPROXIMATE SOLUTION 

A . Results of the Computed Solution 

With the specification of a value for the free-stream Mach 
number , the location of the bow wave and the properties 

of the flow field can be computed for either the magnetopause or any 
of the ionopause shapes shovm in figure 6-6 by application of the two 
highly developed computer programs described in the preceding chapter. 
Although we have only carried out these calculations for solar -wind 
flow past Mars and Venus, we present some of the calculated results 
for flow past the Earth achieved "by Spreiter ^ aJ. (1966a, 1968) and 
Spreitef and Alksne (1969, 1970) in order to provide a background for 
a later discussion about the similarities and differences between 
these two cases. 

1. Flow past the Earth 

Figure 8-1 shows the position and shape of the bow shock wave 
generated by the solar wind's interaction with the magnetosphere 
boundary and calculated by Spreiter et al. (1966a) using .the computer 
programs just described. These results, which' are presented in terms 
of cylindrical coordinates X and p . normalized by the distance D 
as originally denoted by Spreiter et £l. (1966a) but exactly equivalent 
to our Pq , are for = S/3 and Mco= 8 , a representative value 

for the free-stream Mach number Mqq . Also included on this figure 
are several additional solid lines representing streamlines, and 
dashed lines representing characteristic or Mach lines of the 'flow 
which correspond to standing compression or expansion waves of in- 
finitesimal amplitude. As indicated by equation (5-18) with A = 0 , 
Mach lines cross the streamlines at such angles that the local velo- 
city component normal to the wave is always exactly equal to the local 
speed of sound . Mach lines thus exist only where the flow is 
supersonic and, consequently, are absent from the vicinity of the 
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l.O 0 -1.0 

X/D 


magnetopause nose because the 
flow there is subsonic. This 
plot also shoxi7S that the angle 
between the shock wave and the 
incident stream is much larger 
along the entire length of 
the shock wave included in 
the illustration than the asymp- 
totic angle which is just the 
angle between the streamlines 
and Mach lines of the incident 
undisturbed flow upstream from 
the box'? shock \<?ave. We can 


thus see that knox^ ledge of the 

Figure 8-1. Streamlines and X'jave j- i 

° ^ . asymptotic direction of xjeak 

patterns for supersonic 

flow past the magneto- shock waves is of limited use- 

sphere. M £V) “ S 3 fulness in the estimation of 

^ — 5/3 (Spreiter et 

al. a 1968) o the location of the box^ shock 

V7ave, except at extremely 
great distances from the planet , 
Contour maps shoxjing lines of constant density ratio P ! 3 
velocity ratio yl^co > temperature ratio T/Tq^ x-?hich x<rere ob- 
tained from the computed solution of equation (7-1) and (7-2) and 
originally presented by Spreiter ^ (1966a) are shoxra in figure 

8-2 for the same conditions of ^^ 0 = 8 and ^ — 5/3 as in figure 8-1. 

These results shox\r that the density ratio P / remains near the 

maximum value (ZT+l) (jf — i ) = -4 a strong shock X'?ave in a gas 

with = 5/3 along nearly the entire length of the portion of the 
box-? X7ave shox'?n. The gas undergoes a small additional compression as 
it approaches the stagnation point at the magnetosphere nose and then 
expands to less than free-stream density as it floX'?s around the flank 
of the magnetosphere. The velocity remains less than in the free 
stream,, hox-zever, throughout the same region. As indicated, a single 
set of contour lines, although X'?ith different labels, serves for both 
V/V^ and T/T^ because, as Spreiter et (1966a) pointed out. 



93 



Figure 8-2. Density, velocity, and temperature fields for super- 
sonic flow past the magnetosphere. 8 , ^ = 5/3 

(Spreiter ^ , 1966a) . 

the temperature ratio T/T^ is closely related to the velocity 
ratio through the expression 

X = 

Tcx> 2 

derived by integrating equation (5-3) with B = 0 , combining with 

relations for enthalpy h and speed of sound c3 given in equation 

(5-2), and rearranging. If values for (1 and' are given, it 

is then a straightforward and simple calculation to determine T/Tco 

as a function of . The results presented in figure 8-2 show 

that the temperature is substantially higher than in the free stream 

throughout the entire region illustrated. Particularly noteworthy is 

the large increase in temperature of the solar wind in the region of 

the magnetosphere nose. If, for example, the temperature of the in- 

5 

cident solar wind is 10 °K , then the temperature at the stagnation 
point is calculated to be 2.23 x 10^ which is of tie same order of 
magnitude as the temperature of the gas in the solar corona before 
it is accelerated to the high velocities characteristic of the solar 




(8-1) 
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wind. This result simply demonstrates that the energy content of 
the solar wind is, to a large extent, dependent on the temperature 
of the solar corona. Although the gas cools considerably as it flows 
around the magnetosphere, it is still about one million degrees K as 
it passes the Earth ( x/D = 0 ) and about three -'quarters of that 
farther do^mstream at x/D — ! 

Spreiter ^ a^. (1966a) also showed that the constant -velo- 
city contours of figure 8-2 can also be used to represent, with a 
suitable relabeling, lines of constant local Mach numbers M = v/c5l 
Their conclusion follows from equation (5-2) and (8-1) which show that 
a is proportional to T for a given pair of values for H 'and 
rn and that T/T<d depends only on V/Vco for given Mffi> and Jf 
There exists, therefore, for flow characterized by given Mco , JT 
and hn , a unique value for M ‘associated with every value for 
V/Voo . Contours for constant local Mach number are thus identical 
to those for constant -velocity ratio. Similarly, since the mean velo- 
city of the particles is proportional, but not equal, to the speed of 
sound according to the simple kinetic theory of a gas, contours of 
constant ratio of directed to random velocity are also identical to 
those fox* constant V/Vco 

Figure 6-1 and the accompanying discussion show that the value 
of 8 for the free- stream Mach number is well centered in the range 
of value's to be expected in the incident solar wind. However, since 
both the velocity and the temperature,, and hence the speed of sound, 
of the solar wind vary substantially in the course of time, it is of 
interest to review some of the results given by Spreiter ^ ad. (1966a, 
1968) for other Mach numbers. Figure 8-3 shows their plots of con- 
tours of constant -density ratio free stream Mach numbers 

of 5 and 12 . Although differences in the density ratio are clearly 
evident, they remain sufficiently small to be of only secondary impor- 
tance for most purposes. It should be noticed, however, that the bow 
wave recedes from the magnetosphere as the Mach number diminishes al- 
though the change is small as the Mach number decreases from 12 to ‘ 8 
The entire portion of the bow wave shorn for M=12 is, j.n fact, very 
near its. asymptotic position for infinite Mach number, and further 
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M=5 M=I2 

P^P<o 



Figure 8-3. Effect- of Mach number on density field for flow past the 
magnetosphere, ^ = 5/3 (Spreiter £t , 1968). 

increases in Mach number produce virtually no effect. This lack of 
dependence on Mach number does not apply, however, far do’wnstream of 
the planet where the bow wave approaches alignment with the asymptotic 
direction of weak discontinuities in the undisturbed incident solar 
wind . 

The corresponding results for the velocity and temperature are 
shown in figure 8-4. Like the density, the velocity field is hardly 
affected by the change from Mach number 5 to 12 . However , the 
temperature fields for these two Mach numbers differ substantially, 
the temperatures for Mach number 12 being about five times larger 
than for Mach number 5 . Such strong effects on the temperature can 

be easily explained by an examination of equation (8-1) for T/Tco 
A lthough we do not display the results graphically here, Spreiter ^ al . 
(1968) showed that the temperature contours can- also be put into a 
form relatively independent of Mco normalizing the temperature 
^y Vto tather than Tk, . The explicit form they chose was the 
ratio 2 Cp (T-Tca)/v CO of tlie change of the enthalpy Cp (X ” Toq) — 
^Cv(T “Too) kinetic energy and derived the expression 

2Cp(T—T„)/v<| = I - vVv CO by combining the equation 
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M = 5 V/Veo 



M = I2 



Figure 8-4. Effect of Mach number on velocity and temperature fields 
for flow past the magnetosphere, ^' = 5/3 (Spreiter 
et al . , 1968) . 

(T~Ta 3 )/Too = , the relation h=CpT defining the 

enthalpy per unit mass, and the expression Mco = Vco/{ 2 (fRTco) 
for the free-stream Mach number of fully ionized hydrogen. Since 
vVvl, tends to be independent of Mach number at large M^o , as 
illustrated in figure 8-4 for Mco “ 5 and 12, and Tod in the same 
region is very much smaller than T , they concluded that the ratio 
2Cp(T-Tco)/v|, also tends to be invariant with changes in Mach 
number at large Mg, . 

a. Distortion of the interplanetary field . The results 
presented in the preceding paragraphs represent solutions of equations 
(7-1) and (7-2). As described previously, once the dynamic flow 
properties- are numerically computed from these gasd}mamic equations 
and the orientation for the magnetic field in the incident stream is 
specified, the deformation of the interplanetary magnetic field B by 
the flow around the planet can be determined by numerically integrat- 
ing equations (7-3), or, equivalently, equations (7-4) which indicate 
that magnetic -fie Id lines deform in exactly the same manner as do 
fluid lines. Spreiter ^ al. (1966a) have demonstrated that these 
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equations lead to a straightforx^ard , but tedious calculation in which 
the vector distance from each point on an arbitrarily selected field 
line to its corresponding point on an adjacent field line in the dox^n- 
stream direction is determined by numerically integrating ^ ^ d T 
over a fixed time interval X • This procedure in general results 
in field lines that are curved in space. 

The intensity of the magnetic field at any point may then be 
determined from the relation |^||jBco i = {p[Aii)/(poo|AUl) 

in which aJJ is the length of a small element of a flux tube. 

Figure 8-5 shox^s the results of such calculations performed for three 
different angles betx^een Bco and Vco by Alksne (1967) and 



Figure 8-5. Direction and intensity of the magnetic field in the plane 
of S 3 ?mmetry for three different angles betxi7een Bco and 
Voo 5 Mco=^ 8 , ^ = 5/3 (Spreiter and Alksne, 1969), 
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Spreiter and Alksne (1969) for the plane that passes through the cen- 
ter of the Earth and contains- both the Voo ^^^d Bco vectors. Be- 
cause of the assumed axial symmetry of the magnetosphere and sur- 
rounding flow, all these field lines are planar, rather than spatial, 
curves. While the calculation of the results for the case shoxim in 
the upper left must be carried out in more or less the method just 
described, the results for the other cases displayed in figure 8-5 can 
be obtained more easily. Those in the lox-^er left for Bco parallel 
to Vco can be obtained immediately from the gasdynamic calculations 
because Imai (1960) found that B is proportional to throughout 

the entire flox-7 field for this case. The field lines thus coincide 
with the streamlines, and the contours of constant B coincide x-7ith 
those of constant . The results shoxm on the right for 45 ° 

angle betx-?een Bco snd Vcb x^ere originally calculated by Alksne 
(1967) in the lengthy manner described above, but recently Alksne 
and Webster (1970) have shoxm that the linearity of the equations 
(7-4) for B 3 x^ith p and V already specified by the axisym-. 
metric gasdynamic calculations, permits the magnetic field at any 
point in the plane of symmetry to be determined by addition of the 
B vector for 0 and 90 degree orientations. Furthermore, re- 
suits for any angle betxx'een Bca Vqd can be obtained by de- 

composing Bco iato components (Bco}[j and {Ba>)j_ parallel and 
perpendicular to Voo , and superposing the contributions. 

Figure 8-6 (a) shows the magnetic -fie Id direction and inten- 
sity for the more complex case in x^hich the set of field lines con- 
sidered in the incident solar wind is contained in a plane parallel 
to the plane of symmetry but offset from it by a distance O/S . 

The lines of intersection of the offset plane xvith the box-/ X7ave and 
the magnetosphere boundary are indicated by the solid and dotted 
lines. Since the field lines are no longer confined to a plane behind 
the boxj wave, the results are presented in the form of a pair of pro- 
jections, These and related results for other cases x-7ere originally 
calculated by Alksne (1967) and Spreiter at aT. (1968) in the laborious 
manner described at the beginning of the preceding paragraph. Alksne 
and Webster (1970) have shoxm that' the magnetic field ac any point in 
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(a) Magnetic -field direction and intensity, Mco = 8 . 



(b) Method of calculation. 


Figure 8-6. Direction and intensity of the magnetic field out of the 

plane of symmetry, and illustration of components of inter 
planetary magnetic field used in 'equation (8-2) 

.(Spreiter and Alksne, 1970). 
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the flow field can be calculated more simply, hov?ever, by vectorially 
summing the contributions of the three component fields indicated 
in figure 8-6(b). At any point P , the magnetic field B c. is thus 
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(8~2) 


In this equation, the three ratios with subscripts 11 , _L , and 

h refer to the contributions associated with the components of 'Bco 
parallel to ^co > perpendicular to Vco in. the plane that contains 
the point P , the center of the Earth, and the vector Vco ; and 
the normal to the that plane, From the previous discussion, we see 
that ( Bp/Boa) II = (/3p yp)/(^^lVcol) and (Bp/B<o)j_ is as illus- 
trated in the upper left part of figure 8-5 for = 8 and 

5/3 . Since the remaining ratio (Bp/Bco)|^ can be shown equal 
to {rp/7p}/{ Ppo, ^£o) where r is the radial cylindrical coordinate 
of the streamline, as indicated in the sketch of figure 8-5 (b), 
values for it. can be determined directly from the gasdynamic solution 
for the flow. 

The results displayed in figure 8-5 clearly show how the 
magnetic -fie Id lines bend discontinuously as they pass through the 
bow wave at any angle except a right angle, and then curve in a con- 
tinuous manner throughout the entire region between the bow wave and 
the magnetosphere. The discontinuous bend at the shock wave is, 
moreover, always in the direction that preserves the sign of the 
tangential component of the field, as is required for all physically 
relevant hydromagnetic shock waves. Interestingly, the field lines 
illustrated in figure 8-5 are all draped around the nose of the mag- 
netosphere. This is characteristic of the results for all relative 
orientations, except perfect alignment, of the magnetic -field and 
velocity vectors for the special plane for which these results are 
presented. Outside of this plane, however, the field lines drift 
around the nose with the flow and deform as illustrated in figure 
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8-6(a) into three-dimensional curves. 

' The results shorn in figures 8-5 and 8-6 indicate that the 
intensity of the magnetic field is usually a few times larger near 
the magnetosphere nose than in interplanetary space, except when the 
magnetic field is aligned or nearly aligned with the flow direction, 
and that it may be either larger or smaller than the incident -stream 
value along the flanks of the magnetosphere. They also show that 
the direction of the field in the region between the bow wave and the 
magnetosphere is highly distorted by the flow and bears little 
direct relation to the direction of the field in the incident solar 
wind . 

The results in these figures are for a free- stream Mach 
number of 8 , but no qualitative changes are anticipated if a 
somewhat different value is selected since the velocity field is 
only slightly influenced by variation of Mach number. 

2* Flow past Mars and Venus 

In exactly the same way as was done for the Earth, once a 
value for the free-stream Mach number is specified, the location of 
the bow wave and the properties of the flow field can be computed for 
any of the Martian or Venusian ionopause shapes shoim in figure 6-6 
by application of the numerical gasdynamic programs. 

Results for the location of the bow. wave for flow with 
Mco— 8 aiiti 5/3 are presented in figure 8-7 for each value for 
H/r^ for which the coordinates of the ionopause are illustrated in 
figure 6-6. Although these results are for a specific Mco 5 it has 
been shown in the previous section relating to the Earth that the 
shape of the bow wave, and also the associated distributions of 
P j 5 s and B/Bq:> for any given direction for Bco are 

relatively independent of free-stream Mach number for all M^o 
greater than about 5 (Spreiter ^ aj. , 1966a, 1968). 

In figure 8-8 are showm contour maps for the density ratio 
PI Pa> > velocity ratio \j/\/co > temperature ratio T/Tqj for . 
flow with H/r^-0.2 . Such a value is appropriate, for example, 
for a Venusian ionosphere composed primarily of atomic hydrogen at 
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Figure 8-7. Calculated location of ionopause for various H/r^ , 
and associated ' location of bow wave for Mco~8 5 

5/3 • ■ The corresponding results for the shape of the 
Earth's magnetopause and bow wave, nondimensionalized so 
that the magnetopause nose is at x/ri> = I , is included 
for purposes of comparison. The dashed line for the bow 
wave for H/r^ = 0.2 has been omitted because it is in- 
distinguishable from the dotted line representing the 
Earth's bow wave (Spreiter, Summers, and Rizzi, 1970b). 
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about 600°K, 

In view of the current tendency for some to interpret the 
Mariner~5 results as indicating that the primary constituent in the 
Venusian atmosphere is molecular hydrogen, or even heavier ions, we 
have also made a set of calculations for the density, velocity, and 
temperature in the flow field for H/r^=0.1 . The results are 

shorn in figure 8-9. Corresponding plots for other quantities of 
interest, such as the magnetic- field strength and direction, or 
-MaMrellian pro ton- velocity distribution, can also be calculated in 
the same way as has been done previously for the Earth by Spreiter 
and Alksne (1969, 1970) and Spreiter ^ (1966b). Such calcula- 

tions’ have not been carried through, hox-?ever. This is partly because 
they are time consuming to do and require considerable space to dis- 
play the results for even a minimum number of cases, but principally 
because an approximate correspondence rule to be described in the next 
section makes it possible to convert quickly any of the numerous con- 
tour plots already available for the properties of the flow around 
the Earth's magnetosphere into that for the flow around an ionosphere 
having any H/Pq between 0.01 and 1 by a simple relabeling of 
the coordinate axes . 

a. Correspondence rule . The close relationship between ' 
the shape of the magnetopause and that of the ionopause for VljCo ~ 
0*2 , and the general similarity of the ionopause shapes for all 

H/rjj , suggests the possibility that a correspondence rule relating 
ionopause shapes for H/Pq other than 0.20 to an appropriately 
scaled magnetopause may be found if the coincidence between the Earth 
and planetary centers be relinquished. Great practical utility would 
result from the availability of such a correspondence rule, because 
it would enable a substantial body of results already calculated for 
solar-wind flow past the Earth to be applied with minor change to 
nonmagnetic planets having a wide range of ionospheric parameters. 

The results displayed in figure 8-10 show that it is indeed 
possible to achieve a fit that is probably sufficiently good for most 
purposes relating to the interpretation of data obtained in space. 
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Figure 8-8. Density, velocity, ’and temperature fields for supersonic flow past 
the ionosphere; H/fo = 0.2 » Mqo— 8 > ~ 5/3 (Spreiter, 

Summers, and Rizzi, 1970b). 
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H/ro = 0.! Mco = 8 /=5/3 



Figure 8-9. Density, velocity, and temperature fields for supersonic flow past 
the ionosphere; HlCo ~ O.i , Mqj =8 , ^ ~ 5/3 (Spreiter, 

Summers, and Rizzi, 1970b). 
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This figure is basically that given originally by Spreiter al ■ 
(1966a) for the magnetosphere boundary and shock wave for the Earth, 
and the associated characteristic or Mach line pattern, for Mcx)=8 
and g" = 5/3 . Superposed on it are the ionopause curves from figure 

8-7 each with its nose retained at x/Qj = ! , 0=^0 , but with its 
scale adjusted so that its ordinate coincides with that for the mag- 
netopause curve at © = TT/2 . The center of the nonmagnetic planet 
is, in general, no longer at the origin of coordinates, but at the 
points >^c/n3 indicated along the axis of S3mimetry for each value 
for H/r^) • The coordinates of these points are listed on figure 8-10 
for accuracy and convenience in applications of the correspondence 
rule. To convert the nondimensional shape for the magnetopause into 
that for the ionopause for a given H/r© > we must thus place the 
center of the planet at ^c/Qj change the labeling of the scales 

so that X/r^ —0 at Xc/f^ , and x/Tq =1 at the -ionopause nose. 

It follows from the close correspondence between the coor- 
dinates of the ionopause 'and the magnetopause that the coordinates of 
the bow wave should also 'display a similar relationship. It may be 
seen from figure 8-10 that this is indeed true. A similar degree of 
correspondence may be anticipated for the coordinates of contour 
lines for constant values of the flow parameters, such as those for 
the density, velocity, and temperature shoxm in figures 8-8 and 8-9. 

Although some of the ionopause curves depart • significantly 
from the curve for the magnetopause somewhat do’cm.stream of the planet, 
the effects of these differences have no influence upstream of the 
rearwardly inclined characteristic line emanating from the point on 
the boundary vzhere the differences first become significant. It 
may be seen from figure 8-10 that this is sufficiently far downstream 
for ionospheres having .H/f^ betxreen 0.01 and 1 that nearly all 
of the results we have presented previously for the properties of the 
flow field about the magnetosphere can be carried over V7ith no other 
change than relabeling the scales to obtain a good approximation for 
the conditions' around Mars and Venus. 
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Figure 8-10. Illustration of degree of coincidence of 

curves representing the ionopause and bow waves 
for various H/r^ obtained by application of 
correspondence rule. The solid lines have been 
omitted where they are indistinguishable from 
the dotted lines (Spreiter, Summers, and 
Rizzi, 1970b). 
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B. Evaluation of the Theory 

Since all of the foregoing theoretical results are founded 
on the numerical solution of the gasd 3 mamic equations which, we argued, 
approximate the hydromagnetic equations in this application, it is 
informative and important to examine these theoretical calculations 
for internal consistency. This examination should assay the three 
principal postulates of our theory; namely, do the numerical calcula- 
tions represent accurate solutions of the gasdynamic equations, how 
precisely do these equations approximate the complete hydromagnetic 
equations, and finally, how accurately does the Newtonian pressure 
predict the flow pressure on the boundary. 

The accuracy of the numerical solution has been verified 
theoretically by the consistency checks built into the computer pro- 
grams, As discussed in chapter VII, in the subsonic region the con- 
servation of entropy along .streamlines and total enthalpy ever5u-7here 
provides an independent theoretical check on the pr-ograms internal 
consistency whereas the constancy of mass flow furnishes a similar 
check in the supersonic region. All of the numerical calculations 
presented in this report accurately satisfies these checks. In addi- 
tion, Inouye and Lomax (1962), Spreiter* aX, (1966a), and others 
have found that these computer programs provide solutions to the gas-' 
dynamic equations that do compare very well with actual aerodynamic 
measurements obtained experimentially . On this question then, we are 
assured that out numerical results do indeed accurately represent the 
solution of the gasdynamic equations. 

One of the simplifications fundamental to the attainment of 
the numerical results presented in the foregoing sections of this 
chapter assumes that the Alfven Mach number is sufficiently 

greater than unity so that the gasdynamic equations satisfactorily 
approximate those of hydromagnetic theory, and subsequently that the 
magnetic field can be calculated by application of the concept of the 
magnetic -fie Id lines moving with the fluid. Comparison of the values 
for the Alfven speed displayed in figure 6-1 with representative values 
for the velocity of the solar wind assures that this condition is 
easily met in the free-stream. A similar comparison for the conditions 
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behind the bow wave is more difficult to make a priori because both 
the flow velocity and the Alfven velocity vary considerably throughout 
the flow field. Once the velocity, density, and magnetic -field in- 
tensity have been calculated for any point, however, the Alfven Mach 
number is easily evaluated, and we can therefore determine to what 
extent the postulate of large is fulfilled. The findings of 

such an examination carried out by Spreiter et (1968) are presented 

in figure 8-11 for two orientations of the interplanetary magnetic 
field both with Moo=8 • The results are given in the form of a 
ratio of the local Alfven Mach number to the free-stream Alfven 



figure 8-11. Alfven Mach number distribution, Mco = 8 and 
^=5/3 (Spreiter et al. , 1968). 


Mach number ■ in order to take advantage of the simple propor- 

tionality betxjeen the two quantities at any point in the flow field. 
IT usual for the incident solar wind, the results 

shov7 that is greater than unity every^-7here except in a small 
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region near the magnetosphere nose. Although, as predicted by Lees 
(1964), some localized departures from the calculated results might 
be anticipated in this region, it would seem that the larger-scale 
features of the flow should be quite adequately represented by the 
gasdynaraic and frozen-field calculations. Moreover, although this 
case is not displayed here, as the magnetic-field orientation approaches 
alignment with the flow direction, this region shrinks even smaller, 
and the validity of the numerical results becomes even more assured. 

As decreases below 10 , however, the critical regions expand 

significantly for all orientations, and increasing differences be- 
tween the approximate numerical results and more exact theoretical 
or observational results should be anticipated. When M^oo becomes 
as small as 2 or 3 , vjhich it does only infrequently, the 

critical region of small local Alfven Mach number would become so ex- 
tensive that substantial departures from the results so far presented 
must be expected. (This ’interesting behavior at low Alfven Mach num- 
ber will later be more thoroughly discussed in chapter X.) 

Lastly, there remains the question of how accurately does the 
simple Newtonian pressure formula p= approximate the 

variation of pressure along the magnetosphere boundary. The answer 
to this question for the present applications is easily arrived at 
by comparing the Newtonian pressure on the magnetopause with the pres- 
sure derived from the detailed gasdynamic calculations. The results 
of such a comparison, originally presented by Spreiter et (1966a), 
are sho™ in figure 8-12 for several values of Mc» ^ 

can easily see that the simple Newtonian expression does indeed 
provide a generally good approximation over most of the magnetosphere 
boundary of interest in the present studies. The agreement is par- 
ticularly good over the portion of the magnetosphere along which the 
flow is subsonic. No significant changes are expected, therefore, 

■in the shape of the nose portion of the bow wave if a more accurate 
and necessarily more complicated calculation of the pressure is car- 
.ried out. The Newtonian pressure formula does, however, underes- 
timate the pressure somewhat along the flanks of the magnetosphere 
which indicates that a more exact calculation using the gasdynamic 
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Figure 8~12. Comparison of exact and approximate 

pressure distributions on magnetosphere 
boundary (Spreiter ^ » 1968) . 

pressure distribution would yield a slightly slimmer shape . This 
change in shape along the flanks of the magnetopause should never*- 
theless be small because the geomagnetic pressure. v?hich 

balances the gas pressure p is nearly proportional to the inverse 
sixth power of the distance from the center of the Earth, A point of 
equal interest displayed by these results is the lack of significant 
variation of the pressure distribution with changes in and 

over the range of cases presented. 

Although we do not display the results here, we have per- 
formed a corresponding examination of the accuracy of the modified 
Newtonian pressure formula 


P = p^~ Bysir)cos^y -i- P^ 


for the ionopause of either Mars or Venus and found the results similar 
to those- presented in figure 8-12. For this case, too, agreement is 
particularly good in the nose region, but the modified pressure formula 
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partially underestimates the pressure along the flanks of the iono- 
pause. The resulting slimmer ionopause shape which a more exact cal- 
culation using the gasd 3 mamic pressure distribution would indicate, 
however, would also not differ very much from the present shapes be- 
cause the ionosphere pressure which balances the gas pressure p varies 
exponentially with distance from the center of both Mars and Venus. 

C . Comparison with Observations in Space 

Although a number of checks on the internal consistency of the 
fluid theory and its subsequent results are given in the, original 
references, some of which we have reviewed in the previous section, 
the final test of any theory of solar -wind flow past these planets 
must be based on comparisons v?ith observations actually made in space . 
Such comparisons have been made by several different authors on several 
levels of detail. These range from simple comparisons of gross fea- 
tures, such as. the existence and location of the magnetosphere boundary 
and bow wave, to more detailed comparisons of the velocity, density, 
and temperature of the flowing plasma and the intensity and direction 
of the magnetic field. Here, rather than review in detail the enormous 
amount of space observations made in the vicinity of the Earth, we 
shall restrict ourselves to only a brief account of near -Earth data 
intended to show the actual validity of one of the fluid assumptions 
and then concentrate on a more thorough comparison between the appro- 
priate theoretical results presented so far and the data observed 
near Mars by Mariner spacecraft 4, 6, and 7 and near Venus by Mariner 
5. In addition, several extensive summaries of the comparison between 
Earth-related space observations and theoretical fluid calculations 
have been given recently by Spreiter et al. (1968) and Spreiter and 
Alksne (1969, 1970). 

A key point in the hydromagnetic theory of the steady-state 
interaction of the solar wind and both the magnetic field of the 
Earth and ionospheres of Mars and Venus is that the magnetopause must 
be represented by a tangential discontinuity. Sonnerup and Cahill 
(1967, 1968) have investigated this by examining the magnitude of the 
magnetic -field component normal to the magnetopause in the data of 
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Explorer 12 satellite. Although the accuracy of those data is mar- 
ginal for the task, and the orientation of the magnetopause is un- 
certain, they conclude that this quantity is sufficiently small for 
most of the crossings to he compatible with the theoretical require- 
ment given in'equations (5-10) that * Fairfield (1967) has 

also investigated this question using data obtained simultaneously ■ 
by the IMP 1 and IMP 2 satellites when they were on opposite 
sides of the bow wave. Within the accuracy of the determination, the 
results showed that the magnetic field is tangential to the magneto- 
pause, which satisfies the condition Bp,=0 , and also that the 
orientation of the magnetic-field vectors throughout the entire re- 
gion between the bow wave and magnetopause is consistent with the 
theoretical predictions provided by the simplified hydromagnetic 
model. This conclusion has been recently confirmed and strengthened 
by the results of a more extensive study by Behannon and Fairfield 
(1969) using 1661 hours of data from four space.craft. Explorers 28, 

33, 34, and 35. 

1. Comparison with Mariner 5 data for Venus 

Due to insufficient data, a similar examination of the ful- 
fillment of the tangential discontinuity conditions cannot be carried 
out for the ionopause of either Mars or Venus. We can, however, make 
a general assessment of the relation between the theoretical results 
shoim in figures 8-7 to 8-9 and those observed by Mariner 5 as it 
flew past Venus. Figure 8-13 shows the time variation of the bulk 
velocity V » ion number density H , and intensity j8[ and direc- 
tion o( and jS of the magnetic field together with a plot of the 
trajectory giving the position of the spacecraft as a function of 
time (Bridge ejt , 1967). The angles o( and /3 are in spherical 
coordinates in which o{ is the longitude measurement, in the RT 
plane, and from the R direction, where R is a unit vector in the 
antisolar direction and T is an orthogonal vector which is parallel 
to the Sun’s equatorial plane and points in the direction of the 
planet's orbital motion; and /$ is the latitude angle of the mag- ‘ 
netic- field vector, considered positive when northward. For ease of 




TIME FROM ENCOUNTER, min 


Figure 8-13. Plasma and magnetic data measured by Mariner 5 in its encounter 
with Venus on 19 October 1967 (Bridge , 1967). The indi- ■ 

cated trajectory has been projected into a plane by rotation about 
the Sun-Venus line. Times are in minutes relative to the encounter 
time E of closest approach to the planet (Spreiter, Summers, and 
Rizzi, 1970b). 
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reference, lines indicating the theoretical location of the ionopause 
and bow wave for H/ro=0-l and 0.25 , for 'Moo- S , S' = 5/3 , 
and pQ = 6500 km are superposed on the plot of the trajectory. 

The corresponding location of the bow wave for Mco = 5 and H/l^ 0,25 
is also included in order to better assess the consequences of the fact 
that the values of about 590 km/sec and 300,000°K reported for the 
velocity and temperature in the solar wind before and after Mariner 
5's encounter with Venus indicate a value of about 6.5 ' for . 

The magnetometer and plasma probe experimenters have con- 
cluded from the presence' of abrupt and easily recognizable changes in 
their data that Mariner 5 crossed the bow wave at the points labeled 
(D and . These points are reasonably close to the theoret- 
ical locations of the shock wave for H/f^ = 0.25 for Mco ~ 5 , 

and somewhat farther from Venus than indicated by the theoretical 
results for Mco — 8 j for H/pj^cr. O.iO . - On the other hand, the 
theoretical results for H/f^ — 0.25 indicate that Mariner 5 should 
have crossed the ionopause, whereas those for H /Pq = 0.10 indicate 
that it should have just skimmed along this surface. Although the 
data display decreases in n/Doo and V/Voo that are at least qua- 
litatively similar to those for either H/Po^O.lO or 0.25 , 
there is no certain evidence that Mariner 5 actually penetrated the 
ionopause. In particular, the intensity of the magnetic field re- 
mained of the order of that for interplanetary space, and V dimin- 
ished to only a modest fraction of Vcs> rather’ than to a value com- 
parable with the speed of the spacecraft relative to Venus. On the 


other hand, the observations made when the spacecraft was near the 
theoretical location of the ionopause indicate values for n/Poo 
and V/VcD that are only about one-half of those indicated by the 
theory for the flow exterior to the ionosphere. 

A possible explanation for these discrepancies is that the 
tail of the ionosphere might taper inward toward the axis of symmetry 
rather than extend straight downstream. The theoretical results give 
no indication of such a trend, but they are of low reliability for 
this feature of the flow because of the deterioration of the quality 
of the modified Newtonian approximation for the pressure of the solar 
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wind as approaches 90° , and the substantial departures from a 
constant scale height of a planetary atmosphere at great heights 
caused by the diminished gravitational acceleration, and also by 
probable increase in temperature. In addition, it is virtually cer- 
tain that the ionosphere in an extended tail would be supported par- 
tially by dynamical interaction with the flowing solar xjind instead 
of entirely by the simple hydrostatic means assumed in the present 
analysis. Although neither theoretical nor observational evidence is 
as yet sufficient to provide a definitive statement regarding the 
nature of the ionosphere tail, we are inclined to explain the dis- 
crepancy betcreen the theoretical and Mariner-5 results in terms of an 
extended tail with a relati.vely thick boundary layer between the plan- 
etary ionosphere and the solar wind. The latter may be contrasted 
with the thin boundary between the solar wind and the Earth's mag- 
netosphere that is- becoming increasingly evident as the time resolu- 
tion of the observations is improved. That the ionopause might be 
thick well back from the nose, whereas the magnetopause is usually 
thin , is plausible in view of the fundamentally different nature of 
the two boundaries. For the Earth, the magnetopause is essentially a 
boundary betTveen the flowing solar plasma and a relative vacuum, and 
would be somewhat like the boundary formed by a free streamline in a 
water flov? v?ith embedded air cavity. For Venus or Mars, the iono- 
pause is essentially a boundary between tv70 different bodies of plasma. 
As such, a significant boundary layer would be anticipated which would 
provide an increasingly thick transition between the ionosphere and 
the solar wind with increasing distance from the stagnation point at 
the nose of the ionopause. It would spread, moreover, both into the 
outer part of the ionosphere and outv?ard into the surrounding flov 7 , 
broadening the transition region in which all properties of the plasma 
including the density, velocity, and magnetic field change from their 
values for the flowing solar plasma to those of the planetary iono- 
sphere, At the location of the ionosphere boundary indicated by the 
present dissipationless theory, the plasma velocity might be expected 
to be substantially less than indicated by the theory. Since the 
plasma velocities observed by Mariner 5 in the vicinity of the 



theoretical location of the ionopause display such a trend, it is 
tempting to conclude that Mariner 5 entered the boundary layer sep- 
arating the ionosphere and the flowing plasma, but did not enter the 
ionosphere proper. A more detailed examination is clearly required 
before a definitive statement can be made, however. 

2 . Comparison with Mariner 4, 6, and 7 data for Mars 

As described in less detail previously, data from the magneto- 
meter (Smith ^ si- 5 1965), plasma probe (Lazarus ^ a 1967), and 
energetic-particle detectors (Van Allen ^ £l* a 1965 and O' Gallagher 
and Simpson, 1965) on Mariner 4 displayed no effects attributable to 
the presence of Mars as it flew past that planet. Although such a 
negative finding does not provide much material for comparison, 
Spreiter et (1970b) found it at least consistent with the indica- 
tion of the present theory for the conditions measured at the time of 
the encounter that Mariner 4 did not cross the bov? wave into the 
region influenced by Mars at any time when signals xsjere being received 
from the spacecraft . 

More recently, hoxjever, it has been reported by Sm.ith (1969) 
and Kavanagh £t (1970) that a closer examination of the magneto- 

meter data indicates that Mariner 4 may have detected a Martian box^ 
wave, although its effects are weak and its presence could not be con- 
firmed by the other instruments. In figure 8-14 are reproduced 16 hr 
of the Mariner -4 magnetometer record (Smith, 1969) shox^ing the inten- 
sity [b[ of the field as the spacecraft passed Mars. This quantity 
is plotted as a function of the time that the data x^as received at 
Earth, and is 12 min later than the time at which the measurements 
were made and transmitted by the spacecraft. The distance of the 
spacecraft from the center of Mars is indicated -by the aerocentric - 
distance scale at the top of the plots. The point of closest approach 
to Mars is denoted by GA . The txTO abrupt changes in jBj signify 
crossings of the proposed Martian box7 x^ave. Those authors see fit, 
hox'/ever, to repeat the cautioning of the experimenters (Smith ^ al . , 
1965 and Smith, 1969) that the field disturbances seen x-7hile the 
spacecrafu xx>as near Mars could have- been interplanetary -fie Id 
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Figure 8-14. Trajectory of Mariner 4 past Mars, projected into a plane by rotation about the Sun-Mars 
line; calculated location of the ionopause and bow wave for Moo=^6 ? ^ =5/3 , and 

H/Po ■= 0,25 rotated 4.5 to allow for aberration effects arising from Mars’ orbital motion 
as viewed along the trajectory of Mariner 4 (Smith, 1969); and magnetic record from the 
latter source showing abrupt changes at virtually the precise times that Mariner 4 crossed 
the calculated position of the bow wave (Spreiter and Rizzi, 1971). 
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fluctuations bearing no relation to the planet. Although data from 
the plasma detector reported by Lazarus (1967) show a slightly broader 
than average energy spectrum from about 300 to 100 eV when the space- 
craft was inside the proposed bow wave, these indications were con- 
sidered indecisive for the confirmation of the shock wave's presence 
because failure of the instrument's high voltage supply made the lower 
energy data unreliable. If, however, it is assumed that the magnetic 
discontinuities do, in fact, signify crossings of the bow wave, the 
locations of the crossings can be interpreted (Dryer and Heclunan, 

1967) as indicating the existence of a Martian magnetic field having a 

-4 

dipole moment about 2.1 x 10 that of the Earth. 

The alternate interpretation, presented by Spreiter and Rizzi 
(1971) that the Martian bov? wave results from the presence of an 
ionosphere, and may develop even in the absence of a planetary field, 
gains substantial support from the remainder of the results sho\ 7 n in 
figure 8-14. These consist of a plot of the trajectory of Mariner 4 
near Mars projected onto the Sun-Mars-spacecraft plane (Van Allen 
et al . , 1965) and superposed on a plot of the calculated position of 
the bow wave for Moo ~ 8 ^.nd S" = 5/3 resulting from interaction 
with the Indicated ionopause dra^m for ~ 0-95 and H/f^ — 0,25 

The distance Pq was calculated to be 5 percent greater than the 

g 

3.4 X 10 cm radius of Mars from the electron density profiles for 
the Martian ionosphere determined from the Mariner 4, 6 , and 7 ob- 
servations together with the present theory. The value H/po=0.25 
is representative of conditions that would prevail if the dominant 
molecule were atomic hydrogen at a temperature T of about 200°K as 
is indicated by Kliore e^ aT. (1965) and Fjeldbo and Eshleman (1968) 
to be the case in the upper ionosphere, the region of greatest im- 
portance to the interaction with the solar wind. The associated bow 
wave is drawn for Moo = 8 " 5 l>ut results of similar calculations 
for the Earth's bow vjave have sho™ that it would be in virtually 
the same location for higher Mach numbers. Significantly lower Mach 
numbers would lead to more remote locations for the -bow wave, but 
are considered to be unlikely for Mars in view of extensive knowledge 
of solar-wind properties in the vicinity of the Earth's orbit and the 
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prediction of solar-wind theory that Mco increases monotonically 
with increasing distance from the Sun. These theoretical calcula- 
tions displayed in figure 8-14 are identical to the results presented 
by Spreiter ^ al. (1970b) for the same conditions except that (a) 
the direction of the solar wind has been rotated 4.5'* from the Sun- 
Mars line to allow for the aberration effects of the motion of Mars 
about the Sun, in the same way as done previously by Smith (1969) in 
a similar -appearing plot made using a scaled-doxm drawing of the 
Earth's bow wave and magnetopause, and (b) the calculated shock- 
wave location has been indicated for greater distances do^«istream of 
the planet. To facilitate comparison with the magnetometer record, 
markers have been added to the trajectory to indicate the time at 
which measurements at each point were received at Earth. The times 
at which the trajectory crossed the calculated bow wave are indicated 
on the magnetometer record by the vertical dashed lines labeled shock. 
Their near -perfect coincidence with the times of abrupt changes in 
I B I supports the proposition that these data do indeed provide the 
signature of a Martian bow wave and, since Mars almost certainly 
possesses an ionosphere having a ratio H/Cq of the order of .25 
in its upper -levels, that its existence arises from the interaction 
betv7een the solar wind and the planetary ionosphere along the lines 
described by Spreiter £t (1970b). 

' A corollary is that there is no need to postulate a planetary 
magnetic field to account for the observations. The value proposed by 
Dryer and Heckman (1967) may still be regarded as an upper limit for 
the strength of a dipole field, but the actual existence of such a 
Martian field should not be inferred from the Mariner 4 observations. 

Mariners 6 and 7 subsequently approached within 2000 km of 
the Martian surface, but neither spacecraft carried a magnetometer, 
plasma probe, or energetic -particle detectors. As a result, the only 
observation of relevance to the present study is the confirmation of 
the previous determination of the ionospheric density by Mariner 4. 

The data from these experiments, and also those of Mariner 5 
for Venus, are of considerable significance to the present theory, 
hovrever, because they demonstrate that the ionospheres of -both Mars 
and Venus are .sufficiently dense to stop the solar V 7 ind at the subsolar 
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point and to deflect it around the ionosphere . Since the margin 
in the product is only about 6 for Mars and 10 for 

Venus, considering the ionospheric properties to remain fixed, the 
knoxm variability of the solar wind suggests that it may be sufficiently 
enhanced at times to proceed directly through the ionosphere at the sub- 
solar point and into the lower atmosphere where it would be absorbed by 
collisions. Under such conditions, the interaction would be quite 
different from that described here, and would probably tend toward 
that described by Spreiter et (1970a) for the Moon. Since such 
conditions can be expected to occur only occasionally, if at all, 
particularly if the ionospheric density and temperature are enhanced 
by a strong solar-i^ind flow, we conclude that the theory given here is 
both plausible and capable of providing a reasonably accurate descrip- 
tion of the conditions that prevail most of the time at Venus and Mars. 
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CHAPTER IX 

NWIERICAL SOLUTION OF THE EXACT EXTERIOR FLOW 

All of the results presented so far in this thesis have been 
based on solutions of the gasdynamic equations which for large Alfven 
Mach number only approximate the complete magnetohydrodynamic equa- 
tions (5-1) of hydromagnetic theory. We have showed in the preced- 
ing chapter that these gasdynamic equations provide a satisfactory 
approximate solution for most cases, but the decoupling of the dynamic 
and magnetic properties which results from this simplification does 
represent a loss of the magnetic character of the solar -wind inter- 
action, In the gasdynamic simplification not only do the terms rep- 
resenting the magnetic force in the flow equations disappear, but 
even the number of parameters required to characterize the flow is 
reduced. 

For the numerical gasdynamic solution we needed to specify 
the Mach number and velocity direction upstream of the ,bow x«?ave while 
for the full hydromagnetic solution, in addition to those tx70 para- 
meters, xje must also specify the free-stream Alfven Mach number and 
magnetic-field direction. However, because of the combined effects 
of M^«t and the strong interaction nature of the flex'?, the 
dynamic terms in the momentum equation of (5-1) will dominate, and 
therefore, for a value of 10 or above the differences be- 

tween our approximate solution and a more exact one are expected to 
be small.' Nevertheless, as the magnetic field becomes stronger and 
decreases, figure 8-11 shox?s that these differences x?ill grox?. 

At the present time the general solution of equations (5-1) 
through (5-5) representing the fluid description of the solar-x?ind 
interaction x?ith the planets Earth, Mars, and 'Venus has not yet been 
XTOrked out. This solution is of considerable interest because it 
XTOuld express the interdependent dynamic and magnetic interaction and 
consequently be a more consistent fluid theory for representing these 
phenomena. Since equatio'ns (5-1) through (5-5) represent a more 
refined theory, the general solution x?ould also provide a good check 
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on the accuracy and applicability of the simpler gasdynamic one al- 
ready presented. Furthermore, a comparison- between the exact solution 
based on magnetohydrodynamics and the approximate one based on gas- 
dynamics would determine whether Dryer's (1970) suggested use of the 
free-stream magnetoacoustic Mach number in an attempt to include 

the effects of both and in the gasdynamic solution ac- 

tually improves the accuracy of that theory. 

A . Equations for Aligned Hydromagnetic Flow 

In our fluid formulation the hydromagnetic equations (5-1) 
through (5-5) represent the flow of the solar wind past the planets 
Earth, Mars, and Venus. Notv, instead of simplifying these equations 
to those of gasdynamics as we did before, we attempt to determine the 
general solution to them. The completely general solution is, how- 
ever, still intractable at this time, but we are able to find the 
solution for the particular case of aligned velocity and magnetic - 
field vectors. Although this hydromagnetic solution is for a par- 
ticular orientation of the upstream magnetic field and cannot show 
the effect of variation of that direction, it will demonstrate, what 
the approximate gasdynamic solution did not, the role that the mag- 
netic field ^ and the parameter pla-ys in determining the 

character of the flow downstream of the bow wave. Moreover, study of 
the Friedrichs, or phase and .group velocity, diagrams .illustrated in 
figure 5-1 for hydromagnetic flow with any orientation of the flow 
"direction relative to the magnetic -field direction shows that perfect 
alignment of the Vco and vectors is not essential so that 

qualitatively similar effects may be anticipated for other alignments 
as well. This solution, therefore, should be representative of the 
character of the general hydromagnetic solution and provide some in- 
sight to the influence that the magnetic field bears on the flow 
properties . 

Compare the second equation of (7-4) and the expression 
~ rate of change of velocity for a fluid 

particle in steady flow. Imai (1960) originally pointed out that the 
rate of change of the vectors B/^ and V are given by identical 
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equations and concluded that if these two vectors are initially in 
the same direction they xi?ill remain parallel ever 3 ?where and their 
lengths along a given streamline will remain in the same ratio. 
Mathematically, this result is expressed by 



P 


= ;i(T)v 


or 


(9-1) 


B = A(y) 


where is a constant along streamlines but can differ on 

different streamlines. Thus, if the velocity and magnetic -field vec- 
tors once become aligned, then they will remain aligned throughout 
the entire flow field, and the magnetic-field lines will coincide with 
the streamlines of the flow. Looking at it from a more physical view- 
point, we see that this result is just a consequence of the frozen- 
field property of 'these perfectly conducting hydromagnetic flows. 

By following and extending the work of Imai (1960) and 
Spreiter et ,al. (1970a), we can derive the explicit form of the general 
equations for aligned hydromagnetic flow and then specialize them for 
our application. The conditions which A(v) must satisfy are easily 
derived by substituting equations (9-1) into the third of equations 
(5-1) which yields the expressions (V-y) X =0 [a] = 0 

requiring AIv) to be constant along streamlines even passing through 
a shock wave. In our application upstream conditions are uniform, 
and A is therefore .a global constant . For this situation equations 
(5-1) and (5-3) reduce to (Spreiter et aX . , 1970a) 
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= 0 




(v-2)s = 0 


S Sco Cy V/yi 


P/P. 


oo 


W 


(9-2) 


'?• [^y{~v^+h)] = /o(v*y){~v2 -j-h) =0 


and the discontinuity conditions (5-5) to 


[^V„] = 


0 


/oV^V + 


(O 


P + 8ir I ^ ^■ 




= 0 


(9-3) 


/0Vn(-^V^ 4- h) 


-I- h 


= 0 


Equations (9-2) and (9-3) along with the auxiliary relations (5-2) 
and (5-4) are the- equations of nondissipative aligned hydromagnetic 
flow which in our fluid theory represents solar -wind flow with 
parallel velocity and magnetic fields. They shovj that aligned flow 
past any of the boundary shapes for either the Earth, Mars, or Venus 
illustrated in figure 6-6 is axisymmetric about a line through the 
center of the planet and parallel to the direction of the solar vjind. 
This is so because the azimuthal components of all the terms vanish 
identically for flows that are uniform at infinity and disturbed 
only by an obstacle that is axisymmetric about a line parallel to the 
direction of flow in the undisturbed stream. Nonaligned 'flows past 
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these planets are not axisymmetric , hotrever, because the azimuthal 
components of the magnetic terms in equations (5-1) and (5-5) are, 
in general, not zero. In the latter case, nevertheless, mirror 
symmetry still remains about the plane through the center of the planet 
that contains the Vc» and’ Boa vectors. 

Since = 0 indicates that entropy is constant along 

a streamline, except where derivatives cease to exist as at a shock 
wave, and S = Sod upstream of the shock, the flow is isentropic 

until the shockwave is crossed. Equations (9-2) shov; that + n 

is also a constant along a streamline equal to -^Voo ho? = hjt 

where subscript S’!" refers to conditions at a stagnation point any- 
x-jhere in the isentropic part of the flow field. Since equations (9-3) 
show, in addition, that -^V^+h is conserved across a shock wave in 
aligned flow, is an absolute constant throughout the entire 

flow field. We should note that the constancy of hgf everyr'jhere in 
aligned flow also holds in the corresponding gasdynamic flow. This 
equivalence results because the Lorentz force BXCUrl B is every- 
where normal to the streamlines in aligned flow, and hence the energy 
equation (5-3) reduces to that of gasdynamics. Moreover, the con- 
tinuity and energy equations of (9-2) are in fact exactly equivalent 
to those corresponding equations (7-1) of gasdynamics. Equations (9-2) 
actually differ from the gasdynamic equations (7-1) only by the Lorentz- 
force term standing on the right-hand side of the momentum equation. 

A corresponding examination of equations (9-3) and (7-2) reveals a 
similar comparison for the discontinuity relations. 

However, the numerical solution to equations (9-2) and (9-3) 
remains formidable because of their nonlinear nature and complicated 
regions of mixed elliptic, parabolic, and hyperbolic character. Fur- 
thermore , the partial derivatives representing the Lorentz force in 
the momentum equation prevents making direct use of any of the vast 
amount of work already carried out in developing the two gasdynamic 
computer programs described in chapter VII because their finite-dif- 
ference schemes approximate the gasdynamic momentum equation and do 
not correspond to the momentum equation for aligned flow. Conse- 
quently, to compute the solution of equations* (9-2) and (9-3) we would 
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have to begin anei:7 and develop entirely different programs for which 
extensive analyses of the stability and convergence would have to be 
carried out iri much the same way as was done for the gasd 3 naamic pro- 
grams in chapter VII. The magnitude of such a task is formidable. 

This problem is very much simplified, however, by introducing 
a transformation of the flow properties which reduces the aligned- 
hydromagnetic-f low equations to those of gasdynamics for a fictitious - 
gas obeying a peculiar equation of state. This reduction of the 
aligned-flow equations to those of gasdynamics has been known for 
sometime (Cowley, 1960). For example, the details of this corres- 
pondence has been worked out for the case of two-dimensional flow by 
Grad (1960) for linearized motion, by Imai (1960) and luriev (1960) 
for irrotational motion, and more recently Spreiter et al. (1970a) 
have developed the transformation for axisymmetric, shock-free ro- 
tational flow. In our application we have extended this concept and 
developed the transformation for more general aligned flow which ad- 
mits shock waves in the flow region. 

For our purposes, then, a mathematical transformation of the 
aligned-fiow properties is introduced, and the resulting transformed 
variables are defined as pseudoproperties, indicated by * and related 
to the actual physical properties by the transformation equations 



h + (I - A"p/^•'iT) 

c;^'- = c; 
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From equations (5-4) and the definition for A given by equations 
(9-'l) we easily derive 


av ^ _!_ = 

w m ^2 • 

Substitution of equations (9-4) into 
and (9-3) yields 





the aligned-flow equations (9-2) 


\7 . 0 

^*(y*- V) y* + Vp* =' 0 
(V*. V) s’" = 0 


(9-6) 


V 


h^)l = V) + h^) -- 0 


and the discontinuity relations 


[/>X] = 0 

+ p*n] = 0 
4- h^] - 0 


(9-7) 
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which are identical in form to the equations of ordinary gasdynamics. 
The definition of V* and shox^s that -f- = 

=■ x^hich for uniform upstream conditions is constant 
everyx-7here and is the integrated form of the pseudoenergy equation. 

An equation of state relating the variables > p* j 

and is needed to complete the system of equations (9-6) and (9-7) . 

Although the equation of state relating p , ^ , and S is the 
same as in ordinary gasdynamics for a perfect gas, the corresponding 
relation between , and is substantially different. 

Explicitly equations (7-1) shoxv- that 


&o 

4 


yO^exp [(S - S<3o)/c^] 


(9-8) 


The derivation of the corresponding relation for the pseudovariables 
begins by using equations (9-1) to eliminate 8^/ 81T from the defi- 
nition of given in equation (9-4) and by introducing the enthalpy 

in place of V through the relation n+ = constant. There 

thus results 


P* - p + = p + pHh,.- h)/p^ 


(9-9) 


in which the stagnation enthalpy CpToo is given by 


hst = 


2 ''CD 


4- 


R 


oo 


' Pc 


y, 


CO 


OD 




2 

CO 


(9-10) 


Substituting h = ^ and equation (9-8) into equation (9-9) and 

rearranging, v?e finally"^ express p^ in terms of density and 
pseudoentropy S’’ by 
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arises frcm the substitution, into equation (9-9) of = 2(hsi' ~ h) 

= -^) ^ 0 , From equation (9-11) we can proceed to 

obtain the desired relation betoeen p*'' , 1 S’" by using 

equation (9-10) and the following relation derived from equation (9-4) 
to express p in terms of p^~ ; 


r 




Hit > 



(9-12) 


Thus the correspondence between the aligned-flow equations and 
those of gasdynamics for a hypothetical gas obeying an unusual equa- 
tion of state is complete. A significant advantage results from ap- 
plying these equations because now the partial differential equations 
(9-6) for the pseudogas are exactly equivalent to those of gasdynamics, 
and therefore we can make use of much of the numerical analysis and 
program development that has been carried out for gasdynamic applications. 
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I'Jhat these transformation equations (9-4) have accomplished, in ef- 
fect, is to absorb the Lorentz-force term of equations (9-2) into 
the inertial and pressure-gradient terms while leaving the form of 
the continuity and energy equations unchanged. However, this simpli- 
fication of the momentum, equation by the effective elimination of the 
Lorentz force comes at the expense of a more complicated equation of 
state. Nevertheless, the net result remains a significant advantage 
because v?e have reduced the number of partial derivatives which must 
be approximated by finite -difference expressions while only increasing 
the complexity of an alegebraic equation which is simply used in an 
auxiliary step to relate flow variables . 

B. Adaptation for Numerical Solution 

Once the foregoing transformation is introduced and the 
aligned-flow equations (9-2) and (9-3) are reduced to the gasdynamic 
equations (9-6) and (9-7) for a pseudogas obeying the equation of 
state p"" — , the two computer programs described in 

chapter VII and originally designed for gasdynamics can be modified 
to the flow field for. aligned .hydromagnetic flow. We have indicated 
in the foregoing section that the basic numerical approach (inverse 
iteration method for the subsonic region and method of characteris- 
tics for the supersonic region) as well as the finite-difference 
'schemes all remain intact, the only changes required' are to the aux- 
iliary relations derived from the properties of the pseudogas. Thus, 
the explicit forms of the auxiliary equations which are used in each 
of the two computer programs must noxi? be derived. 

1. Inverse iteration method 

For the inverse iteration program xvhich computes the floxj-in 
the subsonic region, the dependent pseudovariables are the two 
pseudovelocity components U* , 'h* along xfl'ith the pseudo pressure 

and density p^ and . In precisely the same xjay as outlined in 

chapter VII, the derivation of the alternate form of the pseudoenergy 
equation 
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(9-13) 


follows from the constancy — 0 of pseudoentropy along 

streamlines and the definition of the pseudosound speed 




3 P^j 3 ^ 


(9-14) 


By performing the indicated differentiation of expression (9-14) on 
equation (9-11) by means of the chain rule and equation (9-12) , we 
obtain 



or, equivalently, by a similar operation on the definition of p* 
given in equations (9-4) we obtain 



is the Alfven speed of sound- Therefore, the continuity equation and 
the two components of the momentum equation from equations (9-6) 



13-3 


together with the energy expression in the form of equation (9-13) 

supplemented by the auxiliary relation c3'^ •= P^) provide 

the four ‘partial differential equations required to determine the 

four unknowns U''^ , , p , and • In. this form we see 

that the pseudogas transformation has in effect transported the Lorentz 

force from the momentuin equation' to the energy equation (9-13) through 

the influence of the auxiliary relation for , 

Although it is not intrinsically required for the numerical 

solution of these equations, an explicit form for h" = 

can also be derived from the definition = h + f i — 

M-TT ' o 'it' ' 

or, from what is easily derived from that by use of equations ' 
(9-5), ' h* = h T { I — ) . In a manner similar to the deri- 

vation of p^' , we replace V by means of the relation h + 

= hji = constant, h by its calorically perfect form ^ » 

and p by equation (9-12), and after considerable rearranging, obtain 













1 

ar-1 



-1 

+ 


(9-17) 


It should also be pointed out now that the conservation of 
entropy S along streamlines and of stagnation or total enthalpy 
hgi every^-?here in the aligned flow carries over and is maintained 
for the pseudoentropy S and pseudo total enthalpy in the 

pseudogas flow. This is an important result because the internal 
consistency checks built into the computer program of Lomax and 
Inouye (1964) will remain valid independent checks on the accuracy of 
the computed solution of the pseudogas. Furthermore, since and 

y are parallel vectors, the streamlines for both flows will coincide, 
and similarly the definition of the stream function = 

/f(/iudx — /3TTdr) = f r dx - ) is 

equivalent for both flows. In this way the criterion of ''fl vanishing, 
which is used to determine the calculated body shape in the subsonic 
region, can still be maintained. 
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2 , Method of characteTistics . 

This close equivalence between the actual aligned flow and 
the pseudogas also extends into the supersonic region of flow where 
the pseudo Mach number M* has a value greater than unity. The 
relation between the pseudo Mach number and the two Mach numbers M 
and Ma of the aligned hydromagnetic flow is derived .from the defini- 
tion and equation (9-16) and rearranged into the form 


M 


^2 


-h “ I 


(9-18) 


From this relation we can readily see that when M = 1 also 

has a value of unity, and consequently the sonic line and transition 
region from subsonic to supersonic flow coincide for both these flows 
(see next section for further discussion about this transition). Be- 
cause of this equivalence we can use the me thod-of- character is tics 
program to compute the pseuddflow in the supersonic region. 

Since the flow equations (9-6) for the pseudogas are exactly 
the same as ordinary gasdynamics, the equation derived from them and 
applicable along characteristic lines is exactly equivalent to the 
corresponding equations (7-22) for ordinary gasdynamics although the 
shape and location of the characteristic lines will be different. To 
compute the pseudoflow by use of the present characteristic program, 
in addition to the pseudoenergy equation in integrated ' form 

+ n = = constant, and the constancy of pseudoentropy 

along streamlines, an auxiliary equation ) is 

required to complete the mathematical description. This relation, 
however, cannot be obtained in explicit form from the relation 
== ) because the equation of state P* ~ 

cannot be inverted analytically to yield P'^^ S^’) 

This finding causes no difficulty, hov?ever, because given numerical 
values for p'^ and in the computer algorithm, it is a simple 

matter to compute the corresponding value for from equations 
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(9-11) and (9-12) by a standard root-finding technique. Then with 

S^') effectively determined, the mathematical descrip- 
tion for the pseudogas corresponds exactly to that for gasdynamics 
and has precisely the form required for the method-of-characteristics 
computer program. Furthermore, since » the independent 

mass-balance check on the accuracy of this program is still maintained. 

Although the discontinuity conditions (9-7) relating pseudogas 
properties across a shock wave have exactly the same form as those of 
gasdynamics, they do not reduce to the simplified relations of or- 
dinary gasdynamics for a perfect gas because of the unusual pseudo 
equation of state. This result stems from the fact that even though 
the aligned hydromagnetic flow is a perfect gas, the bow wave in that 
flow is a fast hydromagnetic shock wave and not a gasdynamic one. 
Therefore, because both of the above numerical programs require ex- 
plicit equations relating the flow conditions across the shock wave, 
the discontinuity relations (9-7) for the pseudogas must be expressed 
in such a way that the pseudoflow variables > "^ 1*5 5 

on the downstream side of the shock wave can be determined from 
the conditions and upstream. For this purpose 

the four equations of the discontinuity conditions ( 9 - 7 ) together with 


equation (9-17) provide a determinate system to calculate the five 
do^stream unknowns , "tr,^ , , p* , and . By com- 
bining equations (9-7) and (9-17) and rearranging, we form the more 


convenient expressions 




where 
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The first equation of (9-19) is a quartic in which can be solved 

numerically if the upstream conditions and shock angle are specified. 
Once p^ is determined, then and V,^ are easily calculated 

from the remaining three equations of (9-19) . This quartic equation- 
in P^ - has four roots, the trivial solution pp ~ and three 

others which must be determined by numerically solving the remaining 
cubic equation for p^ . Over the range of values for parameters 

5 V^co i Pa , and M/ao? in our application, we found the 
solution of this cubic to consist of one real and two complex conjugate 
roots. Since equations (5-14) state that a hydromagnetic shock wave 
must satisfy the condition -jr- < we find that the nontrivial 

/^\ Cco 

real root is the physically relevant solution. With yO" determined, 
the calculation of the remaining do^rastream properties follows directly 
from equations (9-19). 

3 . Identification of flow regimes 

In the preceding paragraphs V7e have sho^jn that except for the 
equation of state the analogy between the pseudogas equations and 
those of gasdynamics is exact. Ho\7ever, if we expect to solve -these 
pseudogas equations by use of the inverse-iteration and method-of- 
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characteristic computer programs, then the character of these equa- 
tions must also correspond, that is, the pseudogas equations must be 
elliptic in the subsonic region and hyperbolic in the supersonic one. 

. The approach to gasdynamics as Bco~^ 0 and hence 

00 is clearly evident from equations (9-4), (9-5), (9-15), 
and (9-18) . But although the mathematical analog holds for all , 

it becomes less suitable when Bco increases and Mp, < 1 since /O* 
then becomes negative and the two velocity vectors V^. and V are 
of opposite sign.- As in the corresponding gasdynamic case, an ex- 
amination of the local pseudo Mach number indicates the character 

of the pseudogas equations and in turn the aligned hydromagnetic 
equations also. From an investigation of equations (9-18) we deter- 
mine the following relationship between the actual Mach number M and 
the pseudo one M ^ : 

if Mp^ > I then ^ 1 implies ~ | (9-21) 

and = I for all M when. = I . Thus for > 1 throughout 
the flow field the aligned hydromagnetic equations are elliptic in 
the subsonic flow region and hyperbolic in the supersonic region. This 
situation corresponds exactly to that for ordinary gasdynamics, arid 
therefore the application of our two computer programs to the ficti- 
tious flow is possible. A closer inspection of equations (9-18) 
reveals, hox^ever, some unusual behavior when Mp^ < 1 and consequently 
when /S' ^ becomes negative and has the opposite sign of V 

For the case when both < 1 and M < I , equation (9-18) shows 

that > i when > 1 and. <0 when -P < I . 

The first of these results indicates that even though the aligned 
hydromagnetic flow is subsonic and subAlfvenic, the flox-7 equations 
are hyperbolic when -P while the -latter indicates that the 

wave speed becomes imaginary and the equations are elliptic x<?hen 
-P < 1 1 For the case xjhen > ! , subsonic aligned 

hydromagnetic flox-7 is related to supersonic ordinary gasdynamic flox-7 
in the opposite direction. As indicated by the Friedrichs diagrams 
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in figure 6-2, hydromagnetic disturbances would then propagate along 
upstream rather than dotnistreain characteristics, and slow hydromag- 
netic shock waves would extend upstream from the body rather than 
downstream as in ordinary gasdynamics. For' ^ ^ 

is negative, has the opposite sign of and supersonic aligned 

hydromagnetic flow is related to subsonic ordinary gasdynamic flex'? in 
the opposite direction past a similar body. With this peculiar be- 
havior x\?e cannot expect our numerical programs x^hich x^ere designed 
for ordinary gasdynamics to successfully compute the pseudogas floxj. 
Nevertheless, some applications for < I flox^s past obstacles 
have been carried out by Seebass (1961), Tamada (1964), Geffen (1965, 
1966) and others . 

This behavior of M'" is illustrated in figure 9-1 x^hich 
shox7S the range of values for for various pairs of values for 

M and . Also included in that figure is an isentrope (dashed 

line) indicating variation of as a function of M along a 

streamline. The explicit form of that function is given by the ex- 
pression 




2 hsi 


Vooexp 


S~Sco 

Cv 


e(ir-o 


(9-22) 


derived from the definition of Mach number M , equations (9-5) and 
(9-8), and the relation + h = • From this equation it is 

clearly evident that along a given streamline is a monotonically 

increasing function of M . Since entropy reaches a maximum value 
and M a minimum value at the stagnation point, equation (9-22) 
indicates that also reaches i-ts minimum value there. 

The value of , x-7hich depends on Mco and Mp^co > therefore, 

determines x\?hether the condition x-zill occur in any region 

of the flox7. For example, as illustrated in figure 9-1 if — 1 

then both and increase monotonically from the stagnation 

point to the free stream along the stagnation streamline just as for 
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Figure 9-1. Variation of M for various ranges of values 
• fox M and • The isentropes display the 

behavior of as a function of M according 

to equation (9-21) along the stagnation stream- 
line , 
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ordinary gasdynamics, and never becomes less than unity at any 

point in the flow field. On the other hand if < 1 , Ma 

■Still increases monotonically , but M" varies in a complicated 
manner: it is imaginary at the stagnation point, then decreases to 

- OO j jumps up to + CO , and then decreases doT-m to a minimum which 
is less than 1 , and finally increases monotonically up to the free- 
stream value. The condition, then, that M* increases monotonically 
and that the pseudogas equations have the same character as those of 
ordinary gasdynamics is easily determined from equations (9-5) to be 


which for hypersonic flow would usually be satisfied if > 2 . 

If condition (9-23) is satisfied, then both the supersonic and sub- 
sonic regions of the pseudogas flow will correspond directly with 
those of ordinary gasdynamic flow, and therefore the inverse- itera- 
tion-method and method-of-charactexistics computer programs will be 
suitable for computation of the pseudoflow. 

Actually, difficulties arise even before the free-stream 
Alfven Mach number reaches its lower limit. An examination 

of equation (9-15) for and (9-17) for reveals that when 

input values for Maoo become lower than about VTo , the denom- 
inators of both these equations vanish at some point in the flow field, 
and both and n would be infinite at that point unless the 

numerators also vanish simultaneously, , If the critical values for 
the flow properties at this point are designated by the subscript 
cr , we see from the denominators of those equations that 
Pec ~ {? ” I ) ^ or equivalently after use of equation (9-12) 

that Pz.'c~ 05-0/5 ^ similar analysis of the .numerators 

yields Ptr = , and they also vanish. Both 

quantities c3^’ and , then, are of indeterminate form at this 

point , This situation actually stems from the fact that if the 
pseudoentropy is considered as a surface described by equations 
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( 9 - 11 ) and ( 9 - 12 ) over the — plane, then S'^r is multivalued 
at the critical point . However, an analysis of the limit 

process for these equations shows that as the point p^^ , is 

approached in the plane, a continuous and finite limit for 

which depends on the direction of approach to that point is also 
reached, and in addition, that the indeterminate form of equations 
( 9 - 15 ) and ( 9 - 17 ) actually leads to a continuous finite value for 
and ’ at the critical point. Some limit values for S^r 

corresponding to certain approach directions in the p^— plane, 
however, are not physically realistic and incite the growth of in- 
stabilities in the numerical computation of the flow field. To con- 
trol the instability we have developed an algorithm that restricts 
the path of approach to the critical point to only those that lead to 
a physically realistic limit for were able to get suffi- 
ciently smooth values for , and in the vicinity of 

the critical point so that the computation scheme was stable. 

Even with these added procedures the lowest value of M/\co 
for which we could compute an accurate solution was for = 2.5 

Uncontrollable numerical instabilities associated with the 1" co- 
ordinate line becoming duly inclined and consequent deterioration of 
accuracy as indicated by the consistency checks on the conservation of 
S^' and as well as increasing discrepancy between the calculated 

and desired body shape prevented solutions being computed for lower 

. In fact in retrospect it appears rather remarkable .that even 
for values of 2.5 and higher, the shock-shape equation 

( 7 - 21 ) which produced ellipsoidal body shapes to a high degree of 
accuracy for gasdynamic flow of a perfect gas also produces similar 
body shapes to the same accuracy for aligned hydromagnetic flow. 
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ClIAPTER X 

FEATURES OF THE EXACT SOLUTION 

We have shown in the last chapter that with the introduction 
of the pseudogas transformation the equations for aligned hydromag- 
netic flow can be reduced to those of gasdynamics supplemented by an 
unusual equation of state. Once the required auxiliary relations 
derived from the pseudo equation of state have been determined and 
written into the inverse iteration and method of characteristics 
computer programs, all that remains is the specification of the input 
parameters Mco and then the solution to the pseudogas 

equations (9-6) and (9-7) can be computed. It is then a simple matter 
to invert the transformation equations (9-4) and convert the pseu- 
doflow properties back to the actual aligned-flow properties which 
satisfy equations (9-2) and (9-3). As discussed before this solu- 
tion now depends on the two input parameters and j iia- 

stead of the single parameter Mqj as was the case in the gasdynamic 
theory. Thus, once we specify 1 j we can 

compute the hydromagnetic solution to equations (9-2) and (9-3), and 
therefore determine what effect variation of the parameters 
and Mco > representing different conditions in the oncoming flow, 
has on the aligned-flow properties. And even though this investiga- 
tion is based on the solution to the aligned hydromagnetic flov? 
equations, study of the Friedrichs diagrams suggests that qualita- 
tively similar effects may be anticipated for other alignments as 
well. 

A . Broad Features of the Aligned -Flow Solution 

Together with the input Mach numbers, we must also specify 
an axisymmetric planetary boundary. For all the calculations presented 
in this chapter, we have chosen to compute these hydromagnetic solu- 
tions for aligned-flow interaction with the magnetosphere boundary of 
the Earth partly because of the interest in and enormous amount of data 
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measured in near -Earth space, but mainly because application of the 
correspondence rule quickly converts any of the properties of the 
flow around the Earth's magnetosphere into those for flow aroun’d the 
ionosphere of either Mars or Venus by a simple relabeling of the co- 
ordinate axes. Using the magnetopause shape illustrated in figure 6-6 
for the magnetosphere boundary, we have computed aligned-flow solu- 
tions for various sets of input Mach number Moo and Alf-ven Mach 
number in the range 5 < iMqo < 20 and 2.5 ^ ^ 20 . A 

comparison of the different solutions obtained for differing values 
of input Alfven Mach number but the same value of Mach number 

immediately reveals qualitative variations in these solutions. 

The most apparent distinction noticed among these solutions for dif- 
fering values of is the large diversity in the shapes and loca- 

tions of the fast hydromagnetic shock v?ave that forms upstream of the 
planet. To illustrate this effect we have plotted in figure 10-1 
the calculated bow-wave positions derived in the computed solutions 
for aligned hydromagnetic flov? using values for of 2.5 , 

3 , 5 , 10 , and 20 but while maintaining = 5/3 and the Mach 

number Mco constant at a value of 10 , v/hich is representative of 
normal conditions in the solar wind. For ease of comparison the cor- 
responding bow-wave position determined by using Mco ~ bhe 

approximate gasdynamic theory described in chapter VII is also in- 
cluded in figure 10-1. Easily noticed in this figure is the trend 
that as Maco decreases the position of the box? wave for aligned flox? 
further departs from that for the gasdynamic boxv x^ave. The departure, 
is such that as ^p^co diminishes the hydromagnetic shock-x-jave. loca- 
tion becomes closer to the magnetosphere boundary in the nose' region 
than does the gasdynamic one but moves substantially further axvay 
from the magnetosphere tail dox-rnstream from the Earth. The nature of 
this variation in box-j-x^ave position as Mao) decreases is almost as 
if the bow x-?ave xjcre rotated, counterclockx-rise about the point of in- 
tersection of all the box-j x^aves illustrated in figure 10-1. No partic- 
ular significance, hoxjever, is attributed to this point. For Macs 
= 20 the hydromagnetic and gasdynamic shock xvaves nearly coincide, 
and the approach to gasdynamics as f'^^Aco further increases is clearly 













145 


evident . 

This migration of the hydromagnetic bow v?ave from the posi- 
tion of the gasdynamic bow wave as decreases can be understood 

by comparing the relative values for the pseudo sound speed In 

the hydromagnetic solution and the speed of sound c5 in the gasdy- 
namic one. Such a comparison was’ carried out, and it revealed that 
throughout the flow <3* was greater than <3 except in the nose re- 
gion of the magnetosphere where c3 was less than cl . These rela- 
tive differences become greater as diminished. Since the 

pseudovelocity and the gasdynamic velocity V are also about 

the same magnitude, the higher value of c3'^ means that infinitesimal 
disturbances in all regions of the aligned flow except near the nose 
of the magnetosphere travel faster, relative to the flow speed , 

than do similar disturbances in gasdynamic flow and, consequently, 
propagate further upstream in a given time interval before coalescing 
into the bow wave. Conversely, in the nose region of the magneto- 
sphere the lov;er value of c3'^ indicates that these disturbances 
travel slower relative to than do corresponding ones in gasdy- 

namics and, therefore, do not propagate as far before coalescing into 
the bow wave. As diminishes, this propagation effect .increases, 

and in this way the aligned hydromagnetic shock wave reaches the posi- 
tions illustrated in figure 10-1 for the various values of 

In relation to the magnetopause the position of the hydro- 
magnetic bow v?ave in the nose region is best characterized by the 
standoff distance defined as the distance along the axis of symmetry 
between the magnetopause and bow wave. Not only does figure 10-1 
distinctly show that this distance diminishes with decreasing M/^co 
but also the numerical calculations indicate that the extent of the 
subsonic-flox'j region likewise contracts with a diminution of M^co 
We should point out that a decrease in just the opposite 

effect as does a corresponding decrease in the Mach number iu 

ordinary gasdynamic flow. However, this behavior connected with 
diminishing actually aided the stability of the computation 

scheme because a decrease in the standoff distance means fewer steps 
in the marching technique and a shrinking subsonic region aided 
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control of the associated numerical difficulties in the neighboring 
transonic region. 

All o£ the hydromagnetic shock-wave shapes presented in 
figure 10-1 have been computed for varying values of input parameter 
while holding the other input parameter Mco fixed at a con- 
stant value. The combination of the definition of Mach number and the 
calorically perfect relation for internal energy 6 leads to the 
equation 

^V'"= vg. y fe 

m 

which indicates that the square of the Mach number expresses the ratio 
of the directed energy to the internal energy. Thus, the free-stream 
Mach number specifies the relative gasdynamic energy in the fluid of 
the oncoming flow but not its magnetic energy. Heretofore, in gas- 
dynamic theory the shock-xjave position for a given magnetopause shape 
was solely determined by the free-stream Mach number Mco > which 
in turn is set by the values of Vco and Tco . Figure 10-1 shows 
that for- > 10 this is still true for hydromagnetic shock waves. 

However, since for < 10 this figure indicates that the observed 

variation in bow -wave positions is also a function of Maco > the 

I ** 

)/(B^/ST )J for input parameter 

M/\(x> suggests that it is the ratio of the kinetic to the magnetic 

energy of the incoming flow as well as its relative gasdynamic energy 

/ 

that determines the hydromagnetic shock-x7ave position for low values 
of M^(X) . For a given free-stream velocity Vca this ratio, and 

hence MyCi^jj) , can only be made small by either unusually .small 
values of density or unusually high values of magnetic field 

Bco in the undisturbed flow. If we hold the free-stream velocity 
Vco and temperature Tco constant, and hence specifying Moo j 
we conclude on the basis of the results presented in figure 10-1 
that either a low value of or a high value of Bco or a 
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combination of those two effects can cause the position of the hydro - 
magnetic shock wave to substantially differ from that of the cor- 
responding gasdynamic shock wave calculated for the same value of 
Moo • 

Although most of the time the flow properties of 'the solar 
wind are steady and are approximately the values indicated in figure 
2-3, several instances of unusually low density or high field strength 
have been observed. For example, for nearly the entire day of July 31, 
1967 Lyon ^ (1968) reported that Explorer 35 in the vicinity of 

the Moon observed an anomalously low solar -wind flux of about 0.2 
protons/cm flowing at 200 km/sec and that these results are cor- 
roborated by data from Mariner 5 in flight to Venus at 0.84 a.u. 

Both magnetometers on Explorer 35 (Colburn ^ al. , 1967; Ness et al. , 
1968) sho^? that the intensity of the interplanetary magnetic field 
was about 7 ^ during this time. Substitution of these values into 
equations (5-4) leads to a value for of about 0.6 . Unfortun- 

ately, no data on the Earth's box-? x^ave x^as reported on this date so 
that the effect on the position of the boxv' x^ave remains unknoxm. 

Hox7cver, Ungstrup (1971) has recently analyzed data that show 
a very strong correlation betx>7een increasing magnetic -fie Id strength 
in the solar x^ind and shock-wave movement out from the magnetosphere 
tail. He has examined magnetometer and plasma-detector data measured 
simultaneously by the satellites OGO 5, Pioneer 9, and Explorer 33 and 
35 covering a x^ide region of near -Earth space during a 36 -hour period 
spanning November 17-18, 1968. The magnetometer data from Explorer 33 
revealed the remarkable occurrence of the satellite folloxjing the 
movement of a segment of the boxj x^ave doxvnstream from the Earth for 
a distance of more than 18 Earth radii during the 36-hour period. 

From these data and simultaneous observations made by the three other 
satellites, Ungstrup (1971) x?as able to determine that not only did 
the density of the solar x^ind remain at a practically constant value 

3 

of about 2 protons/cm but the velocity and temperature, and hence 
the Mach number , also remained relatively constant throughout the 
entire 36-hour period. The da.ta indicated that the only varying flox'7 
property of the solar viind x-?as the magnecic field X’/hich linearly 
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decreased from 13 to 3.5 J during that 36-hour time interval. 

Since the magnetometer data from Explorer 33 indicated that the rear 
portion of the bow wave downstream of the Earth traveled steadily 
inward from its unusually large distance from the magnetosphere tail 
towards its more usual position as predicted by gasdynamic theory 
during that same time interval, Ungstrup (1971) concluded that the 
large magnetic-field strength in the solar wind was responsible for 
the unusual position of the bow wave and that the decreasing field 
strength brought about the movement of the bow wave back towards its 
usual position. Since a linear decrease in magnetic-field strength 
from 13 to 3.5 ^ corresponds, for the steady plasma conditions 
'of the solar wind reported by Ungstrup (1971), to a linear increase 
in Alfven Mach number . from about 2 to around 8 , the space 

observations which he analyzed as well as his conclusions are in 
complete agreement with the trend displayed in the theoretical curves 
presented in figure 10-1 for varying positions of a hydromagnetic ^ 
shock wave for differing values of * Such space observations 

confirm the validity of our hydromagnetic fluid model and the con- 
clusions which are dra^wi from it . 

With the applicability of this refined model substantiated, 
we can now proceed to examine the effects of input Alfven Mach number 
bn the aligned -flow properties which are calculated as part of the 
numerical solution. For each individual solution using a given set 
of input parameters Mco and main flow properties of 

interest are the density ratio p! poD » velocity ratio V/Vco > 
and magnetic-field-strength ratio B/Bco well as the streamlines 
of the flov7 which indicate the direction of both the magnetic-field 
and velocity vectors throughout the flow. The temperature ratio is 
also of interest, and since the energy equation for aligned hydro- 
magnetic flow integrated to exactly the same form as in gasdynamics, 
equation (8-1), \?hich is derived from it and relates velocity to 
temperature, still holds. Thus the contours of constant -temperature 
ratio T/Tco coincide with those for velocity ratio V/Vco » ^ud 
specific values for T/T^o can be -determined by a simple application 
of equation (8-1) . We present in figures 10-2 through 10-5 these flow 
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properties of four individual solutions for aligned liydromagnetic 
flow computed with input conditions of Me© = 10 and ^ =5/3 for 
each of the four values for M/\co 2.5 , 5 , 10 , and 20 . 

For each of the four solutions presented in figure 10-2 to 10-5 four 
-subfigures have been dra™ to illustrate lines of constant -property 
contours for (a) density ratio P/P^ s (b) velocity V/Vcd and 
temperature T/T^o ratios, and (c) magnetic-field- strength ratio 
B/Bco 5 well as (d) direction of magnetic -fie Id and velocity 
vectors. 

A comparison betvjeen constant density and velocity contours 
from the hydromagnetic solution for Maco “ ^0 presented in figures 
10-5 (a) and (b) and the corresponding contours from the gasdynamic 
solution in figure 8-2 shows that the differences between these two 
solutions are practically indistinguishable, and clearly indicates 
that the hydromagnetic solution approaches that of gasdynamics as 

further increases. A similar comparison between these gasdy- 
namic properties and the corresponding hydromagnetic properties in 
figures 10-4 (a) and (b) computed for shows that the 

differences between these two solutions are also relatively small. 
Hox-yever, an examination of all four solutions presented in figures 
10-2 to 10-5 reveals some general qualitative trends in these flox-? 
properties as the free-stream Mach number decreases to the value 2.5 
For example, figures 10-2 (a) to 10-5 (a) shoxj that as dimin- 

ishes the zone of high-density ratio in the nose region of the mag- 
netosphere shrinks xihile in the remaining part of the flox^ the density 
ratio approaches unity. Hence the overall effect of decreasing 
is a decrease in both the compression and expansion of the flowing 
plasma. Mearaihile, figures 10-2 (b) to 10-5 (b) indicate that a 
decrease in produces a slight increase in the velocity ratio, 

and consequently a small decrease in the temperature ratio, uniformly 
throughout the entire flox<r field. On the other hand, figures 10-2 (c) 
to 10-5 (c) manifest that diminishing causes the magnetic- 

field-strength ratio to increase slightly in the flox-? region near the 
nose of the magnetosphere but to approach unity in the remaining 
part of the flow. Finally, figures 10-2 (d) to 10-5 (d) show that as 
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Figure 10-3. Various field properties for aligned hydromagnetic 
flow past the magnetosphere for MAoo = 5 , 

Mco = 10 , and X = 5/3 represented by lines of 
constant-property contours for; 
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Figure 10-3 (d) , Direction of magnetic -field and velocity vectors. 
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Figure 10-5. 


Various field properties for aligned hydromagnetic flow past 
the magnetosphere for == 20 , Moo = 10 , and 2T = 5/3 

represented by lines of constant property contours for: 

(a). Density ratio ' 
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^■^Aco 'iecreases the angles that the streamlines make with the axis 
“of S 3 mimetry also decrease somewhat throughout the entire flow except 
for that portion of the nose region where they remain relatively con- 
stant. Since a decrease in implies an increase in the mag- 

netic-energy content of the oncoming flow relative to its kinetic 
energy, v?e thus conclude from the foregoing results that as 
decreases it becomes increasingly more difficult for the flow to 
stretch and distort magnetic -field lines. Consequently, the density 
and field strength are less disturbed and tend to remain closer to 
their free- stream values while the velocity tends to increase some- 
what because of the smaller flow deflection. Although these four 
solutions presented in figures 10-2 through 10-5 are all for the same 
value of Mach number , we have also computed similar solutions for 
this same set of values for but for different values of 

and found the same variation in flow properties for decreasing 
as we did for the case with Mg 3 ~ 10 . These effects, then, appear 
to be controlled by the value for f'^Aco ’ hence and Bco 

and to be relatively independent of M ui . 

B. Comparison Betoeen Approximate and Exact Solutions 

In the foregoing section we have discussed in a rather 
qualitative way some of the broad features of the aligned hydromag- 
netic solution and its behavior for decreasing ^Aco ' • this 
point it is appropriate to carry out a more thorough comparison of 
these solutions for varying M^co ultimately a collation of them 

with the more approximate gasdynamic results. 

A quantitative comparison of the shock-wave position is most 
easily and adequately performed by plotting characteristic distances 
of the shock position from the 'planetary boundary. The two measure- 
ments we chose are the shock standoff distance Xg to characterize 
the nose region and the radial coordinate measured from the axis 

of symmetry to the shock wave at a distance 2.5 Pq doimstream of 
the Earth, Since, for approaching infinitely large values, 

the hydromagnetic solution reduces to that of gasdynamics, we have 
plotted in figure 10-6 two groups of curves depicting the two 
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distances and versus order that values on the 

ordinate represent the distances for a gasdynamic shock wave . In 

this figure three curves in each group are drawn representing three 

sets of solutions over the given range of values 5 , 10 

and 20 for . For decreasing values of both these 

curves clearly indicate a decrease in Xs ^.s the bow wave moves 

closer towards the nose of the magnetosphere and an increase in the 

shock coordinate T as the downstream segment of. the shock moves 

farther from the magnetosphere tail. Furthermore, the similarity of 

these curves for different values of Mco shows that the relative 

shock movement for decreasing is independent of the value of 

Mco although the absolute shock position does vary v?ith Moo 

Finally, we should point out that these cur'ves indicate the smooth 
/ 

approach to gasdynamics for increasing ^Aco » that as the value 

of decreases below 10 the position of the hydromagnetic shock 

departs significantly from that of the gasdynamic shock. 

To compare this effect with the variance of position of the 
gasdynamic bow shock for differing values of M co » in figure 10-7 
we have plotted versus I/Mco corresponding values for Xg 

and for the shock calculated in the gasdynamic solution. This 

figure indicates that the distance for the gasdynamic shock 

varies with M^o in roughly the same manner as the corresponding 
distance Pg for 'the hydromagnetic shock does with f^Aoo figure 
10-6. This equivalence, however, does not carry over to the shock 
standoff distance Xg . For this distance figure 10-7 shox^s that 
Xg- varies with Mco in just the opposite sense that the corres- 
ponding for the hydromagnetic shock does X7ith ^A'® ' i ®3 

decreasing Mco causes the gasdynamic shock standoff distance to 
increase while decreasing brings about a decrease in the 

hydromagnetic X^ . This result is just another manifestation of 
the differing natures of the txro input parameters Mco nnd Maoo • 

Figure 10-7 further provides a suitable examination of the 
effect of using, as advocated by Dryer (1970) in an attempt to in- 
clude some effect of the second input parameter value 

of the free -stream pseudo Mach number M'co for the input Mach 
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number in the gasdynamic theory. For example, if the free-stream 
parameters Mco and aligned flow both have a value of 

10 , then equation (9-18) determines Mco to be about 7.09 • 

With dashed lines we have indicated this value on the abscissa for 
the curves in figure 10-7 so that the values for X5 and re- 
sulting from the gasdynamic solution using instead of 

as the input parameter can be easily read off. The dashed lines in 
this figure show that for this case the standoff distance would have 
a value about = 0,29 Pq while the shock ordinate would be about 
fg = 3.7 To . These values together with those obtained by using the 
gasdynamic theory in the usual way with Mco instead of Mco can 
nov7 be compared to the exact values determined by the hydromagnetic 
solution. For Pg figure 10-7 shows that gasdynamic theory using 
Mco = 10 calculates a value fg = 3,57 To as compared to 3.7 To 
obtained from the same theory using Mco =7.09 . Comparing these 
tv70 values with the exact value =3.66 To displayed in figure 
10-6 for Mcd " 10 and = 10 , we see that use of in the 

gasdynamic theory does improve the accuracy of the shock ordinate 
calculated downstream from the Earth by that theory. However, 
in a similar comparison of the value for Xg figure 10-7 indicates 
the value to be Xg = 0.26 Pq for the gasd37namic theory using 
Mcd = 10 and = 0.29 P© using M® = 7.09 while figure 10-6 shows 
the exact value to be Xg = 0.25 Pq . Therefore, the use of 
in place of Mco actually lessens the accuracy of the gasd3Tiamic 
theory in calculating the standoff distance. Because such use of 
Mlio in the gasdynamic theory increases the accuracy of that theory 
only for that part of the flov7 downstream from the Earth while re- 
ducing it for the flow upstream, Mco appears to be the best overall 
choice of input parameter to use in the gasd37namic theory for approx- 
imating the hydromagnetic 'solution. 

C . Synopsis of the Foregoing Results and Conclusions 

The comparison between the solution of the gasdynamic equa- 
tions and the solution to the more exact aligned hydromagnetic equa- 
tions representing our fluid model of the solar-v7ind interaction has 
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been presented in this chapter , and it has verified the accuracy and 
applicability of the simpler gasdynamic theory for modelling the solar- 
v?ind interaction for conditions in the solar V7ihd when a 

value around 10 or higher . When this provision is fulfilled, all 
of the previous theoretical results' based on gasdynamic theory that 
has been derived in an analysis of the solar-wind interaction with 
the Earth, Mars, and Venus remain accurate and valid. However, when 
the value of lower than 10. , the aligned-flow solu- 

tion indicates that substantial discrepancies exist between the 
exact hydromagnetic theory and the simpler gasdynamic one. These 
differences manifest themselves in not only the values of the flow 
properties but even the gross features of the theory such as the 
shape and location of the bow shock wave. This result shows that the 
position of the bow wave is not only determined by the values of 

and Tco undisturbed flow as is the case in gasdynamic 

theory, but also by the values of and Bco which fix the value 

of indeed the validity of this conclusion based on the 

aligned-flow solution has been supported by the satellite measure- 
ments’ made in space and reported by Ungstrup (1971). Eor the condi- 
tions in the solar wind at that time, he concluded that the unusual 
position of the bow wave was determined by the strong magnetic field 
and that it varied with decreasing field strength in precisely the 
manner indicated by our -theoretical hydromagnetic analysis . 

Finally, we should point out that although these results and 
conclusions are all based on the solution of the hydromagnetic equa- 
tions (9-2) and (9-3) for aligned flow, the study of the Friedrichs, 
or phase and group velocity, diagrams in figure 5-1 as vjell as the 
good agreement with space observations shows that perfect alignment 
of the Vco and Boo vectors is not essential, and qualitatively 
similar effects can be anticipated for other alignments as well. 
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CHAPTER XI 

COMPARISON WITH ALIGNED -FLOW OBSERVATIONS MADE IN SPACE 

In the last two chapters we have developed the nondissipative 
hydromagnetic theory. for the case of parallel velocity and magnetic 
fields and presented a numerical procedure for solving the resulting 
align'ed-flow equations and determining the position of the bow shock 
wave and the hydromagnetic flow properties of the shocked plasma. In 
the previous chapter we have qualitatively examined these computed re- 
sults for a variety of free -stream conditions, compared them with cor- 
responding results from gasdynamic theory, and have recognized trends 
in the general character of the hydromagnetic solution for decreasing 
free-stream Alfven Mach number . Since this preliminary com- 

parison indicated very encouraging agreement with previous theoretical 
and experimental investigations for large with measure- 
ments observed in space for low now appropriate to 

carry out a more detailed, quantitative comparison with space obser- 
vations in order to further verify the accuracy and validity of this 
aligned hydromagnetic model. 

For such a comparison space observations made at a time when 
the velocity and magnetic -fie Id directions in the solar wind were 
reasonably aligned is obviously most suitable. One instance of such 
experimental data is the measurements of the plasma and magnetic - 
field properties of the solar wind made' by Mariner 5 as it flev 7 past 
Venus. The magnetometer and plasma-probe experimenters (Bridge ^ al. , 
1967) found that the interplanetary plasma conditions in the solar 
wind x<rere quite steady, and although the direction of the magnetic 
field showed some fluctuations its magnitude was nearly constant.. 

They do not attribute these fluctuations in field direction to changes 
in overall interplanetary conditions, but rather to a series of local 
solar-wind structures, perhaps filaments, which were convected past 
the spacecraft. For our purposes this localized character can be dis- 
regarded. The nature of the average magnetic field encountered' along 
the spacecraft’s trajectory is illustrated in figure 11-1 as originally 
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Figure 11-1. Mariner-5 trajectory and magnetic-field vectors. 

The three panels contain aphrodiocentr ic- solar - 
ecliptic projections of the Mariner trajectory and 
of the average, measured field at specific points. 
(Bridge ^ , 1967). 
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presented by Bridge et (1967). The three parts of this figure 
show the usual orthographic projections of both the Mariner trajectory 
and the measured fields. The aphrodiocentric-solar-ecli'ptic (ASE) 
coordinates are defined by a plane parallel to the ecliptic passing 
through the center of Venus and by j the projection in this 

plane of the direction to the Sun. The axis points toward 

the north ecliptic pole, and completes the orthogonal, right- 

handed system. The circled numbers through correspond to 

the labels added by the experimenters to the plasma and magnetic- 
field data in figure 8-13 and denote features of special significance. 
The field vectors are not shown at equally spaced intervals but at 
times that best illustrate the character of the average magnetic 
field within the various regions along the spacecraft's trajectory. 

That the magnetic -field direction is reasonably parallel to 
the velocity direction which is radially outward from the Sun can be 
determined by study of the orthographic projections in figure 11-1. 
The projection, which is viewed from the cross-section 

edge of the plane parallel to the ecliptic passing through the center 
of Venus and along the line between the Sun and Venus, shows that the 
magnetic field remains nearly parallel to the Sun-Venus line through- 
out the time interval displayed. Furthermore, the small components 
of the magnetic field illustrated in the T^gg " ^ASE P^^j^^^tion, 
which is viex^ed from the Sun toward the planet, also indicate that 
the field direction is primarily parallel to the ecliptic and in a 
radial direction from the Sun during the interval. And finally, the 


X, 
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projection, viewed from the north ecliptic pole, con- 

ASE 


firms that the magnetic field lies in the ecliptic plane, and a com- 
parison with the theoretically calculated field direction displayed 
in figure 10-3(d) for aligned flow doxmstream of a shock wave shows 
that the observed direction of the field is indeed consistent with 
that for aligned flo\7. We "then conclude that at the time these data 
were measured the average magnetic -field direction was reasonably 
aligned with the velocity direction, and thus these data are very 
appropriate for comparison with the theoretical calculations made with 
our aligned-flow model. 
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Although there have been other examples of measurements made 
in space which show the solar wind to be more steady and the velocity 
and magnetic-field directions to be more perfectly aligned, we choose 
to make a detailed comparison vjith these near -Venus data from Mariner 
5 for two particular reasons. The first is that in chapter VIII we 
have already carried out a comparison between these data and the 
theoretical positions of gasdynamic shocks arising from interaction 
with several of the ionopause shapes illustrated in figure 6-6 and 
found general agreement. An additional comparison of these same data 
with the calculated results of a more refined interaction model would 
therefore enable us to also contrast the theoretical results cal- 
culated using the simpler gasdynamic theory with those from the more 
refined hydromagnetic theory with reference to these Mariner~5 ob- 
servations. But the primary reason is, however, that the values of 

3 

about 590 kra/sec , 4 protons /cm , and 8 ^ reported by Bridge 

et al . (1967) for the velocity, ion number density, and magnetic- 
field strength in the solar wind before and after Mariner 5's en- 
counter with Venus indicate a value of about 6.75 for ^/:\co ’ 

our analysis of the aligned hydromagnetic theory presented in chapter 
IX showed that when substantial differences between the 

theoretically calculated hydromagnetic and gasdynamic results begin 
to appear . Thus on this theoretical argument alone we should expect 
better agreement to exist between the hydromagnetic results and the 
observations of Mariner 5 than existed for the corresponding gasdy- 
namic results . In fact at this point V7e can already predict better 
agreement with the hydromagnetic theory on the basis of the following 
argument. Bigure 8-13 shows that the gasdynamic shock wave calculated 
for Mco ~ 5 and associated with'the ionopause shape for H/r^ = 0.25 
indicates that the theoretical shock is positioned somewhat closer to 
Venus than was observed. However, since we have learned in chapter X 
that for 10 the position of a hydromagnetic shock wave is, 

in all regions except near the nose, somewhat farther away from the 
planet than the corresponding position for a gasdynamic shock, hydro- 
magnetic theory V7ill predict a shock position that better agrees with 
the observations than did the gasdynamic theory. This good agreement 
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between the calculated position of the shock wave associated with the 
ionopause shape for H/t^ == 0.25 and that observed by Mariner 5 
suggests that the upper ionosphere consists primarily of ionized 
atomic hydrogen at a temperature of about 700 °K corresponding to a 
scale height of around 1500 Ion . On the other hand, however, the 
theoretical results displayed in figure 8-13 for H/i^ = 0.25 indicate 
that Mariner 5 should have then crossed the ionopause, but there is 
no certain evidence that Mariner 5 actually did penetrate the iono- 
pause. As discussed in chapter VIII, a possible explanation for this 
discrepancy is that the tail of the ionosphere might taper inward 
toward the axis of symmetry rather than extending straight downstream. 
The theoretical ionopause shapes displayed in figure 6-6 give no 
indication of such a trend, but they are of low reliability for this 
feature of the flow because of the deterioration of the quality of 
the Nex^tbnian approximation for both the gas and magnetic pressure of 
the solar wind as approaches 90 degrees. This aspect of our model 
can be refined, however, by using the exact pressure of the solar wind 
calculated in the numerical solution instead of the approximate New- 
tonian expression and thus determining the exact shape of the iono- 
pause as part of the solution. Such a calculation would then accurately 
represent the shape of the ionopause tail for the model nresented. 


A . Calculation of the Exact Ionopause Location 

The ionopause shape is defined as the locus of those points- 
at which the tangential discontinuity conditions (6-12) are exactly 
satisfied. For all the boundary shapes illustrated in figure 6-6, the 
calculated results presented in figure 8-12 imply that although the 
discontinuity relations (6-12) are not satisfied exactly on any of 
these shapes the discrepancy is larger on the flanks than on the nose 
section of each boundary. In fact we found that the necessary incre- 
ment in the radial coordinate r of the present ionopause shapes 
needed to satisfy the discontinuity relation 



= P.«T= (Rr PJ P( - + Pv/ 
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on the f onward part of the boundary is imperceptibly small on the 
scale of figure 8-12 and would not change the calculated results for 
the exterior flow in this region. On the other hand a similar pro- 
cedure for the coordinates of the ionopause flanks shows that large 
increments are necessary, thus indicating serious errors in the present 
shape of the ionosphere - tail . This error can be corrected by calculat- 
ing a new shape for the ionopause. Since the error lies in the region ■ 
where the flow is supersonic and the equations hyperbolic, the posi- 
tion of the boundary is not required to be specified a. priori , . but 
rather can be calculated by the method of characteristics as part of 
the solution of the exterior flow. A simple numerical procedure to 
calculate the position of these nev? boundary points as well as the 
exterior-flow properties such that equation (11-1) is satisfied at 
those points was developed. A representative boundary point G is 
indicated in figure 11-2, and the unluiown flow properties at point C 
are to be determined from the know flow properties at points B and 
D by using the aligned-flow equations (7-22) to (7-25) in the follow- 
ing V7ay. We use a crude predictor and assume that the ionosphere 
boundary and right-going characteristic lines are straight and then 
solve for the location f of the point C . With this value of F 
the ionosphere pressure calculated by equation (11-1) 

and then {p *i" therefore be determined at 

point C by use of equation (6-12) . Thus v?ith p* and now 

known, we then solve the right -going characteristic equation (7-22b) 
for the flow angle 0 ^ at point C . Finally, using this value of 
0Q , we can compute the average of the flow angles at points G and 
D and the average of the angles at points G and B , and 

then correct the position of point C using these averaged angles to 
specify straight lines from points B and D . This corrector process 
can be repeated till the new position of point C differs from the 
old by less than 0.01 percent. (For all of our solutions this pro- 
cedure alxvfays converged in less than 20 iterations.) In this way, 
then, the precise shape of the ionopause, on V7hich the tangential 
discontinuity conditions (6-12) are satisfied exactly, is determined 
simultaneously as part of the solution for the exterior flow. We 
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can thus solve the exact free-houndary problem and no longer need to 
decouple the calculation of the ionopause shape from the exterior-flow 
solution as was done before. 



Figure 11-2. Sketch of characteristic lines for determination 
of the location of, and flow properties at, a 
representative mesh point C on the ionosphere 
boundary . 

B • Comparison with Mariner -5 Data 

Since the pressure of the external .hydromagnetic flow is 
balanced by the ionosphere pressure at the interface, the exact shape 
of the ionopause depends, on the input parameters of the exterior flow 
and the expression (6-11) for the ionosphere pressure. In order to 
make a comparison of the calculated aligned -flow properties around 
this new ionopause shape with those observed by Mariner 5, values for 
these input parameters which correspond to conditions in the solar 
V7ind at the time of observation must be determined. For the exterior 
flow we have already found = 6.75 while values of about 590 

Icm/sec and 300,000 "K reported for the velocity and temperature 
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in the solar wind before and after Mariner 5’s encounter with Venus 
indicate a value of about 6.47 for Mco • For the interior pres- 
sure we choose a value for H/Pq of 0.2 corresponding to a Venusian 
ionosphere composed primarily of atomic hydrogen at about 600 " K 
because it agrees reasonably well with both the dual-frequency-occul- 
tation and Lyman-alpha measurements made by Mariner 5 on the dayside 
as well as nightside ionosphere. Similarly, the pressure in the 
ionosphere wake Bl/S'JT is about 5 x 10 dyne/cm 

which, in view of the uncertainties attendant with these measurements 

made in space, approximately agrees with the value of 800 electrons/ 

3 

cm at a temperature of about 675 °K for the electron density and 
temperature in the ionosphere tail reported by Fjeldbo and Eshleman 
(1969). 

Figures 11-3 (a) to (d) illustrate the aligned hydromagnetic 
solution for the exterior flow represented by lines of constant- 
property contours for (a) density ratio pf Pa> » (h) velocity ratio 
V/Vto > (^) magnetic-field-strength ratio B/Boo j 

rection of the velocity and magnetic-field vectors as well as the bow- 
wave position and the exact ionopause shape all of which were computed 
for Mco = 6.47 , 6.75 , ^ = 5/3 , and H/f^ = 0.2 . The 

exact ionopause shape illustrated in these figures is seen to be 
somewhat slimmer than the approximate shape for H/i^ = 0.2 presented 
in figure 6-6. The ionosphere tail in the exact calculation also 
tapers slightly inward toward the axis of symmetry rather than ex- 
tending straight domstream as did the approximate shape. This inward 
tapering just indicates that the solar wind closes in behind the 
ionosphere plasma some distance doxmstream from the planet as might 
be expected since the pressure in the ionosphere decreases with in- 
creasing altitude and Venus is a nonmagnetic planet. Also included 
in figures 11-3 is a plot of Mariner 5‘s trajectory near Venus. The 
plotted trajectory has been obtained by transforming the ASE coor- 
dinates of the trajectory given in figure 11-1 by Bridge £t a^. (1967) 
to a second cartesian coordinate system having its X axis along the 
Sun-Venus line, its y axis pointing in the direction of Venus's 
motion, and its Z axis completing the right-handed system. To take 
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Figure 11-3. Trajectory o£ Mariner 5 past Venus rotated into a plane about the wind -Venus line and superimposed 
on the theoretical location of the ionopause and bow wave as well as various field properties for 
aligned hydromagnetic flow represented by lines of constant-property contours for; (a) Density ratio p/p 
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into account the effect of the aberration angle of 3.4 degrees 
arising from the relative motion of Venus in its orbit about the Sun, 
this new system was then rotated 3.4 degrees about its 2 axis. 
Finally, since our flow solution possesses axial syinmetry, the y and 
Z coordinates were formed into a radial coordinate r = (y^ -5- , 

and these x and r coordinates of the trajectory were then plotted 
in figures 11-3. 

The magnetometer and plasma-probe experimenters have eon- 
cluded from the presence of abrupt and easily recognizable changes in 
their data that Mariner 5 crossed the bow wave at the points of the 
trajectory labeled and (Bridge ^ al., 1967). Point 

is in virtually perfect agreement with the theoretical location of the 
hydromagnetic shock wave although point is somewhat farther up- 

stream from the theoretical results. But figures 11-3 also indicate 
that -Mariner 5 should have crossed the ionopause. These results, how- 
ever, are best interpreted by comparing the observed density, velocity, 
and field strength of the solar wind with the corresponding calculated 
properties. Figure 11-4 shows the time variation of the bulk velo- 
city V , ion number density H , and field intensity |b| observed 
by Mariner 5 as it flew past Venus originally as presented in figure 
8-13 only now enlarged. Superposed on this plot are dashed lines 
representing the flow properties theoretically calculated at the posi- 
tion of the spacecraft’s trajectory illustrated in figures ll-3(a) to 
(c). 

Although conditions in the solar wind were reasonably steady, 
study of figure 11-4 shows that all the data presented display some 
fluctuations just before' time 1 . After a discontinuous increase' at 
time 1 , the agreement between the theoretical and observed value of 
the intensity as they both decrease somewhat up to time 2 is quite 
good. Agreement is also good between the theoretical and observed 
velocities as they likewise decrease slightly to time 2 . The ob- 
served ion density, however, does not rise discontinuous ly through the 
shock as does the theoretical density, but the agreement is nearly 
perfect from then till just before time 3 when the spacecraft passed 
through the calculated ionopause and out of the theoretical flow field. 



km/ sec 








Figure 1,1-4. 


-ISO -120 -60 0 

TIME FROM ENCOUOTER, min 

Plasma and magnetic -fie Id, data measured by Mariner 5 in its encounter with Venus on 
19 October 1967 (Bridge ^ aT. , 1967) . Superposed on these- are dashed lines repre- 
senting the theoretically calculated flow properties that occur along the Mariner 5 
traiectory as illustrated in figures 11-3 (a) to (c) . 
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Meanwhile after time 2 the measured velocities begin to be somev7hat 
lower than their calculated counterparts as they both decrease to' 
just before time 3 when the spacecraft entered the theoretical iono- 
sphere and the computed value falls discontinuously to zero in our 
hydrostatic model of the ionosphere. The observed values also decrease 
sharply but only to a value about one-half the interplanetary level. 

At the same time fluctuations in the field intensity have substantially 
increased after, time 2 , although the average intensity remains at 
about interplanetary level till after time 4 , while the theoretically 
calculated intensity continues to decrease slightly until falling dis- 
continuously to zero just before time 3 at the ionopause calculated 
in our theory for a nonmagnetic planet. After emerging from the 
theoretical ionosphere after time 4 , the spacecraft observed a 
rapid and large increase in field intensity which corresponds fairly 
well with the displayed theoretical values. Just before time 5 the 
theoretical intensity falls discontinuously across the shock wave to 
the interplanetary value while at nearly the same time the observed 
intensity also decreases sharply to the interplanetary value, Simi- ■ 
larly, after the spacecraft re-enters the theoretical flow field 
after time 4 , it observed a sharp increase in density closely cor- 
responding to the large compression in the theoretical density v?hich 
then falls discontinuously through the shock wave to interplanetary 
values just before time 5 . Although during this interval the cal- 
culated solution indicates a discontinuous drop following a larger 
rise in density than that observed, this discrepancy is not serious 
since it involves only one data point which is a 5 -minute average and 
cannot sharply represent such a large compression rapidly followed 
by a discontinuous decrease . Just after time 4 the theoretical 
velocity returns to about three-fourths of the interplanetary value 
and roughly agrees with the measured velocities. At the theoretical 
shock position the calculated velocity jumps discontinuously to its 
full interplanetary value while the observed velocity continuously 
increases up to that level. 

Overall the theoretically calculated flox^ properties agree 
with those observed by Mariner 5 very well indeed. In particular. 
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this good agreement appears to suggest that the upstream crossing of 
the bow wave was encountered >by Mariner about 5 minutes sooner than 
the experimenters had originally indicated with their label 
Support for this belief increases if the fluctuations superposed on 
the interplanetary value of intensity just before time 5 are dis- 
regarded and the time-averaged nature of the plasma data i^s taken into 
account. If so, the observed bow-wave crossings and the calculated 
position of the bow wave are in virtually perfect agreement. These 
theoretical results also indicate that Mariner 5 should have crossed 
the ionopause into the ionosphere where- the plasma velocity and mag- 
netic-field intensity would be zero and the density would be that of 
the ionosphere plasma. However, although the data display decreases 
fit n/rico 3,nd V/Vco that are at least qualitatively similar to 
the theoretical values, there is no certain evidence that Mariner 5 
actually penetrated the ionosphere. In particular, the intensity of 

■ r 

the magnetic field remained of the order of that for interplanetary 
space, and V diminished to only a modest fraction of Vco rather 
than to a value comparable V7ith the speed of the spacecraft relative 
to Venus. 

These results for the more refined hydromagnetic solution 
using the improved calculation for the tapering ionosphere tail give 
further support to the concept, originally suggested in chapter VIII 
as an explanation of a similar discrepancy with the simpler gasdynamic 
model, that a relatively thick boundary layer grows along the iono- 
pause shape betv?een the planetary ionosphere and the solar wind. Be- 
cause the ionopause is essentially a boundary between these tvjo fun- 
damentally different bodies of plasma, it is highly probable that the 
ionopause would provide an increasingly thicker transition between 
the ionosphere and the solar wind with increasing distance from the 
stagnation point at the nose of the ionopause. As such, it would 
spread both into the upper part of the ionosphere and outward into 
the surrounding flow, broadening the transition region in which all 
properties of the plasma including the density, velocity, and magnetic 
field change from their values for the flowing solar plasma to those 
of the planetary ionosphere. Thus at the location of the ionosphere 
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boundary indicated by the present dissipationless hydromagnetic 
theory, the plasma velocity would probably be substantially less than 
indicated by the theory. Since the plasma velocities observed by 
Mariner 5 at the times when it penetrated the theoretical location 
of the ionopause display such a trend, we conclude that just before 
time 3 Mariner 5 entered the boundary layer separating the iono- 

V 

sphere and the flowing plasma, but did not enter the ionosphere proper, 
and finally emerged from it shortly after time 4 . 

One peculiar aspect of these data is the substantial increase 
in fluctuations of the magnetic -fie Id intensity at time 2 . Even 
though there is little change in average field intensity, the experi- 
menters (Bridge et 1967) chose to interpret these fluctuations 

as changes in the Interplanetary medium. Another possible interpre- 
tation is that these fluctuations are generated by disturbances in 
the ionosphere boundary, possibly caused by turbulence in the boundary 
layer. To give some idea as to the location of such boundary dis- 
turbances, we have plotted in figure 11-5 solid lines representing 
characteristic lines of the flow which correspond to 'standing hydro - 
magnetic compression or expansion waves of infinitesimal amplitude. 

Also included in this figure is a dashed line representing the charac- 
teristic line passing through point of the spacecraft's trajec- 

tory. The segment of the ionosphere boundary downstream from the 
intersection point A of this line with the ionopause is then' a 
possible region for the generation of disturbances which propagate 
out to point . We might thus infer that a turbulent boundary 

layer exists downstream of point A . Actually, though, the fluc- 
tuations in the magnetic-field data measured by Mariner 5 persist 
after emerging from the theoretical boundary layer near point , 

which suggests that this boundary layer may be turbulent much farther 
upstream than point A . Such interpretations , however , must be 
regarded as purely speculative at the present time. 

The overall good agreement betoeen the observations made by 
Mariner 5- near Venus and the corresponding calculated properties of • 
our hydromagnetic model offers excellent support for the validity of 
our hydromagnetic theory for the interaction of the solar- wind with 
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nonmagnetic planets, such as Mars and Venus, which have a sufficient 
ionosphere to deflect the solar plasma around the planet and its 
atmosphere. In particular, it has verified the supposition that the 
current sheath which bounds the solar wind away from the planet is 
formed by interaction with the ionosphere rather than with a planetary 
magnetic field as in the case of the Earth and that continuum hydro- 
magnetic theory sufficiently represents this interaction of the solar 
plasma with the ionosphere and the consequent development of a shock 
wave as well as the variation of the average plasma and magnetic -field 
properties throughout the surrounding flow. Furthermore, the compari- 
son of these theoretical results with spacecraft observations also 
demonstrates the validity of treating the subsolar segment of the 
boundary between the solar vrind and the upper planetary ionosphere 
as a tangential discontinuity even though on the remaining segment 
of this theoretical boundary certain differences suggest the presence 
of a thick transition region or boundary layer. 

We therefore conclude that the theory given here is both 
plausible and capable of providing a reasonably accurate description 
of the conditions that prevail at Mars and Venus. 



CHAPTER XII 


COOTEKDIUM 

To provide a basis for the theoretical approach of the major 
topic of this thesis, we have reviewed the data acquired in recent 
years in. the vicinity of the Earth, Moon, Mars, and Venus that have 
disclosed three essentially different types of interaction with the 
solar wind. In addition, the basic properties and idealizations of 
the solar wind itself, which have been predicted theoretically and 
also actually observed in space, were outlined in some detail in 
order to illustrate the basic unity of the entire theoretical descrip- 
tion of the solar wind and its ‘interaction with these planets. 

The fundamental assumption underlying our theoretical analysis 
of the large-scale features of the interaction of the solar wind with 
the Earth, Mars, and Venus is that the average bulk properties of the 
flow can be adequately described by a continuum -flu id theory. Al- 
though this assumption has not been theoretically justified, we re- 
marked that the presence of a weak and irregular magnetic field in 
the incident solar -wind plasma appears to couple the motions of par- 
ticles even in the absence of collisions. More convincingly, v?e 
argued that the real support for the use of the continuum- flu id 
model is provided by the outstanding agreement between results cal- 
culated in this way and those actually measured in space. 

Our fluid model is based on the differential equations of 
magnetohydrodynamics for the steady flow of a perfect gas having 
infinite electrical conductivity as well as both zero viscosity and 
thermal conductivity. Study of these equations reveals five classes 
of discontinuities that can exist in. the flow because of the omission 
of dissipative terms from the differential equations. T^to of these, 
the fast hydromagnetic shock wave and the tangential discontinuity, 
are of concern in our analysis. The first is relevant for relating ■ 
conditions on the two sides of the boT? wave which forms in the solar 
wind upstream of the planets, and the latter is the only .one that has 
properties compatible with those required to describe a boundary 
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surface that deflects the solar wind around the planet. For non- 
magnetic planets which have a highly conducting ionosphere, such as 
Mars and Venus, this boundary is formed by interaction of the solar 
wind with the ionosphere rather than with the geomagnetic field as 
in the case of the Earth. 

With the problem so formulated, there follows an account of 
the manner in which these equations can be simplified to obtain a 
tractable mathematical problem without undue loss of realism or ac- 
curacy of representation of principal features of the interaction. 
First, by using the Newtonian pressure formula the approximate shape 
and location of the ionosphere boundary is calculated for selected 
values of the parameters that characterize the ionosphere just as in 
the analogous problem of the Earth. Following that, the magneto- 
hydrodynamic equations are then simplified to those of gasd 3 m.ainics by 
dropping the magnetic terms on the basis of -the large magnitude of 
the Alfven Mach number. The location of the bow wave and the proper- 
ties of the flow field are then determined for a substantial range of 
values for the parameters describing conditions in the solar wind by 
solving these simpler gasdynamic equations using numerical techniques 
specifically developed for that purpose. It is found that the results 
for a wide range of ionospheric parameters can be brought into close 
correspondence with those for flow past the Earth's magnetosphere by 
application of a simple geometric transformation of the coordinates. 
This correspondence rule enables a substantial body of theoretical' 
results available for solar -wind flow past the magnetosphere to be 
applied to Mars or Venus with only minor modification. Finally, these 
approximate theoretical results are compared with the observations 
made by Mariner 5 as it flew past Venus, and implications of the points 
of agreement and disagreement are discussed, A similar comparison 
with data from Mariners 4, 6, and 7 indicates that Mars has a suffi- 
cient ionosphere for the theory to be applicable and that the observed 
position of the Martian bow wave is exactly predicted by the theory. 

A more refined solution is carried out for the complete mag- 
neto hydro dynamic equations for the special case of parallel velocity 
and magnetic -fie Id vectors. The explicit hydromagnetic equations for 
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aligned flow are derived, and the numerical procedures developed 
originally for gasdynamics are modified so that solutions to these 
equations are able to be computed for a substantial range of values 
for the parameters describing conditions in the solar wind. Study 
of these exact hydromagnetic solutions discloses a strong dependence 
on the content of the magnetic energy relative to the kinetic energy 
in the' solar wind. When this ratio is low, or equivalently when the 
Alfven Mach number is high, this refined solution shows that the bow- 
wave location and flow-field properties are very similar to those of 
the simpler gasd^mamic solution. However, when this ratio is high, 
the bow-wave location and flow properties are strikingly different. 

This unusual phenomena for low-Alfven -Mach -number flow past the Earth 
has recently been confirmed by spacecraft measurements. A detailed 
comparison of these two models is also carried out to determine the 
accuracy and usefulness of the simpler gasdynamic model. 

The calculation of the exact shape of the ionopause is also 
performed by determining as part of the flow solution the points at 
which the tangential discontinuity relations are satisfied exactly. 
Lastly, these refined theoretical results are found to agree better 
with the observations made by Mariner 5 as it flew past Venus than 
do the simpler gasdynamic results. However, certain differences near 
the theoretical location of the ionopause suggest the presence of a 
thick boundary layer . 

Overall, we conclude that the theory presented here is both 
plausible and capable of accurately describing the phenomena associated 
with the interaction of the solar wind with the Earth, Mars, and Venus. 
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